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SOME PHOTONEUTRON THRESHOLDS! 
B. G. Cuipvey,? L. Katz, AND S. KOWALSKI 


ABSTRACT 
The photoneutron thresholds for a number of elements which have 100°, 
abundance (or nearly so) in one stable isotope have been measured. The following 
values in Mev. were found: 


Na (12.47+0.05), Al (12.98+0.08), P (12.50+0.05), V_  (11.16+0.05), 
Mn (10.14+0.05), Co (10.44+0.05), As (10.24+0.05), Y (11.82+0.05), 
Nb (8.86+0.05), Rh (9.46+0.08), I (9.14+0.05), Cs (9.11+0.05), 
La (8.81+0.05), Pr (9.46+0.05), Tb (8.16+0.05), Ho (8.10+0.05), 
Tm (8.00+0.05), Ta (7.66+0.05), Au (7.96+0.07). 


INTRODUCTION 
In the course of a recent investigation on the photoneutron cross sections 
of elements which have 100% abundance (or nearly so) in one stable isotope, 
we also measured their photoneutron thresholds. These measurements were 


particularly desirable since fairly large samples of the rare earths were available 
and because of recent improvements in the accuracy to which such thresholds 


could be determined. 


EXPERIMENTAL TECHNIQUES 

The samples used were in their elemental form except for the rare earths, 
which were oxides. They varied somewhat in weight depending on availability 
but were mostly about 50 g. each and they were over 99% pure in all cases. 

The (y, 2) reaction yield was measured by detecting the emitted neutrons 
in an apparatus previously described (Montalbetti e¢ a/. 1953; Rybka and 
Katz, to be published). This apparatus consisted of four BF; counters embedded 
in a paraffin cylinder with the sample placed in a hole passing through the 
center of the cylinder. About 10% of the emitted neutrons were detected. 

Since only the threshold for the (y, 2) reaction was of interest, no correction 
was applied for photon self-absorption in the samples or for variation of 
detector sensitivity with energy. In any case these corrections were sensibly 
constant over the limited energy regions involved. 
a ‘Manuscript received November 21, 1957. 
_ Contribution from the Department of Physics, University of Saskatchewan, Saskatoon, 
ag UO Institute of Technology, Department of Physics, Pittsburgh, Pa., U.S.A. 
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Thresholds were determined by plotting the neutron yield per X-ray dose 
of irradiation against the peak bremsstrahlung energy, which was varied in 
steps of about 50 kev. The threshold was taken as the energy at which the 
yield curve crossed into the background. This was difficult to determine to 
the desired accuracy in some of the cases, particularly in cobalt, iodine, thulium, 
and gold. Over the limited energy region involved it is possible to write 
(Sher et al. 1951; Birnbaum 1954; Axel and Fox 1956) 


Y ~ (E-E,)" 
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Fic. 1. The photoneutron thresholds in sodium (12.47+0.05 Mev.) and aluminum (12.98 
+0.08 Mev.). The straight lines shown for aluminum were fitted by least squares. The break 
at 13.34 probably represents an average over the 0.23 Mev. and 0.42 Mev. levels in Al?®. 

Fic. 2. The photoneutron thresholds in phosphorus (12.50+0.05 Mev.) and vanadium 
(11.16+0.05 Mev.). 

Fic. 3. The photoneutron thresholds in manganese (10.14+0.05 Mev.) and cobalt (10.44 
+0.05 Mev.). 

Fic. 4. The photoneutron thresholds in arsenic (10.24+0.05 Mev.) and yttrium (11.82 
+0.05 Mev.). 

Fic. 5. The photoneutron thresholds in niobium (8.86+0.05 Mev.) and rhodium (9.46 
+0.08 Mev.). 

Fic. 6. The photoneutron thresholds in iodine (9.14+0.05 Mev.) and cesium (9.11+0.05 
Mev.). 
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where Y is the yield, - is the maximum bremsstrahlung energy, ” is a constant, 
and E, is the desired threshold energy. 

In any one threshold determination we measured Y asa function of £. From 
these data it was possible to evaluate both n and £, in a log—log plot as origin- 
ally used by Katz et a/. (1950) to determine 8-ray ranges and applied by Sher 
et al. (1951) to thresholds. The log—log plots tend to give thresholds which are 
somewhat lower than would be obtained by visual inspection of the data, the 
average difference being 0.02 Mev. with a maximum of 0.06 Mev. (for gold) 
in one case. 

Plots of Y asa function of & are shown in Figs. | to 10. 
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Fic. 7. The photoneutron thresholds of lanthanum (8.81+0.05 Mev.) and praseodymium 
(9.46+0.05 Mev.). 

Fic. 8. The photoneutron thresholds of terbium (8.16+0.05 Mev.) and holmium (8.10 
+0.05 Mev.). 

Fic. 9. The photoneutron threshold of thulium (8.00+0.05 Mev.) and tantalum (7.66 
+0.05 Mev.). 

Fic. 10. The photoneutron threshold of gold (7.96+0.07 Mev.). The threshold of bismuth 
was used to calibrate the betatron energy scale. It is shown here to indicate how precisely it 
could be located. 


BETATRON ENERGY CALIBRATION 

It has been pointed out in a number of publications (Katz et al. 1954; 

Penfold and Spicer 1955) that it is now possible to hold the peak bremsstrah- 
lung energy in a betatron to better than +5 kev. for hours. 

Extensive tests have shown the setting, /, of our energy-controlling helipot 

to be linear in the momentum of the electrons which give rise to the brems- 

strahlung. That is 


(1) T = at+bp = atdV|E(E+ 1.022)|, 
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where a and 0 are constants, p is the electron momentum, and E£ is its kinetic 
energy less moc, which is also the peak bremsstrahlung energy. In calibrating 
the energy scale a and 6 were evaluated by noting the value of J corresponding 
to the photoneutron thresholds of H? at 2.226+0.001 Mev. (Mattauch et al. 
1956), Bi??? at 7.430+0.050 Mev. (Huizenga 1955), and Cu® at 10.826+0.020 
Mev. (Quisenberry et al. 1956). 

The linearity of the energy scale was established by checking it against 
the photoneutron thresholds of Be’, H?, Pb?®’, Bi?®, Cu®, K%9, O'8 and C? 
(IKXatz, unpublished). All measurements were consistent to within the known 
accuracy of their energies with this equation. If we use the other thresholds 
as standards we find the bismuth threshold at 7.40+0.03. This serves as a 
check on the value of 7.480+0.050 accepted by Huizenga and used in the 
present work. Bismuth was adopted as a standard because its neutron yield 
rises rapidly above threshold, Fig. 10, resulting in good reproducibility. 

The position of the deuterium threshold could be established to within 
+10 kev. on repeated measurements, bismuth to about the same accuracy, 
and copper to within +20 kev. Hence the calibration of the energy scale 
in the region of our measurements was known to better than +40 kev. Thresh- 
olds, as will be seen from the graphs, were determined to better than +30 
kev.* Thus we estimate the over-all accuracy to be +50 kev. During the 
course of the experiments, the energy scale was checked periodically. This was 
accomplished with a monitor which was able to detect an energy drift of 
+3 kev. (Katz et al. 1954). 


RESULTS AND DISCUSSION 

The threshold energies found are tabulated in Table I. In most cases several 
determinations were made on each threshold and averaged. While all these 
measurements were completed some 18 months ago, they were repeated 
recently on Na, Al, P, and As. In this case both sets of results are given for 
comparison, with the older value in brackets. In particular, arsenic was 
remeasured to check its threshold, which was not too well established origin- 
ally. The new value is the average of two separate determinations and falls 
somewhat higher than the original one. The aluminum measurement was 
repeated to see if the 7 = 1 level at 0.23 Mev., J = 0+ in AP® could be 
located. The average of five separate runs is shown in Fig. 1. A break was 
found at 13.34 Mev., but probably represents an average over the 0.23 Mev. 
and 0.42 Mev. levels. Our statistics are just not good enough to resolve these 
levels. 

To compare our results to known nuclear masses in combination with 
B-decay energies it will be useful to use the following notation. Irradiation 
of isotope X“ with photons gives 


(2) X4+hv — n+ X4-1, 
The residual nucleus X“~' then decays by either 8~ emission 
(3) x 4-1 — B-+Q0-+Y4- 


*Except Al, Rh, and Au, where we estimate the determination to be within +70 kev. 
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MEASURED PHOTONEUTRON THRESHOLDS 


Measured 


Other Q values, 








Reaction Q value, Mev. Mev. Method Reference 
Na*4(y, n)Na® 12.47+0.05 
(12. 4640.08) 
12.05 +0.20 Threshold Sher et? al. (1951) 
12.42 +0.05 Massdata  Wapstra (1955) 
12.417+0.014 Review Mattauch et al. (1956) 
Al7(y, »)AL® 12.98+0.08 
(12.93+0.08) 
14.4 +0.3 lhreshold Becker et al. (1947) 
14.0 +0.1 lhreshold McElhinney et al. (1949) 
12.75 +0.20 Threshold — Sher et al. (1951) 
13.06 +0.06 Massdata Wapstra (1955) 
9 Q* value Laubitz (1955) 
12.98 +0.11 \ Mass data Mattauch et al. (1956) 
P31(, n) P30 12. 5040.05 
(12.45+0.08) 
12.35 +0.2 Threshold McElhinney et al. (1949) 
12.4 +0.2 Threshold Katz and Penfold (1951) 
12.05 +0.20 Threshold Sher et al. (1951) 
12.40 +0.08 Massdata  Wapstra (1955) 
12.391+0.030 Review Mattauch et al. (1956) 
V5l(+, n)V% 11.16+0.05 
11.15 +0.20 Threshold Sher et al. (1951) 
11.24 +0.23 Massdata Wapstra (1955) 
11.05 +0.12 Mass data Duckworth (unpublished) 
Mn®(y, 2)Mn*4 = =10.14+0.05 
10.15 +0.20 Threshold Hanson et al. (1949) 
10.00 +0.20 Threshold = Sher et a/. (1951) 
10.14 +0.26 Massdata Wapstra (1955) 
~ { Massdata Duckworth (unpublished) 
10.28 +0.07 § Q* value Jung and Pool (1956) 
{ Mass data) Duckworth (unpublished) 
10.09 +0. 10 \ O- value Way et al. (1955) 
Coy, 2)Co 10.44+0.05 
10.25 +0.20 Threshold — Sher et al. (1951) 
{ Mass data Duckworth (unpublished) 
10.491 +0.012 \ Q* value Way et al. (1955) 
10.490+0.012 Mass data Quisenberry et al. (1956) 
As75(y, n)As74 10.24+0.08 
(10.14+0. 10) 
10.18 +0.19 Massdata Wapstra (1955) 
10.3 +0.2 Threshold Ogle et al. (1950) 
10.10 +0.20 Threshold — Sher et a/. (1951) 
-, {Mass data Duckworth (unpublished) 
( ‘ 
9.92 +0.07 \ Q* value Way et al. (1955) 
~»{ Mass data Duckworth (unpublished) 
i * 
ae Gen io value Way et al. (1955) 
Y*%y, m)¥88 11. 82-40.05 
ie - { Mass data Duckworth (unpublished) 
11.76 £0.15 Q* value Way et al. (1955) 
11.80 +0.33 Massdata Wapstra (1955) 
Nb%(4, 2) Nb”? 8. 86+0.05 
8.70 +0.20 Threshold — Sher et al. (1951) 
8.69 +0.40 Mass data Wapstra (1955) 
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Measured 


Q value, Mev. 


9.46+0 


9.14+0 


9.112:0. 


8.8140 


».46+0 


8.1640. 


8.10+0 
8 .00+0 
7.6640 


7.9640 


O8 


05 


05 


O07 


TABLE I (concluded) 


Other Q values, 


Mev. 


+0.20 
+0.30 


9.35 
9.45 


9.41 +0.34 


+0.2 
+0 .2 
+0.20 


+0.12 


9.45 
9.10 


9.46 


9.12 


+0.20 
+0.19 


9.05 
9.03 


+0.20 
+0.20 


+0.19 


8.80 
8.73 
9.07 


+0.10 
+0.3 
+£0.06 


4 +0.2 


+0. 20 


8.00 +0.15 
+0.1 
+0.20 


+0.40 


7.90 
7.78 


£0.17 | Q- value 


Method 


Threshold 


{ Mass data 


Q> value 


Mass data 
Q* value 


Threshold 
Threshold 
Threshold 
Mass data 
Qt value 

Mass data 


Threshold 
Mass data 


Threshold 
Mass data 


{ Mass data 


Q- value 


Threshold 
Threshold 
Mass data 


Threshold 
Chreshold 
Vhreshold 


Threshold 
Threshold 
Threshold 
Mass data 





MEASURED PHOTONEUTRON THRESHOLDS (concluded) 


Reference 


Sher et al. (1951) 
Duckworth (unpublished ) 
Kochendorfer and 

Farmer (1954) 
Duckworth (unpublished ) 
Kochendorfer and 

Farmer (1954) 


McElhinney et a/. (1949) 
Ogle et al. (1950) 

Sher et al. (1951) 
Duckworth (unpublished ) 
Koerts et al. (1955) 
Duckworth (unpublished ) 
Koerts et al. (1955) 


Sher et al. (1951) 
Johnson and Nier (1957) 


Sher et al. (1951) 
Johnson and Nier (1957) 
Duckworth (unpublished ) 
Glover and Walt (1957) 


Hanson et al. (1949) 
Ogle et al. (1950) 
Johnson and Nier (1957) 


McElhinney et al. (1949) 
Johns et al. (1950) 
Sher et al. (1951) 


Hanson et al. (1949) 
Parsons and Collie (1950) 
Sher et al. (1951) 

Hay (1955) 


*Errors given are from mass values only, those of Q* and Q~ values are not known. 


or 8+ emission and/or orbital capture 


(4) 


or 


X4-1 — B+4+Qt+Z4— 


> Qt4+Z4-, 


Where possible the neutron binding energies were calculated from atomic 
masses or from the masses and the Q- or Q* energies whose sources are quoted 
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in the last column of Table I. Previously measured thresholds are also quoted 
for comparison. 

Axel and Fox (1956) have discussed at length the validity of various extra- 
polation methods for determining photoneutron thresholds. In particular they 
emphasize that the equation Y ~ (E—E,)" can only be used with some 
assurance over a very limited energy region. They also point out that Birn- 
baum’s careful work shows that a precise apparent threshold can be defined 
by this equation but question whether it corresponds to the physically real 
one. 

An indication of the reliability of our method for finding the true photo- 
neutron threshold can be obtained by comparing our results to those cases 
in which the thresholds are known to at least the same accuracy from mass 
data. In Table II we list these cases and compare them. No systematic differ- 
ence is apparent and, except for praseodymium, the differences fall within 
their combined error. 


TABLE Il 


COMPARISON OF THRESHOLDS FROM MASS DATA AND FROM PHOTONEUTRON REACTIONS 


Photoneutron Mass data Difference, 


Reaction threshold, Mev. threshold, Mev. Mev. 
Na®(y¥, 2)Na® 12.47+0.05 12.417+0.014 —0.05+0.05 
Al?7(4, 2 )AI?6 12.96+0.06 13.03 +0.06 +0.07+0.08 
P31( +, 2) P29 12.48+0.05 12.29 +0.08 —0.09+0.06 
Coy, n)Co%* 10.44+0.05 10.49 +0.01 +0.05+0.05 
Pri4(y, n)Pr? 9 46+0 05 9.30 +0.06 —0.16+0.08 
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PHOTONEUTRON THRESHOLD AND 
CROSS SECTION FOR ;,Lu'! 


H. J. Kinc? ann L. Katz 


ABSTRACT 


The neutron yield resulting from photoneutron reactions in Lu!”® has been 
measured as a function of peak bremsstrahlung energy up to 23 Mev. The 
threshold energy for this reaction was found to be 7.77+0.05 Mev. The giant 
resonance cross section has a peak value of 225 millibarns at 16 Mev., a half- 
width of 8.4 Mev., and an integrated cross section to 23 Mev. of 1.9 Mev-barns. 


INTRODUCTION 


The determination of the lutetium!”® photoneutron threshold energy and 
cross section from threshold through the giant resonance was undertaken as 
part of a larger program in this laboratory (to be published). The purpose 
was to investigate the (y, 7) activity of all nuclides in which one isotope 
is very abundant (over 97%), particularly the rare earths, in an attempt to 
discover any regularities which may exist in their photonuclear properties. 


EXPERIMENTAL TECHNIQUE 

Except for minor improvements in circuitry the apparatus used to detect 
neutrons resulting from the (y, 2) reaction has already been described (Montal- 
betti et al. 1953). It consists of four BF; counters in a paraffin moderator, 
coupled with a gating circuit and scaler. The sample was placed at the center 
of the paraffin tank, and the X-ray beam from the betatron was collimated 
in such a manner that no X-rays struck the paraffin. The neutron yield was 
measured as the peak bremsstrahlung energy was changed in 400-kev. steps 
to define the shape of the (y, 2) vield curve, and by 20-kev. to 40-kev. steps 
to define the threshold energy. The peak bremsstrahlung energy could be 
controlled to better than +5 kev. for periods of several hours. The gating 
circuit was triggered by the appearance of the X-ray pulse (180 p.p.s.) and 
defined a neutron counting period which started 20 usec. after each X-ray 
pulse and lasted for 250 usec. For threshold measurements all four BF; 
counters were used to increase the neutron detection sensitivity (~10°%); 
however, to measure the neutron yield from 8 to 23 Mev. only one counter 
tube was used. This tube was positioned by trial and error in the paraffin 
moderator until the response was independent of the initial neutron energy. 
The method used to achieve this was previously outlined (Montalbetti et 
al. 1953). 

The relative dose was measured with a Baldwin Farmer* dosemeter and 
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the absolute yield determined by comparing the observed activity to that 
of a copper disk irradiated in the same geometry at 20 Mev. The absolute 
yield for Cu was taken as 2.16108 n./100 r./mole, a value 10% less than 
that given by Montalbetti et a/. (1953) because of a systematic error in the 
Victoreen dosemeter used at that time. This gave 1.93109 n./mole/100 r. 
as the neutron yield from Lu at 20 Mev. 


ANALYSIS 


1. Threshold 
The results are plotted in Fig. 1, which shows a constant background and 


a sharp threshold at 7.77+0.05 Mev. 


COUNTS 
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INTEGRATOR SETTING 


Fic. 1. 


2. Cross Section 

The activation curve, Fig. 2, corrected for background due to cosmic rays 
and due to the aluminum container (less than 2% correction), was normalized 
to absolute disintegration rate, smoothed, and analyzed by the Modified Total 
Spectrum Method (Leiss and Penford, to be published). The resulting cross 
section calculated at }-Mev. intervals is shown in Fig. 3. 


RESULTS 
The values found were: 


1. Threshold 7.77+0.05 Mev. 
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2. Cross section 
Area to 23 Mev. 1.9 Mev-barns 


















Peak value 225 millibarns at 16 Mev. 
Half-width 8.4 Mev. 
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The threshold results are in agreement with a value of 8.06 Mev. calculated 
from Cameron's (1957) modified semiempirical mass formula. No experi- 
mental mass data or other measuremenis are available for comparison. 
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The cross section is of the same general shape as the others measured in 
this mass region and is in agreement with a trend to wider half-widths between 
magic neutron numbers reported by Nathans and Halpern (1954) and apparent 
from the afore-mentioned survey in this laboratory. 
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ANGULAR ENERGY FLUX OF SECONDARY GAMMA-RAYS 
IN MATTER: SMALL-ANGLE SCATTERING 
FROM A POINT ISOTROPIC SOURCE! 


W. R. Dixon 


ABSTRACT 


The distribution in energy and angle of secondary gamma-rays in matter is 
studied, with particular reference to small angles and to energies near the 
incident energy. The ‘‘angular energy flux’’ arising from a point isotropic source 
in matter is calculated, using small-angle approximations, and compared with 
experimental data obtained for a Cs"? source in concrete. It is concluded that 
for small angles and for energies near the incident energy, only a few orders 
of scattering contribute to the angular energy flux out to considerable thick- 
nesses, whereas for large angles and low energies many orders of scattering 
contribute even for small thicknesses. 


I. INTRODUCTION 

The passage of gamma-rays through matter is accompanied by the pro- 
duction of secondary gamma-rays through the processes of Compton scattering, 
the annihilation of positrons, and the emission of characteristic X-rays. The 
general problem of gamma-ray transmission in a medium is to take into 
account in detail the secondary gamma-rays, 1.e., to determine the distri- 
bution function J(r;@, £) such that at the point r in the medium J(r; w, £) 
@wdE gives the energy (in the form of photons) which crosses unit area 
normal to w with photon directions lying in the element of solid angle d*w 
about w and photon energies lying in dE about E. (Here @ is a unit vector.) 
Following Goldstein and Wilkins (1954), we shall refer to J as the angular 
energy flux. We shall also use the same name for the function 

I(r;@,\) = I(t; @, E) dE/dx, 
where \ is the wavelength measured in units of the Compton wavelength. 
(A = 0.511/E, where £ is in Mev.) 

Detailed calculations of the energy flux (defined as the integral of the 
angular energy flux over all angles) have been reported by Goldstein and 
Wilkins (1954) for infinite homogeneous media. These calculations are based 
on the moments method of Spencer and Fano (1951). It is difficult to obtain 
the angular distribution as well by this method, and complete numerical 
calculations of the angular energy flux are not available. However, the feature 
of greatest interest in the angular energy flux is perhaps the rapid fall-off in 
the intensity with increasing angle and decreasing energy for small angles 
and for energies near the incident energy. This feature can be investigated 
theoretically by the methods of Foldy (1951) and of Ogievetskii (1956); such 
calculations for a point isotropic source are outlined in Section II. 
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In order to check these calculations, measurements of the angular energy 
flux arising from a Cs'*? source buried in concrete have been carried out at 
small angles. This work, described in Sections III and IV, is essentially a 
continuation of a study initiated by the author in 1952. In the intervening 
years three similar experiments with Co® sources have been reported, those 
of Whyte (1955) carried out in this laboratory, of Scofield (1957), and of 
Peelle et al. (1957). Other experimental work has been concerned with the 
measurement of the energy flux, and with the determination of build-up 
factors, which are defined as integrals taken over the energy flux with suitable 
weighting factors (Hayward 1952; Weiss and Bernstein 1953; Theus ef al. 
1955; Garrett and Whyte 1954). These integrated quantities in general have 
shown good agreement with calculations based on the Spencer-Fano method. 


II. CALCULATION OF THE ANGULAR ENERGY FLUX 
In the energy region in which only Compton scattering is important in the 
production of secondary gamma-rays, the angular energy flux in an infinite 
homogeneous isotropic medium is determined by the following integrodiffer- 
ential equation: 


(1) [w.grad + u(A)]Z(r;@, 4) = S(r; 0, A) 
2 od 
+f dw! \ dN I(r; 0’, \’)K(N, A, w-@’). 
4a WJ ry 


Here u(A) is the narrow-beam absorption coefficient, S(r;, A) is a source 
term giving the intensity emitted per unit solid angle per unit wavelength 
at and J respectively, and A(X’, A, w.w’)d’wdd is derived from the Klein 
Nishina cross section. It gives the relative intensity of photons at exactly w’ 
and \’ which are scattered into d’@ at and into dd at X: 
(2) K(X’, A, @.@')@eadd = (nre?/2)(N’/d)3[N/A+A/N — 1+ (@.0’)?] 
X 6(A—N —1+0.0’)dwdd. 

In equation (2), 2 is the number of electrons per unit volume, 7o is the classical 
electron radius, and the delta function expresses the Compton relation between 
the change in wavelength and the angle of scattering. 

We consider the diffusion equation (1) with the source term 
(3) S(r; @, AN) = (Eo 27)6(1 —@ .1)6(A — Ao) d(r) dir’, 
where frp is a unit vector in the direction of r. This represents a point source 
located at r = 0, emitting one photon per second of energy /». The primary 
flux is given by 
(4) I (rs0, ) = (Eo/27)8(1—w.18o)6(A— Ane” /4ar’. 
The angular energy flux of gamma radiation which has been scattered only 
once is given by 


(5) I (rs, ) = far’ f d’e’ fan (r':0', N K(X’, A, @-@’) 


exp(— u|r — r’ re _ w.(r —r’)| 
Qalr — rr’)? | r—r’ 
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This reduces to 


Eye**” rK(do, 4) SS ( sin :) 
. ()D 5 2 NEOy Ee oe _ sit 
(6) I'(tie, A) = 4er sin@sina > ye I sina 


sin (a — 6) 


— pr “-u(a — 6), 


sina / 
where cos@ = w.f, cosa = 1—(A—Xo), A (Ao, A) is given by the integral of 
equation (2) with respect to (@.@’), and u(x) is the unit step function: 
u(x) = 1 for x > 0 and O-for x < 6. 

Equation (5) represents a general relation between scattering of orders n 
and (n—1), if J and J® are replaced by J” and J“~-" respectively. It can 
be seen, however, that the resulting expression for even the second-order 
scattering is very complicated indeed, so that other methods of solution are 
desirable. In problems where only the energy flux is required, the method 
devised by Spencer and Fano (1951) is particularly useful. In this method 
the distribution in w is expanded in Legendre polynomials, and the spatial 
distribution in a system of polynomials which is chosen to give rapid con- 
vergence; the energy flux is just the coefficient of the zero-order Legendre 
polynomial in this expansion. However, the method is not useful without 
modification for the calculation of the angular energy flux at small angles 
and at energies near the incident energy, because many terms are required 
to represent the sharp forward peak. Nevertheless, some sample calculations 
of the angular energy flux (at energies considerably below the incident energy) 
have been given by Spencer and Stinson (1952) and a simplification in the case 
of a plane isotropic source has been exploited by Berger (1955). 

On the other hand, approximate solutions have been given by Foldy (1951) 
and by Ogievetskii (1956), which are suitable for small angles and for energies 
near the incident energy. These authors have considered only the plane 
monodirectional source with normal emission. We outline below the similar 
treatment of the point isotropic source, which turns out to have a somewhat 
more complicated solution. 

We first observe that if equation (2) can be approximated by 


(7) K(n’, d, @.0’) = (a/2m) {f(r')/f()} 8(A—0’ —1+0.0’), 


where a = 2rnr,? and f(A) is any function of the wavelength, then the diffusion 


equation for 
(8) g(r; @, A) = {f(A)/f(Ao)} T(t; @, A) 


is simply 
r 
(9) [w.grad + w(A)]g(r;@, A) = (a/2m) { d’w! dy’ g(r; 0’, d’) 
4r @ ro 


X 6(A—A’—1+ 0.0’), 


where the source term (3) has been omitted. The approximation f(A) = ° 
has been used by Foldy (1951) (with 8 = 2), by Ogievetskii (1956) (with 
8 = 2.8), and will be used below (with 6 = 3.0). 
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In order to render the problem mathematically tractable, we make the 
small-angle approximations 


w.grad = cos0d/dr—sin 0(1/r)d/00 = 0/dr—(6/r)d/00 
and w.o = 1—6/2—6/2+06’cos(¢—¢’). 


Furthermore we assume that yu is constant. A plausible argument can be made 
for this assumption, particularly if the incident gamma-ray energy lies close 
to that value for which » has a minimum (Foldy 1951). It is actually found 
that this approximation is useful for moderate depths of penetrations, but 
not in the asymptotic limit of very great depths. This point will be discussed 
later. 

With the above approximations the diffusion equation (9) assumes the 
form 


ad ®a 
(10) dg/dr—(6/r)dg/00+pog = (a/27) | dn’ \ 6'd6’ { 
J ro Ju Jo 


x g(r; 0’, ’)6{A—A’—0°/2—0'"/2+-00'cos(g— $’)}. 


We now make a Fourier—Bessel transformation in @ and a Laplace trans- 
formation in X. (It is well known that the Fourier—Bessel transformation is 
equivalent to a Legendre expansion for small angles.) Defining the transform 
A by 


(11) = EyA(r:p, s) = 42° “iT grr en J 2x0Jo(on)e(r+6 ddd, 
it may be shown that A satisfies the equation 

(12) OA /Or+(p/r)dA/dp = (a/s).exp(—p?/2s) .A, 

for which the solution is 

(13) A303) = exp) (ar/s) J exp(—p'x" 25a : 


The inverse transformation is 


(14) 4ar°e”°’g(r; 8, X) 


a fas. -. ,8(A—Ao) J pdpJo(p@) 


Qn) ie 0 


x > (ar/s)"(1, nn f exp (— px /2s)dx ( 


n= 


= (Eo/27) | barge’ 2)+ar{2(A—Av)O°} 2u(A—Ap—O"/2) 


oo n al al 
os ae J dry... | dx, Ry? (X= o—0"/2Rn*)* 
0 e 


“a n!(n — 2)! 0 


X u(A—Ao—O" 2k.) | . 


where R,? = xP-+xe+ .. . 4-2”. 
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The corresponding solution for the plane monodirectional source has been 
given by Foldy (1951). In our notation it is 


(15) e”°?g(z; 0, ) = (Eo 2m] 80-ra6/2) bas 8(A—Ao—0"/2) 


' 0 
+ », ae a - (X= Ay—6/2n)" (A= 96/2) | 
Successive terms in (14) and (15) represent successive orders of scatterings. 

The solution for a plane monodirectional source for non-constant absorption 
coefficient has been given by Ogievetskii (1956). 

The spectra of the first, second, and third order scattering terms for a 
point isotropic source in concrete are shown in Fig. | for the following con- 
ditions: ar = 10, Xo = 0.772 (Cs!37 source), and for angles @ = 4°, 10°, 20°, 
30°. For the first-order scattering, both the exact term (6) (with wor = 5.14) 
and the approximate term from (14) are shown. The intensity distribution 
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Fic. 1. Calculated values of the angular energy flux for a point isotropic Cs’? source in 
concrete. ar = 10, wor = 5.14, Ao = 0.772. The ordinates have units of Mev. per steradian 
per Mev. interval per second, and are normalized to a source emitting one photon per second. 
(la) Approximate first scatter. (1b) Exact first scatter. (2) Approximate second scatter. 
(3) Approximate third scatter. 
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in & rather than in \ is shown, and we have taken f(A) = A*. The normalization 
is for a source which emits one photon per second. 

For comparison the first five orders of scattering for the plane mono- 
directional case are shown in Fig. 2. The angular dependence here is much 
less than for a point isotropic source, the curves for 4° and 10° being nearly 


@Po-2 (256°) 
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Fic. 2. Calculated values (approximate) of the angular energy flux for a plane mono- 
directional Cs'*7 source. az = 10, Xo = 0.772. The ordinates have units of Mev. per cm.? per 
steradian per Mey. interval per second, and are normalized to a source emitting one photon 
per cm.2 per second. The numbers attached to the curves indicate the order of scattering. 


indistinguishable. The first scatter for a plane monodirectional source is a 
delta function, which is represented in Fig. 2 by a rectangle of the appropriate 
area. The normalization is for a source which emits one photon per cm.* 
per second. 

The asymptotic form of (14) is 


(16) bar” eX" g(r; 0, X) ~ Eo(12a°r’/ me") (246) * exp{¢—36°/2(A—o)}, 
where °° far(\—Xv), while the asymptotic form of (15) is 

(17) e°? 9 (2:0, ¥) ~ Eo(da’2"/me°)(29t) * exp{e—0°/2(A— Do) }, 

where ¢° faz(\—Ay). These forms can be derived from the leading terms 


in the expansions of (14) and (15) in Laguerre polynomials. The expressions 
(16) and (17) do not appear to correspond to physical reality however, because 
they do not show the shape of the distribution in energy and angle to become 
independent of the depth of penetration. This independence was first estab- 
lished for the energy flux by Bethe, Fano, and Karr (1949) by means of a 
calculation involving the approximation that all scattered gamma-rays 
travelled in the forward direction, and verified by Fano, Hurwitz, and Spencer 
(1950), who took into account the effect of small angular deflections. This 
deficiency in (16) and (17) probably can be attributed to the assumption 
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that the absorption coefficient is constant, since there is then no means in 
the calculation of suppressing the pile-up of lower-energy scattered com- 
ponents, and hence no possibility of ever attaining an equilibrium distribution. 


Ill. EXPERIMENTAL ARRANGEMENT 

Whyte (1955) has given an excellent account of the experimental arrange- 
ment and the measurements carried out in this laboratory up to 1954, so it 
will not be necessary to describe the setup in detail. Two major improvements 
in technique have since been introduced: the use of a “total absorption” 
scintillation spectrometer in place of the two-crystal Compton spectrometer, 
and the use of a 20-channel kicksorter in place of the photographic recording 
system. The high efficiency of the present detector makes possible the use of 
weaker, and hence smaller sources, and therefore makes it possible to measure 
the scattered flux at smaller angles than could be used previously without 
interference from the primary radiation. 

In the present experiment two Cs'*7 sources were used of strengths about 
10 curies and 1 curie, and in containers of lengths about 12 in. and diameters 
about in. The sources were imbedded in a concrete barrier composed of 
blocks of dimensions 12 in. by 8 in. by 4 in. and of density 2.37 g./cm.’ 

The gamma-ray spectrometer consists of a large Nal(Tl) crystal, a photo- 
multiplier, and associated electronic circuitry. The Nal crystal, packaged by 
Harshaw Chemical Company, is 4 in. long and 4 in. in diameter, with a well 
of diameter 2 in. and depth Lin. into which a collimated gamma-ray beam 
can be directed. The well was obtained for another purpose, but has the inci- 
dental effect of reducing the escape of backscattered photons from the crystal, 
particularly at low gamma-ray energies. The crystal is mounted on a 5 in. 
diameter Du Mont phototube type 6364. A regulated high-voltage supply 
(Atomic Model 312) is used to provide a positive high voltage of 1100 volts. 

The Nal erystat and photomultiplier are shielded by at least two inches 
of lead on all sides except the rear, where one inch of lead is used. 

An ideal spectrometer would give a unique response for monoenergetic 
gamma-rays, with a small spread due to statistical or other causes. In a 
scintillation spectrometer, a continuous pulse-height distribution arises from 
Compton interactions in the crystal. In large Nal crystals, however, some 
of the scattered photons are absorbed by the photoelectric process, so that 
the corresponding pulses which would have appeared in the Compton distri- 
bution actually appear in the photopeak. Experimental response curves tor 
Co® (1.33 and 1.17 Mev.), Cs"? (0.662 Mev.), and Hg?® (0.279 Mev.) are 
shown in Fig. 3. The photofractions, defined as the areas under the photopeaks 
relative to the total areas, are found to be 0.54, 0.71, and 0.89 for Co®, Cs!87, 
and Hy? respectively, in good agreement with the values 0.54, 0.72, and 
0.96 interpolated from calculations of Berger and Doggett (1956) for a crystal 
of length 8.0 in. and radius 2.0 in. The values from Berger and Doggett are 
plotted in Fig. 4. The corresponding phototractions for a very small crystal 


are just the ratios of the photoelectric to the total cross sections, namely 
0.04, 0.10, and 0.48 respectively for the three sources cited above. 
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Fic. 3. Response of the “‘total-absorption” spectrometer to Co® (1.33 and 1.17 Mev.), 
Cs37 (0.662 Mev.), and Hg? (0.279 Mev.) sources. The small peaks (a) and (6) are due to 
backscattered gamma-rays, while (c) is the K X-ray of TI?%. 
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Fic. 4. Calculated efficiency Y(E) and photofraction p(E) for a Nal crystal of length 
3.0 in. and radius 2.0in. The values of p(E) are interpolated from values given by Berger 
and Doggett (1956). 


The efficiency of the crystal for a collimated beam of high-energy gamma- 

rays is simply 

Y(E) = 1—exp{—z(E)L}, 
where L is the length of the crystal, and u(£) is the narrow-beam absorption 
coefficient. For low energies it is also necessary to take into account absorption 
in the aluminum case and MgO reflector. The calculated efficiencies Y, thus 
modified at the low-energy end, are shown in Fig. 4. 

The “resolution”, defined as the width of the photopeak at half the peak 
value divided by the pulse height for which the peak occurs, is found to be 
13% for 0.662 Mev., and 16% for 0.279 Mev. The linearity of the spectro- 
meter can be seen from Fig. 3 to be adequate. The effects of electron escape 
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and A X-ray escape are not large enough to warrant discussion. The back- 
ground counting rate was found to be about 4000 counts per minute, in 
agreement with a value reported by Foote and Koch (1954); the background 
spectrum did not show any large peaks. The large background rate reduces 
the usefulness of this type of spectrometer, especially with extreme collimation 
of the incident gamma-rays. Because of the dependence of the pulse height 
on the temperature of the crystal and photomultiplier it was found necessary 
to calibrate the spectrometer at least once a day with a Cs!*7 source. Drifts 
due to the electronic equipment were negligible compared with this temperature 
effect. 
IV. RESULTS 

The angular energy flux has been measured for a Cs!*7 source buried in 
concrete to depths of 4in., 8in., 12in., and 16 in. Spectra at angles of 3°, 
{°, 6°, 8°, and 10° for the latter three thicknesses are shown in Fig. 5, together 
with the experimental arrangement used for these particular measurements. 
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Fic. 5. Experimental distributions for a Cs’ source buried to various depths in concrete. 
lhe experimental arrangement is also shown. 
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The collimator was 9 in. in length, and had a hole of diameter § in. With this 
collimator the angular energy flux could be measured at angles as small as 
3° for concrete thicknesses of 8 in. or more, without interference from the 
primary radiation. For a concrete thickness of 4 in., the detector probably 
could receive primary radiation out to angles of about 10°. For this reason 
the results for 4 in. are not shown, although for 6°, 8°, and 10°, the distributions 
have essentially the same shapes as those shown in Fig. 5. 

The sets of curves shown in Fig. 5 are quite similar for the different concrete 
thicknesses, showing either that the shape changes very slowly with thickness, 
or that an equilibrium condition has already been attained for a thickness of 
8 in. Whyte (1955) and Scofield (1957) have also observed that the distribution 
shape changes slowly, if at all, with thickness. 

A comparison of the experimental data with the calculated distributions 
has been attempted by calculating the response of the spectrometer to the 
theoretical angular energy flux, rather than by attempting to unfold the 
experimental curves to obtain the actual spectra. The calculated response 
of the spectrometer (counts per minute multiplied by pulse height) was taken 


to be 


2Eo 
KY ).= j dE I(E)(V/E)Y(E)p(E) (22kE)  exp{ —(V—E)?/2kE} 


ekG 
i‘ | dE I(E)(V/E)Y(E){1—p(E)} j dE! H(E, E’)(24kE’) 


X exp} —(V—F’)* /2kE'}, 


where the first term represents the photoelectric response, and the second 
term the Compton response. Here I’ is the pulse height in energy units, / 
the photon energy, V(/) the efficiency, p(/) the photofraction, k the Gaussian 
parameter, and //(£, £’) the relative probability that a photon of energy 
gives rise to a Compton electron of energy /’. The upper limit of the integral 
over Lt’ is (2E?/moc?)/(1+2E/myc*). In the calculations we have taken H(F, E’) 
to be uniform in #’, in accordance with the experimental results of Fig. 3 
and with the calculations of Berger and Doggett (1956). The value of k 
chosen was k = 0.002 Mev., which corresponds to a resolution of 13% for 
0.662 Mev. The values of /(£) used in this calculation are those shown in 
Fig. 1 for ar = 10, Ao = 0.772. (The value ar = 10 actually corresponds to 
a thickness of 11.1 in. of concrete, which was the distance to the center of 
the source for a nominal thickness of 12 in.) The calculated response to the 
first scatter (exact) is shown in Fig. 6. The dashed curves in Fig. 6 represent 
the response to the sum of the second and third orders of scattering. In order 
to extend the energy range beyond that for which the calculations could be 
expected to be valid, the intensities of the second and third orders of scattering 
(as shown in Fig. 1) have been multiplied by the ratio of the value of the 
first scatter calculated exactly to the value calculated approximately (Fig. 1). 
This ratio is unity for / = 0.662 Mev., and is approximately 0.3 for 4 = 0.450. 

In Fig. 7 the solid curves are the calculated response curves discussed 
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above. The experimental points were obtained for a concrete thickness of 
12in. and with a collimator hole of length 10 in. and diameter } in. The 
counting rates were multiplied by the energy corresponding to the channel, 
and the whole set of experimental curves was then normalized to the set 
of calculated response curves. 


30 
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Fic. 6. Calculated response of the spectrometer to the angular energy flux shown in Fig. 1. 
The ordinates give relative values expected for the product of counting rate per pulse-height 
channel times pulse height. The solid lines show the response to the first scatter (exact), while 
the dashed lines show the response to the sum of the second and third scatters modified as 
explained in the text. 
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Fic. 7. Comparison of experimental distributions for a Cs"? source buried in 12 in. of 
concrete with calculated response curves of the spectrometer. The solid lines give the sum 
of the responses to the first, second, and third scatters as shown in Fig. 6. The experimental 
distributions have been normalized as a set to the calculated response curves. 
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It can be seen from Fig. 7 that only a few orders of scattering suffice to 
account for the observed angular energy flux for angles up to 20° and for 
energies near the incident energy. Although a detailed comparison has been 
made only for a concrete thickness of 12 in., it is clear from the slow variation 
with thickness (as shown in Fig. 5) that this conclusion holds over a wide 
range of thicknesses. For large angles and low energies it is of course necessary 
to take into account many orders of scattering even for very small thicknesses. 

It is concluded that the analytic calculation of the angular energy flux 
outlined in Section II has validity out to thicknesses which are commonly 
met in shielding problems, for small angles and for energies near the incident 
energy. Since the calculation is not valid for large angles, it does not appear 
possible to obtain a useful estimate of the energy flux by integrating over all 
angles. It has also been noted in Section II that the calculation does not 
have the correct asymptotic behavior, a defect which presumably could be 
remedied if it were possible to carry through the calculation using a non- 
constant absorption coefficient, as has in fact been done by Ogievetskii (1956) 
for the plane monodirectional case. Keeping these limitations in mind, there 
is still a wide range of conditions for which such calculations should prove 
useful when a detailed knowledge of the angular energy flux from a point 


isotropic source is required. 
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STUDIES OF SOME INFRARED FRAUNHOFER LINES! 


G. GODOLI 


ABSTRACT 


The apparatus function of the Dominion Observatory solar infrared spectro- 
meter was determined by laboratory observations of the CH, band at 1.7 u 
and checked by observations of the telluric CO band at 1.6 4. The equivalent 
widths and the profiles of four Fe solar lines in the 1.6 4 region were determined 
with the same instrument at seven different points on the solar radius. The 
equivalent widths agree better with those determined by Pagel than with those 
determined by Mohler. The profiles, corrected for the apparatus function, 
suggest a very high ratio y/ye1 and a dependence of the kinetic temperature 
on the line strength. 


A. INTRODUCTION 

The necessity for accurate profile studies of solar spectral lines in the 1.6 4 
region, where the continuous absorption coefficient of the solar atmosphere 
reaches a minimum and where the excitation potential of the lines is very 
high, has been pointed out by several authors. 

Pagel (1955a, 6) suggested the study of the profile of the strong Mg line 
\ 17,108.8 A and estimated that even for such a strong line a resolving power 
of 10° would be necessary. We have, however, attempted to study the profiles 
of four weaker Fe lines in the region 1.6 4 with smaller resolving power and 
we believe that some deductions may be drawn concerning the damping wings 
at least. 

The properties of the lines studied are shown in Table I. These lines have 
been chosen for several reasons: (1) they are free from blends (line 3, only, 
has a weak telluric absorption on the short wave-length wing); (2) they have 
approximately the same excitation potential; (3) the equivalent widths of 
these lines have been measured by Pagel (1955), who found, for the center 
of the sun, results in disagreement with those of Mohler’s table (1955). 


TABLE I 


LINES OBSERVED 





Excitation 


N r Trans. potential 
1 15531 .97 e’D ,—u5 F,° 5.64 
2 15534. 40 e®'D,—u® P»® 5.64 
3 15621 .69 eD,—tD,° 5.54 
4 15662. 10 x8D;°—f PFs 5.60 


B. APPARATUS FUNCTION 
The solar spectrograph of the Dominion Observatory, fed by an 80 ft. 
horizontal solar telescope, consists of a plane grating with collimator and 


1Manuscript received November 21, 1957. : ; 
Contribution from the Dominion Observatory, Department of Mines and Technical Sur- 


veys, Ottawa, Canada, Vol. 2, No. 24. Published by permission of the Deputy Minister. 
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camera mirrors of 20 ft. focal length. The grating has a ruled surface of 
5X3 in. with 14,400 lines per inch. For this investigation the spectrograph 
was used in the first order, double pass (i.e. with the radiation diffracted 
twice by the grating). The solar radiation was chopped at 1080 cycles per 
second before entering the first slit and the output amplified, rectified, and 
recorded on a Brown strip-chart recorder. 

In Table II are reported the dispersions for different grating and chart 


speeds. 


TABLE II 


DISPERSION (IN mm./A) FOR DIFFERENT GRATING AND CHART SPEEDS IN THE 
FIRST ORDER 


Chart speed 








Grating ——- = —— -- ———__ 
speed I II Ill IV V 

I 10.1 13.3 20.2 30.3 40.4 

Il 1.0 5.3 8.1 12.1 16.2 

ITI 2.0 2.6 4.0 6.0 8.1 

3 2.0 4.0 


IV 1.0 bio 


For determining the apparatus function by laboratory observations of the 
CH, band at 1.7 an absorption tube 100cm. long was used. About 150 
tracings of the P(1) line of the 2v; band \ 16,683 A were recorded with different 
pressures and slit widths, the entrance and exit slit widths always being the 
same. According to the results of Boyd, Thompson, and Williams (1952) on 
the splitting of the P branch of the v3 band, we deduce that the inherent width 
of the P(1) line of the 2v; band is negligible in comparison with the presumed 
width of the apparatus function. This experiment was carried out with the 
second grating speed and the second chart speed. 

On each tracing we have determined the widths h, of the profile of the 
line at the depths n(/)>—J,), where n = 0.1, 0.2,...,0.8 and I)—J, is the 
central depth of the line contour. 

All the widths of the 78 tracings obtained under the conditions 


slits widths < 150 yn, CH, pressure < 30 mm., 


could be represented by the same normalized Voigt function with parameters 
B, = 3.6+0.5, B. = 1141.0 (10-® Xd units). 


We have assumed this function as the apparatus function of the solar 
spectrometer of the Dominion Observatory in the first order, double pass, 
second grating speed, in the 1.6 u region. 

In order to check the apparatus function so determined, we have taken, 
with slit widths of 150 yu, 13 tracings of the solar spectrum region between 
\ 1.569 w and X 1.574 w at various altitudes of the sun between 56° and 18°. 
The following combinations were used: 


Icomb. II comb. III comb. 


Grating speed IT] I] I 
Chart speed Il Il II 
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The widths h, of the profiles of the CO: telluric band lines 
15,698.5, -15,709.9, -15,715.8, 15,718.94 
were measured. 
The widths obtained under the conditions 
sun's altitude > 35°, grating speed < II, 
were the same within the experimental errors and may be represented by a 
Voigt function with parameters 
B, = 9.6+0.5, B, = 12.0+1.0 (10-6 \ units). 
These parameters, corrected for the apparatus function, give, according to 


the method suggested by Bell (1951), a negligible Doppler width and a ratio 
Y Yea = 1620 corresponding to 


y/2e = 2.0K 10* sec.—, 
which is identical with the value obtained by Goldberg (1954). 


C. SOLAR OBSERVATIONS 
Tracings of the solar lines of Table | were obtained with slit widths less 
than 150 4 and with the combination 


grating speed I, chart speed I. 


At least three tracings of each line were recorded at distances from the center 
of the solar disk corresponding to cos @ = 1.0, 0.8, 0.6, 0.5, 0.4, 0.3, 0.2. On 
each tracing the following measurements were made (see Fig. 1): 


(1) the height of the continuum Jo; 

(2) the central depth of the line profile Jo—Z,; 
(3) the widths h, of the profile of the line at the depths n(/)—/,), where 
n= GE OZ... RS: 

(4) the area of the line profile. 







y 


Io 
oa 
I. 
ities . ile 


Fic. 1. Diagram showing measurements made on each tracing. 
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The errors in the measurements of the central depths expressed in units of 
the continuum (J)—J,)/Io, widths h,, and equivalent widths W (measured 
with a planimeter) on each tracing were small in comparison with the errors 
due to the estimation of the height of the continuum and to differences 
between tracings. 

In Table III are reported the measurements of central depths and equivalent 
widths. The average deviations are of the order of 1% for the central depths 
and 5% for the equivalent widths. 


TABLE III 


OBSERVATIONAL RESULTS 


In each case the top number represents the equivalent width, W, in A and the bottom number, 
the central depth, (Zo—J.-)/Jo. 


Line cos@ = 1.0 cosé=0.8 cos? =0.6 cos@é=0.5 cos?é=0.4 cosé?é=0.3 cosé = 0.2 





0 
sA 0.135 A 0.130 A 0.133 A 


1 0.105 A 0.102 A 0.120A 0.13 
ons Ue a 0.19 | 0.17 0.19 O17 0.18 

2 0.138A 0.148 A 0.155 A 0.155 A 0.150 A 0.156 A 0.145 A 
0.22 | 0.21. . 0.23 2 0.22 0.20 0.21 

3 0.328 A 0.354 A 0.364 A 0.360 A 0.357 A 0.346 A 0.350 A 
0.35 0.35 0.32 0533 0:32. . 0.29 O.31 

4 0,230 A 0.230 A 0.235 A 0.240 A 0.222 A 0.240 A 0.232 A 


0.27 0.26 0.26 0.26 0.26 0.26 0.24 


We have tried to represent the measurements h, by Voigt functions, choosing 
the form parameter }; by comparing the observed values 


(h, ho.s)ovs 
with those calculated by Elste (1953) 
(hn/No.s) cate: 


The parameters 6; so determined and the corresponding parameters 62, p, 
(oh,,/ho.s) are given in Table IV. The error in the estimation of b; should be 
not more than +0.02. The parameters (1, 62, and A of the Voigt function 
are related to }), be, ho.5, and (I9—T,)/Io by the relations 


Bi = dyho.s, B. = behos, A= bho.s(Lo—I-)/Lo. 


The widths ho. were deduced from all the measured values h, according to 
the formula 


ho.s = LOwoen Do (hin/ho.s) cate: 


In Table IV the three parameters 6), 62, and A are expressed in angstroms 
and also in 10-* \ units for comparison with the apparatus function and with 
the results obtained in different spectral regions. The observations are well 
fitted by the Voigt function and we have not found the systematic deviation 
pointed out by Bell (1951) for medium-strong lines. 

The comparison of our measured equivalent widths, W, and calculated 
equivalent widths, 4, gives an average value 


A/W = 1.10+0.4. 
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TABLE IV 
VOIGT FUNCTIONS FITTING THE OBSERVATIONAL RESULTS 








Bi Be A 
cos 6 db, by Pp (do hn/ho.s)eate hos A 107®X A 107® A 10-®r 
Line 1 
1.0 0.38 0.29 1.44 10.48 0.48 0.18 11.6 0.14 9.0 0.125 8.1 
0.8 0.36 0.31 1.42 10.38 0.44 0.16 10.3 0.14 9.0 0.106 6.8 
0.6 0.36 0.31 1.42 7.83 0.50 0.18 11.6 0.16 10.3 0.135 8.7 
0.5 0.40 0.26 1.46 7.90 0.65 0.26 16.7 0.17 10.9 0.162 10.4 
0.4 0.40 0.26 1.46 7.90 0.54 0.22 14.1 0.14 9.0 0.150 9.6 
0.3 0.30 0.37 1.35 10.12 0.58 0.17 10.9 0.21 13.5 0.133 8.5 
0:2 0.30 0.37 1.35 7.74 0.57 0.17 10.9 0.21 13.5 0.1389 8.9 


Line 2 


1.0 0.28 6:39 1.23 10.04 0.50 0.14 9.0 0.19 12.2 0.146 9.4 
0.8 0.28 0.39 1.33 10.04 0.51 0.14 9.0 0.20 12.9 0.143 9.2 
0.6 0.36 O38 . 1:42 10.38 0.54 0:15. 12:2 O.17 16.9 0.176 11.3 
0.5 0:26 0.41 2.31 9.96 0.64 0.17 10.9 0.26 10.7 0.176. 11.3 
0.4 0.32 0:35 1.38 10.21 0.53 G.17 16:9 O@:19 12.2 0.161 10.4 
0.3 0.30 0.37 1.35 10.12 0.63 0.19 12.2 0.23 14.8 0.170 10.9 
0.2 0.26 6.41: 1.31 9.96 0.58 0.15 9.6 0.24 15.4 0.160 10.3 
Line 3 
1.0 0.40 0.26 1.46 10.58 0.71 0.28 17.9 0.18 11.5 0.364 23.4 
0.8 0.40 0.26 1.46 10.58 0.75 0.30 19.2 0.19 12.2 0.395 25.4 
0.6 0.38 0.29 1.44 10.48 0.85 0.32 20.4 0.25 16.0 0.392 25.2 
0.5 0.42 0.23 1.48 10.68 0.80 0.34 21.8 0.18 11.5 0.381 25.2 
0.4 0.36 6:31 1.4 10.38 0.86 0.31 19:8 0:27 17.3 0.391 25.2 
0.3 0.42 0.23 1.48 7.94 0.80 0.34 21.8 0.18 11.5 0.344 22.2 
0.2 0.40 0.26 1.46 10.58 0.85 O.34 21.8 O22 4.) 0.384 24.7 
Line 4 
1.0 0.30 0.37 1.35 10.12 0.68 0.20 12.8 0.25 16.0 0.248 15.9 
0.8 0.34 0.33 1.40 10.30 0.72 0.24 15.3 0.24 15.4 0.264 17.0 
0.6 0.32 0.35 1.38 10.21 0.72 0.23 14:7 6:25 16.0 0.259 16.6 
0.5 0.30 0.37 1.35 10.12 0.72 0.21. 13.4 @.27 17.3 0.252 16.2 
0.4 0.30 0.37 1.35 10.12 0.73 O22 14.0 O.27 233 0.256 16.4 
0.3 0.36 0:31 1:4 10.38 0.74 0.26 16.6 0.23 14.7 O23 V5 
G:F 1 7 


0.34 0.33 40 10.30 0.73 0.25 16.0 0.24 15.3 0.244 15. 


This discrepancy must be ascribed to the wings at depths less than 0.1(/o—,) 
since for greater depths we have found good agreement between observation 
and calculation. The comparison of our measured equivalent widths, We, 
with those of Pagel, Wp , gives, on the average, 


We ‘Wp = 1.05+0.06. 


This agreement is considered satisfactory. We find between our measure- 
ments for the center of the sun and those of Mohler’s table (1955) the same 
discrepancy which Pagel found. We obtain 


Wo/Wm = 1.381+0.12. 


In Fig. 2 are shown the center-limb variations of log W according to our 
measurements and those of Pagel. The agreement is reasonably good but the 
precision of our measurements is insufficient to enable the most suitable model 
for the solar atmosphere to be chosen from those calculated by Pagel (1955). 
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Pagel 


This Investigation 


a ——* = Te ant as *? line 3 


= : as eat ‘ _tline 4 


i eee ee eT Co ~~ jline 2 
— ; ee ee ee nen 2 E ~}line 1 


el i ! ! peat emer a 
oO 09 0-8 O7 0-6 0-5 0-4 0-3 O-2 
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Fic. 2. Comparison of the equivalent widths obtained by Pagel with those of this investigation. 


D. CONCLUSION 


A treatment of our profiles in terms of anisotropy (Allen 1950; Righini 
1952) or of the formation of the Fraunhofer lines (Houtgast 1952) seems to be 
premature. However, we wish to point out that, correcting the parameters 6, 
for the apparatus function and averaging all the values obtained for the lines 
1 and 2 which, according to the results of Rogerson (1957), are sufficiently 
weak to permit the analysis suggested by Bell (1951), we find for the ratio 
y/Ya the value 210%. This value is greater by a factor 5 than that found 
by Rogerson (1957) for the same excitation potential but for a spectral region 
with larger continuous absorption coefficient. The parameters 82 are so nearly 
the same as those of the apparatus function that, at least for the lines 1 and 
2, no analysis of the Doppler broadening effect is possible; however, their 
trend suggests a dependence of the kinetic temperature on the line strength 
as found by Rogerson (1957). 
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FIELD OF A DIELECTRIC-LOADED, INFINITE CORNER 
REFLECTOR ' 


A. W. ADEY 


ABSTRACT 


A study has been made of the effect on the radiating properties of an infinite 
corner reflector of loading its apex with a section of a dielectric cylinder of 
polystyrene and of a material with a relative permittivity of 4. Calculations 
of the radiation resistance and of the far-field amplitude have shown them to 
be sensitive to the presence of the loading, particularly for small spacings of 
the feeding element from the apex. 


INTRODUCTION 


Several authors have discussed the properties of the corner reflector antenna 
(Kraus 1950; Wait 1954; Moullin 1949). The last author has considered the 
effect on the antenna properties of loading the apex with a metal structure, 
i.e., of truncating it, because of the resulting decrease in cost and wind 
resistance. 

A dielectric loading of the apex could be in the nature of a strengthening 
member and could serve, perhaps, to modify the radiating properties advan- 
tageously.* 

In this report an expression is obtained for the field and the radiation 
resistance in terms of the reflector angle, the position of the feeding element, 
and the parameters of the loading cylinder. Calculations are made for a 90- 
degree reflector, for an element on the bisecting line with spacings from the 
apex up to several wavelengths, and for cylinders both of polystyrene and 
of a material having a relative permittivity of 4. 


THE RADIATED FIELD 

The reflector arrangement is shown in Fig. 1. A current element, of strength 
I amperes, is situated at the position (a, 0) in a cylindrical coordinate system 
(r, @, 3) with the z axis coinciding with the reflector apex and parallel to the 
current element. There is thus only a z component of electric field. The 
reflector angle is 20. The apex is loaded with a section of a dielectric cylinder 
of radius 6. The cylinder material is of permittivity e, permeability un, and 
wave number K (K = 22/wavelength). The corresponding parameters of the 
external medium are €o, uo, Ko. A time variation of e~! is assumed in the 
analysis and rationalized m.k.s. units are used. 

Since the field is confined to the angular space —@ < ¢ < @, one anticipates 





1Manuscript received in original form May 15, 1957, and, as revised, January 14, 1958. 
Contribution from the Radio Physics Laboratory, D.R.T.E., Defence Research Board, 
Ottawa, Canada. Issued as Report 38-3-4, Feb. 1956, of Project PCC No. D48-28-01-03. 

*Some recent experimental results for a dielectric-loaded reflector indicate the broad- 
banding effect of the loading (Natl. Research Council Rept. ERA-296, by J. Y. Wong and 
E. V. Jull, Nov. 1955). 
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Fic. 1, Arrangement of the loading cylinder at the apex of the corner reflector. 


that the term index m for a series expansion of the field will not necessarily 
be integral. Rather, each term for the total field will contain cos v@, where 


(1) cos v9 = 0 or vy = (n+4)7/0 
to satisfy the boundary condition at the reflector surface. 


One can therefore write for the electric field of the current element in an 
unloaded reflector 


— Shel 


(2) Evnc(?, ¢) = F(@,a,r, ) 
(3) a ols x H."(Kea)J.(Kur) coare (r <a), 
(4) = ut > H,(Ker)Jo(Kea) cos vo (r >a), 


n=—a 


where the J, are the Bessel functions of the first kind, the H,“” are the Hankel 
functions of the first kind, and » is given by (1). The factor 7/6 is due to the 
fact that v is not necessarily integral, and reduces to 2 for @ = 2/2, as it 
should (see Appendix I). 

In the usual way, the scattered and internal fields are set up as series of 
Bessel functions of the kinds appropriate to the different regions, with the 
amplitudes to be determined to satisfy the boundary conditions at the surface 
of the cylinder. Thus 


(5) E,(r, ¢) = >. BH, (Ker) cos v¢, 
(6) Em(r, 6) = >> A,J,(Kr) cos v¢. 


na=—ap 


From the Maxwell equation for the ——r of the magnetic field H¢ 


(7) Ht, = —= o i (r, $) 


there results 
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(8) Hie(r, @) = —22 9 3 (Kua) Ji(Kor) cos ¥¢ (r <a), 
(9) ae ~ {nike . HS?"(Kor)J,(Koa) cos vo (r >a), 
o Re () 
(10) H?.(r, ¢) = —— ea B »H,‘”"(Kor) cos vg, 
1) 1K = , r 
(11) Hint (1, ¢) = ae =, A,J;(Kr) cos vd, 


n=—cc 


where the prime indicates differentiation with respect to the argument of the 


Bessel function. 
From the boundary conditions of continuity of E* and H¢ at r = b, one 


obtains 

(12) E ine (0, op) + Beth: >) oe E int (6, >), 

(13) TT ine(b, &) + Hye(b, 6) = Hin (0, >). 

Thus 

(14) A,J,(Kb) —B,H\(Kob) = Hel 5 Hi‘? (Kya) J,(Kob), 


K - Ko — is 


(15) Ay Ji(Kb) —By~ fe (K ob) = HS” (Koa) Ji(Kob), 
yV 


from which 


_ _ enol] J,(0)Hs?"() — Ji@)He ©) a 
ail 46 | a,+ib, H, (w) 


_ Wwpol a ’(w) 


— ~ "260 aybiby’ 
(17) B, = Sholay TS” (w) 
ee ” — 40(a,+ib,) 
where 
(18a) (n= JMO TA ys Ji(u)J,(v), 
ad A. , , 
(18d) b, = J,(u) Yi(v) — J(u) Y,(v); 
By 
(19) u= Kb, v = Kob, w= Koa, 
BS Ky = see) 5 dhe Soe as- Ep. 


The scattered field then becomes, from (5) and (17), 


i ywolr< os ee 
(20) E..(r, @) ae > HT wT (Kor) cos v@ 
a, + 
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and the total field for r > a, from the sum of (4) and (20), 


E(r, ¢) = fee v | ste) 9e Hw) | H}? (Kor) cos vo 

(21) i — feel Y rs =e of) Se Ys | 1K) dobre 

For the far field 

, Dee) 2 y {(Kortre /2—x /4) 
(22) H,’ (Kor) (2. e 
and thus, from (21), 

Be ci os 

(23) E(r, ¢) ~ ~ eho! (=) genre ze T, cos vg, 
where 


b,J,(w) —a,Y,(w) pir 
Gi +ib, : 


For an isolated element the field E,(7, @) is 


; , us ey 
(25) Ejlr, ¢) = a Hy" (Kor) ~ — enol (4) . 


and thus 


( > t ( ’ ) 
9 E ie). 
(26) 3 AP, _—— e T, cos v¢. 


(24) T, = 


i(Kor—m /4) 


RADIATION RESISTANCE 


or an isolated, thin element of radius 6 the field at the surface is given by 
(27) (6) = — 8! HY (Ke) 


and the radiation resistance, R,, by 
Re E(6) _ wno 
I 4’ 


where Re means “real part of’’. This follows since lim Jo(x) = 1. 
x0 


For the loaded reflector the field at the surface of the (thin) line element 
is given by (21) as 


(29) E(a,0) = — Senate , Te" (a), 
Thus the radiation resistance R becomes 

(30) , Re ie: 0) _ onon > W,, 
where 

(31) iy, = (S(w)—a,V(w)]" 


a, +6, 
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and the ratio R/R,, from (30) and (28), 


R rT : 
) eesti es J 
(32) bs Py W,. 


An alternative derivation of this relation will be given in Appendix II. 


CALCULATED RESULTS 
Calculations have been made of the far forward field on the bisector of 
the reflector (@ = 0) and of the radiation resistance, from (26) and (32) 
respectively, for a 90° (@ = 45°) reflector for the cylinder parameters 


Hy = 1; K, = 1.60, 2.00; v = 2, 2.5,3 


for current element separations from the reflector apex up to several wave- 
lengths. 

The results are presented in Figs. 2-5. On comparison of these results 
with those for an unloaded reflector (Moullin 1949), the following points 
are concluded: 

(a) For an unloaded 90° reflector the field amplitude on the bisector varies 
periodically between zero and 4E,, the zero values obtaining for integral free- 
space wavelength separations of the line source from the reflector apex. These 
nulls do not occur for the loaded reflector and the amplitude is not as sensitive 
to the source-apex separation, except for the smallest separations. 

(6) The amplitude minima are shifted from the positions which they occupy 
in the case of the unloaded reflector. 

(c) The periodicity of the amplitude pattern is more complex in the case 
of the loaded reflector, if a periodicity actually exists. Higher values of spacing 
would have to be considered to determine this. 

(d) The radiation resistance fluctuates about that for an isolated element, 
but there is no monotonic falling-off in the amplitude of the fluctuations with 
increasing element—apex spacing, as obtained in the case of the unloaded 
reflector. 


SUGGESTIONS 

The antenna discussed here, while being to a certain extent of academic 
interest, has the same angular pattern as that with a feed element of finite 
length parallel to the apex. In addition, with the recent development of the 
parallel-plate transmission line, an experimental study of this antenna could 
be made. This could include such factors as: 

(a) The angular distribution of the radiated field. 

(6) The effect of covering the conducting plates with, or of replacing them 
by, dielectric sheets. 

(c) The effect of using sections of cylinders of other than circular cross- 
section. 

(d) The use of dielectric loading of the current element as well as of the 


apex. 
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Fic. 2. Far field on bisector of reflector with polystyrene cylinder loading. 


Fic. 3. Far field on bisector of reflector with loading cylinder of relative permittivity e, = 4. 
Fic. 4. Radiation resistance of reflector with polystyrene cylinder loading. 
Fic. 5. Radiation resistance of reflector with loading cylinder of relative permittivity «, = 4. 


APPENDIX I 
Field of an Element in Front of.a Conducting Plane 
By the image principle the field of a tine source of strength J amperes 
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parallel to, and at a distance a@ in front of, an infinite, conducting plane is 
identical with that of the source itself plus that of an image source of —J 
amperes parallel to it and at a distance 2a from it, in the absence of the 
conducting plane. The arrangement of the real and image sources is shown 
in Fig. 6. 






P(r, p) 
ae 


R 
4 ged LINE SOURCE 


I AMPERES 





IMAGE SOURCE 
-I AMPERES 
-@ 


Fic. 6. Arrangement of a line source in front of an infinite, conducting plane. 


Thus the electric field is given by 


ll 


(Al) E(r, 9) — eet (H§(KoR) —Hi(KoS)] 


(A2) = — Hel SHS (Koa) Jn(Kor) 


—(—1)"H\” (Koa) J,(Kor)] cos no (P<), 
(A3) om — esl Sd (A? (Kor) Jn(Koa) 
—(—1)"1I\” (Kor) Jn(Koa)| cos no (r >a), 


which can be reduced to 


(A4)  E(r, 6) = —omol D> HY (Kea)J,(Kor) cos ng (r <a), 
n=1,3,5... 


(A5) = —wpyol 7: H (Kor) Jn(Koa) cos nob (r >a). 


5 


n=1,5,5 


APPENDIX. II 


[lternative Derivation of Equation (82) 
The radiation resistance relation (32) can be derived on a Poynting vector 
basis. From Equation (23) the electric field is given by 


Bs ot ee 
(23) E(r, ¢) ~ nel (2 -) rere > T, cos v¢. 
« 0 n =) 


Thus the @ component of the magnetic field is 
1 0 


(A6) H(r, ¢) = ae E(r, ¢) 


(A7) cad - zy eieormial +i ‘| > T ne o 
ass ~ 9 \8Kor eras ae eer 
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The Poynting vector is defined as 


(A8) S = —}Re(EXH"*), 
where Re means “real part of” and the symbol * indicates the complex 
conjugate. Substituting in (A8) from (23) and (A7) there results 


2 


— 2 
(A9) § =-— Dd W,cos’v¢, 
n=—co 
where, in anticipation of a subsequent integration, only those terms of the 
series product which will give a non-zero integrated value have been retained. 
For the corner reflector the total radiated power is thus 


26 06 
(A10) P= | Srd¢o=r | Sdo 
e 0 


—6 


> W, cos’v¢ do. 


2 6 
_ enor!” 


(All) 


80° J-6 n 

For the line element 
‘ 5 
9 mf oti _ whol [ 2 ) i(Kor—* /4) 

(A12) E,(r) Ls 

are ee 

H,(r) = pom E(r) 

I 2 : i(Kor—m /4) P i | 

(A13) ~ 2) e Kotis ; 
8, = —} Re| Eur; | 
wo l” 
l6mr ~ 


The total radiated power is 


5 _ WHo r 
8 


TR 


(Al4) P, = Zar 
Since for each of the two radiators the radiation resistance is defined as 
(A15) P= [*R 


and we consider J to be the same for the two radiators, we obtain finally, 
from (All) and (A14), 


R rT 
Al6 or aie Ws 
ia R, 0 z 
which agrees with Equation (32). 
REFERENCES 
Kravs, J. D. 1950. Antennas (McGraw-Hill Book Co., Inc., New York). 


Movw.tiin, E. B. 1949. Radio aerials (Oxford University Press, London). 
Walt, J. R. 1954. Can. J. Phys. 32, 365. 











THE TIME-DEPENDENT DIFFUSION THEORY FOR 
CONDENSATION ON SPHERICAL AND PLANE SURFACES! 


J. S. KirKatpy? 


ABSTRACT 


The diffusion field and the temperature distribution are calculated for the sur- 
roundings of a pure solid or liquid phase at constant temperature growing into 
an infinitely extended binary fluid phase of different specific volume. The 
solutions for spherical condensates are applicable to the condensation from air 
of water vapor on ice crystals or water drops. It is demonstrated that the 
generally accepted quasi-stationary diffusion theory for spherical drops leads 
fortuitously to the correct growth law. 


1, INTRODUCTION 
In a consideration of the sublimation of water vapor on ice or condensation 
on water drops, it is usually assumed (Mason 1953; Marshall and Langleben 
1954) that Fick’s first and second laws for the diffusion current and the rate 
of density change, respectively, 


(1) J, = —DVp 

and 

(2) Op/dt = DV*p, 

and the corresponding equations for heat flow and temperature change, 
(3) Jr = —kVT 

and 

(4) aT /at = (k/6Cp)V°T, 


are applicable to the vapor diffusion field and temperature distribution in 
the surroundings. In these, D is the diffusion coefficient and k, 6, and Cp are 
the thermal conductivity, density, and specific heat of atmospheric air, re- 
spectively. All are validly assumed to be constant within the range of applica- 
tion. 

Since the growth rate is observed to be very small, it is assumed further 
that an approximate steady state obtains in the surroundings after a short 
time and that the time derivatives can be set equal to zero (Marshall and 
Langleben 1954). On this basis equation (2) reduces to Laplace’s equation, 
which can be solved in spherical coordinates subject to the conditions 


(5) p=m atr=n, 

(6) P = Pa atr =o 

to give 

(7) P = Pa-(Pa—pi)ro/r. 
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The quantity ro is the radius of the spherical drop or crystal. The rate of 
mass increase of the sphere is written in terms of the ice or water density 
Pp; as 

(8a) dM, a= drr*p,dro ‘dt, 


which by equation (1) also equals 
(8d) 4ar°Dap/dr\ ,=r, = 4arroD(pa—pi). 
Combining the two right-hand sides and integrating, relations are obtained 
for the radius 
Dio. oi 
(9) ro= ee 
V Pr 


and the growth velocity 


dro ‘D (pa — pi) 
10 v= = 4/ eens 


Similarly, the heat transferred at the interface subject to 


(11) P= T; atr = Pro, 
(12) y aay atr =o, 
is 

(13) LdM/dt = 4arnk(T;,—T,), 


where Z is the latent heat of condensation or sublimation. Equations (8) and 
(13) can be combined with a theoretical or empirical interface equilibrium 
relation of the form 

(14) pi = p(T) 


to uniquely determine the growth velocity in terms of the atmospheric con- 
ditions at infinity, pg and 7). 

There are several things wrong with this ‘‘simple’’ approach, both physically 
and mathematically. First of all, for infinite boundaries, a steady state in the 
whole field is never attained since, owing to the large density change at the 
interface, relatively large amounts of solute must be delivered from the 
diffusion field to supply even slow growth. Secondly, the mass motion due 
to the density change at the interface is ignored. Finally, the operations of 
discarding the time dependence to obtain (7) and introducing it again via 
(8) are mathematically untenable. 

In spite of these difficulties, intuition and confirmatory experiments have 
convinced most workers of the numerical correctness of the derived growth 
rate. Accordingly, little effort has been directed toward putting the growth 
law on a rigorous basis. The present paper proposes to do this. 


2. THE LINEAR GROWTH PROBLEM 


As an aid to methodology, consider first the growth of a plane surface of 
pure solid or liquid into a fluid of binary constitution. There is no loss in 
generality if we consider this to be an ice surface advancing into a water 
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vapor — air mixture. Let the boundary conditions be such that the interface 
is positioned at x = 0 at time ¢ = 0 and that 


(15) p(x, 0) = pa, 
(16) p(xo, t) = pa 
(17) T(x, 0) = To, 
(18) T (xo, t) = T;. 


The quantity xo is the instantaneous position of the interface. These con- 
ditions correspond approximately to the case where a slab of ice of temperature 
7; is plunged into the atmosphere at time ¢ = 0 or to the spontaneous nuclea- 
tion of a thin sheet at the time and distance origins. There remains some 
ambiguity as to the interface density p; in the first few instants since the 
initial rapid growth implied by the infinite density gradient suggests in turn 
a large difference of chemical potential for the water in the two phases and 
equivalently, a non-equilibrium density on the vapor side of the interface. As 
is usual in two-phase diffusion problems we will ignore this and assume that 
essentially equilibrium conditions are attained in a time so much smaller than 
the total diffusion time that error is negligible. 

It is desirable to use the number of molecules per unit volume as density 
units. Insofar as the effect of the small temperature change can be ignored, 
this particular choice assures through pressure equalization that the total gas 
density 6 is a constant equal to the sum of the air and vapor densities. This 
is of no consequence as far as the vapor-density distribution and the associated 
D-values are concerned but it does provide that the associated air distribution 
can be obtained from the difference 6—p. The calculation of the general 
solutions of the density equations requires no assumptions about the relative 
magnitudes of the densities involved. However, the subsequent evaluation 
of the growth rates requires appropriate restrictions on the density values. In 
anticipation of this we shall assume that typical atmospheric conditions exist 
such that p; > 6 > pa. 

The first thing to note is that Fick’s laws are not strictly applicable since 
superposed on the mass transfer due to diffusion is a convection towards the 
interface caused by the large density change at that point. This gives rise 
to a velocity term in each of equations (1) to (4). In the next section we 
will solve the n-dimensional diffusion equations containing these explicit 
velocity terms. For the present, however, we will present an alternative 
solution of the one-dimensional problem which is informative. This is obtained 
by first calculating a related diffusion field which is a solution of equation (2). 

Consider a gaseous system with boundary conditions (15) and (16) but 
nucleating and growing on a second pure phase of density equal to the total 
gas density 6 rather than p;. This condition assures that there is no net density 
change at the interface. The solution of equation (2) for such a density field 
has been given by Zener (1949) as 


(19) p= r—(oe—p:) f exp(—A*/4D) dn fexp(—a*/AD) An, 
er a’ 
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where the parameter 
(20) A = «/v1t, 


and a’ is the position of the interface in \-space. If v’ is the instantaneous 
velocity of interface motion, then 


(21) (6—pi)o’ = (6—pi)a’/2Vt = Ddp/dx|-w. 


Evaluating the derivative and assuming at first that a’ < WD we obtain 
, 2 pai 
(22) a Be 


This validates the assumption for the air— water vapor system since the 
coefficient of 1/D is ~10-+. 

It can be seen that the mass distribution given by (19a) is identical with 
that for the desired vapor-to-solid transformation provided we transform the 
interface linearly towards the origin in the ratio of the total density change 
6/p;. The position of the transformed interface in \-space will be at 


(23) a = (6/p;)a’ 
and the velocity of growth will be 

_ 5(pa—pi) = 
24) v 
(  pr(6— pi) 


The associated temperature distribution can be obtained in a similar 
fashion provided the condensing component in the binary phase is dilute 
and provided the specific heat is the same in both phases. This latter con- 
dition can be relaxed if the latent heat obeys 


(25) L> (T1—T,)Cp. 


These conditions are met for the usual atmospheric air-water systems so that 
the temperature distribution after linear transformation is 


(26) r= r.t(T-T) (1 ~ert $+e=2)) / ¢ tse). 


Finally, if we apply the approximate heat balance for the interface, 


(27) k 0dT/dx = —LD dp/dx 

we obtain the psychrometric equation for linear growth 
(28) el ky / Ge 

; Ti-Ta DL L 


This, combined with equations (14) and (24), determines the growth rate. 


3. THE n-DIMENSIONAL DIFFUSION EQUATION WITH DENSITY CHANGE 


The expediency of accounting for the density change by a simple translation 
cannot be invoked for other than linear geometry so we must in general 
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modify the differential equations. This is done by introducing a convection 
term into the current equation, e.g., for linear growth, 


(29) J, = pv’ —Ddp/dax 


subject to the usual continuity equation 
(30) <2 4 2 9, 


These are solved by the following inductive process. Noting that the con- 
vection velocity must be independent of distance, we choose a trial term 





‘ AP us A 2 
where A is an adjustable constant. This leads to the differential equation 
A Op Op °p 

e 2 gaat ~~ = a . 

($2) 2JSt ax + ot Ox” 

Seeking a parametric solution 

(33) p = p(d), 
equation (27) becomes 

: 1=ddp _ dp 

are? 2D dy~ dy? 


with the solution satisfying (15) and (16) 


= Fs iat s | fel x) 
(35) P = pa— (pa pos exp( 1D dn J. EXP\ on Ip dx. 


The quantity a@ is the position of the interface in A-space. Using the transforma- 
tion \’ = A\—A, we can rewrite (35) as 


(36) P = pa— (Pa—pi) : exp(—0" apyan | f exp(—X"/4D) dn 
and comparing this with the linearly transformed solution (19a), 

(196) p = pa-—(pa-—pi) Jo exp(—X" iwyar' / [exp(—n"/4D)an’, 
we see that 

(37) A = —(a’—a) = —a(p,;—64)/6. 


Using the linear treatment as a pattern we can now develop the solution 
for the general case, i.e., the solution of the general diffusion equation derived 
from a current equation of the form 


(38) J, = pv'’—DVp. 


First of all, we know that the convection term in a symmetrical n-dimensional 
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system must decrease as r”~!. We choose, as a trial, a time dependence such 
that the solution has the usual form p = p(A). Accordingly we take 


(39) pv’ = p = aT 


Thence, combining (38) and the continuity equation 


(40) vIn + 2 =0, 

we obtain the differential equation 

ay Ere Rel g + eo 

The substitution p = p(A), \ = r/Vt, leads to the ordinary equation 
2) Lae Sala 

which has for the boundary conditions 

(43) p(r, 0) = pa 

(44) (rns t) = pi 


the solution 
F = ee a ee An 1—n 
(45) P = pa— (pa— pi) i; r exp( 4% fr an—\an / 


co ‘2 
j +" exp( 4s Ef ar—>)an. 


n 


The quantities a, and r, are the interface locations in \ and r space respectively 
and are related by 


(46) t, = aa/t. 


To determine the unknown coefficient A,, reconsider its meaning as defined. 
by equations (38) and (39). By definition, 
fatal 

°° 

is the convection velocity of the whole gaseous mixture relative to the r-origin 
caused by the volume contraction at the interface. Depending on the geo- 
metry this may or may not vary with distance from the origin. However, its 
value right at the interface for a given growth rate a, /2+/t must be independent 
of the geometry. This can be easily seen by comparing the process of con- 
densation of a single interface increment for the various geometries. It follows 
then from equation (37) that 


h(n—2 
A, L An p;—6 


= el? SS ie, oes a 
(47) , Q yl Vt 5 
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evaluated at the interface, or 

—6 
(48) An = —ay" PX, 

6 
As a check, v’’ at the interface can be evaluated from first principles by 
considering in detail the mass transfer at that point. a, is determined in the 
usual fashion by invoking the mass balance at the interface 


(49) (pr—pidon/2Vt = —J,(0:). 


Considering the spherical case in particular we have 


3 2 
5 ee ee as (pr—5) . / 
Se ge nn af, . exo( 206. 4D) 


- 3 2 
©, elaed) 2) 
J exo( 2D6r 4D bis 


and assuming initially that a3 & WD, we obtain as a result of (49) 
/Ree=60, pp 

51 = nT Ey, 

(51) ag V re V 


This validates the assumption since for typical atmospheric conditions the 
coefficient of WD ~ 10-!. The growth velocity is a3/2/t or 


/ D Pa— Pi 
Vx 
ot pr 


(52) v= 
which is identical with the value obtained with equation (10) based on the 
quasi-stationary theory. The reason for this correspondence is not hard to 
see for, since a; K WD, we can approximate p accurately in the region near 
the interface by 


(500) P = pa-—(pa—piasVt/r, 


which is identical with solution (7) except for the explicit time dependence 
of the present result. It is the author’s opinion that this agreement should 
not be used as a justification for the associated quasi-stationary calculation. 
The latter remains essentially incorrect mathematically, and gives the correct 
answer only through a lucky combination of compounded errors and suitable 
geometry. Besides, the correct time-dependent theory is sufficiently concise 
that claims of relative economy for the other are unacceptable. Should an 
even simpler treatment be desired, it is a valid approximation that the con- 
vection term in the current equation be omitted since for a typical atmospheric 
system it has a negligible effect on the magnitude of the density field. It is, 
however, not easy to justify this without reference to the complete solution, 
and even then the phenomenon remains an important item of the physical 
description. 

Our rejection of the quasi-stationary theory as applied to growth problems 
involving infinite boundaries should not be taken as a blanket rejection of 
the method for all boundary conditions. As a matter of fact, the diffusion 
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field about any shape growing very slowly within an enclosure having a 
fixed surface concentration and temperature and the linear dimension / 
attains an approximate steady state within a time after immersion of the 
order of ?/D (i.e., the average time it takes an atom to diffuse a distance /). 
Such systems are fundamentally different from infinite ones since they contain 
a source. These finite systems can be accurately described by Laplace’s 
equation provided sufficient time has elapsed from immersion and provided 
the boundary conditions are well defined. Furthermore, the fields can be 
accurately determined experimentally by analogy with the electrostatic field 
(Jeffries 1918) or by analogy with the flow of electric charge through a con- 
ducting medium in dynamic equilibrium (Marshall and Langleben 1954). 

Mason (1953) has used the electrostatic analogy to predict the growth 
rates of non-spherical snow crystals. No explicit theoretical reference is made 
to the boundary conditions in this work but both the growth equation and the 
associated discussion imply conditions of a semi-infinite medium. It appears 
to the author that the extension of the essentially fortuitous result for spheres 
to non-spherical shapes is not justified at the present state of knowledge, even 
though intuition suggests that such a generalization of the method may lead 
to a good approximation to the truth. 


4. THE TEMPERATURE DISTRIBUTION WITH DENSITY CHANGE 
The associated temperature distribution satisfies an equation analogous to 


(41) provided, as before, that the condensing component is dilute in the 
gaseous phase and provided either or both of 


(53) Cp(phase 1) = Cp(phase 2) 
and 
(54) L > Cp(Ti— To) 


are obeyed. For a typical atmospheric air-water system we can obtain the 
solution 


- ede: wy oe Sei oe J 
(55) T = T,t(Ti-Tha) J nN exp( Axe r a dn 
18 exp AnCrbp 1-1 osx) 
‘ . exp Ok dh “dr 4k dx. 
Applying condition (27) to solutions (45) and (55) with = 3 we obtain 
the approximate psychrometric equation for spheres, 


‘ Pa- pi _ ke 
- ii=-%, 


which of course is identical with that obtained from the quasi-stationary 
theory (Marshall and Langleben 1954). 


5. EVAPORATION FROM SPHERICAL AND PLANE SURFACES 


The current theory for the evaporation rate of water drops in an infinite 
medium is based on the same quasi-stationary theory and equations as used 
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for condensation (Johnson 1950; Best 1952). Unfortunately, the time-depen- 
dent theory, which should be used as a check, is much more difficult than for 
condensation since for spheres we cannot specify at time zero a density dis- 
continuity at the origin, an initial condition which usually leads to simple 
parametric solutions. Fuchs (1934) has obtained an approximate time- 
dependent solution which indicates that the quasi-stationary theory gives 
the correct evaporation rate after sufficient time following immersion (~ one 
second). [t appears, however, that his approximation assumes a stationary 
interface (if not a stationary diffusion field) so the calculation still lacks 
complete rigor. 

The related linear evaporation problem can be solved subject to conditions 
(15) to (18) provided the dense phase has a uniform temperature 7; throughout 
at time zero. In fact, the solutions for both the density and temperature 
distributions are identical with equations (45) and (55) evaluated for n = 1 
provided the signs of the quantities a, and 4, are changed to correspond to 
the reversal of the interface motion and associated mass flow. Since there are 
no quasi-stationary solutions for the linear problem with infinite boundaries, 
it can shed no light on the validity of stationary solutions for spherical 
shapes. 

Luchak and Langstroth (1950) have found time-dependent solutions to 
Fick’s équations for linear and spherical evaporation in enclosures with fixed 
concentrations at the surface. They show that the solutions approximate to 
a high degree of accuracy those given by the quasi-stationary theory provided 
the initial condition is also a quasi-stationary solution. This correspondence 
is not surprising since the initial condition states in effect that at ¢ = 0, 
dp/dt = 0, and this is not going to change unless the interface is moving 
rapidly away from the fixed boundaries. It is clear that any diffusion field 
in an enclosure having fixed surface conditions and deriving from a very 
slowly evaporating surface will attain an approximate steady state and there- 
fore be describable by the quasi-stationary theory. As in the case of con- 
densation, approximate dynamic equilibrium is obtained through the action 
of the sink at the boundary. Fuchs (1934) appears to have been the first to 
work out the stationary solutions for finite boundaries. They have been used 
more recently by Langstroth et a/. (1950) in an experimental check of the 
predicted evaporation rates of droplets in air. 

Considering again in detail the problem of evaporation from a sphere 
immersed in an infinite medium at time zero, we find that there are several 
a priori arguments against the applicability of the quasi-stationary theory. 
First of all, the theory has to be used without reference to the actual boundary 
or initial conditions since the initial condition is already specified by the radius 
of the drop at time zero through equation (7). Secondly, the vapor content 
of the surroundings as indicated by (7) actually decreases (from an infinite 
to a smaller infinite value) instead of increasing as it should. Finally, as for 
condensation, the mathematical procedure used is in question. 

It is still quite conceivable that the correct solution at a sufficient time 
after immersion attains the same space derivative at the interface as given 
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by the quasi-stationary equation (7). The results of Fuchs, abetted by experi- 
ment, physical intuition, and analogy with the corresponding growth problem, 
suggest that this is actually the case. However, it is the opinion of the author 
that the evaporation phenomenon will not have an acceptable theoretical 
description until the calculation is carried out entirely without reference to 
stationary states. 


6. CONCLUSIONS 

The quasi-stationary diffusion theory cannot be validly applied to the 
problem of condensation or evaporation at the surface of a sphere immersed 
in an infinite or semi-infinite medium even though the procedure leads to 
the correct rate of growth for condensation of water vapor on spherical drops. 
It is maintained that all such problems are inherently time-dependent and 
that the latter result is essentially fortuitous. On the other hand, the diffusion 
field about a very slowly condensing or evaporating shape which is enclosed 
by a boundary at fixed composition (thereby acting as a source or sink) must 
rapidly attain an approximate steady state and therefore be amenable to the 
quasi-stationary calculation. 
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ABSORPTION IN A DIELECTRIC SLAB! 


M. P. BACHYNSKI 


ABSTRACT 


The amount of incident energy transmitted, reflected, and absorbed by a 
parallel slab of lossy, high-dielectric-constant material is derived. The results 
are general for any angle of incidence, and calculations are presented for some 
lossy, high-dielectric-constant materials of one-half wavelength electrical thick- 
ness. At angles of incidence greater than 60°, the reflected energy is found to 
be the smaller the lossier the material. The absorbed energy is a maximum at 
certain angles of incidence which depend upon the polarization and dielectric 
constant. 


I, INTRODUCTION 


The amount of incident energy which is reflected and transmitted by a 
slab of dielectric material has in recent years become of increasing interest 
because of its importance in the design of radomes for aircraft and missiles. 
The effects of such a slab of dielectric are here derived. Although in most 
practical situations the transmitted energy is of prime importance, it is 
interesting to note the effect of increasing the loss tangent of the dielectric 
material on the reflected and absorbed power as a function of angle of inci- 
dence. ‘It is found that the reflected energy can in some cases be reduced by 
a lossy material. In addition, the maximum absorption occurs at angles of 
incidence which depend upon the polarization and dielectric constant but are 
independent of the loss tangent of the material. 


Il. THEORY 


Consider a plane wave incident at any angle 6) upon a uniform dielectric 
slab of thickness d (Fig. 1). The angle of refraction within the dielectric is 





Fic. 1. Geometry for calculation of energy transmitted, reflected, and absorbed by a lossy 
dielectric slab. A plane wave in medium of electrical constants €9, uo, Ro is incident at an angle 
6, on a slab of dielectric material of thickness d and electrical constants «1, 41, ki. The angle 
of refraction within the dielectric is 6. 





1Manuscript received November 13, 1957. . 
Contribution from the Research Laboratories, RCA Victor Company, Ltd., Montreal, 
Canada. 
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denoted by 6;. The permittivity, permeability, conductivity, and propagation 
constant of the slab are denoted by «1, ui, o1, and k; respectively. The medium 
surrounding the slab is considered lossless and characterized by the constants 


€o, Ho, and ko. 
The electrical constants of the media are related by 


ke — WE oo, k? = wu +Jwor. 


Henceforth, both media will be assumed non-magnetic so that wo = mw; then 


] €1 j » O1 \ Fi . e 
22) ee —— f= 1 ¢ = . 
(2: €q | itd We} f e {1+ tan 4} ete, 
where w = angular frequency, 
tan 6 = loss tangent, 
jz v1. 


The angles of incidence and refraction are related by Snell’s law: 
ky sin 0; = ko sin Oo. 


The complex transmission factor (7) resulting from a summation of all 
the transmitted rays (Bachvnski 1957) (see Fig. 1) is given by the geometric 
series 


T = exp(jhkid/cos0:) (1 —Tyo?) {1+T oe +Tiote#F+ ...}, 
whose sum is 
(1) T = exp(jkid/cos@;) (1 — Tio?) /(1 — To?e*). 
Similarly the reflected rays are given by 
R = Mati —PMyo*)e#{ 1 +Tire + iote#F+ .. .}. 
The sum of all the reflected rays is 


s (1 oe Pa’ )e*' 
9 oR Soe 
(2) s Pon i. 1—Tyoe”* 
where 


) ky alt 
> — y pe ee a 
é 2d bona ky tand, SIND. ( : 
(Subscript 1 refers to the slab, 0 refers to the region surrounding the slab; 
thus subscript 01 refers to incidence from the outside to the dielectric; 10 
refers to incidence from within the dielectric to the outside.) 

For the electric vector orientated perpendicular to the plane of incidence, 
the I'’s give the ratio of the reflected to incident electric vectors (Stratton 
1941, p. 492). Thus: 
ko CosOo — ky cos; 


(3) I “i = —I| 10 = + cote ky cos6) 3 
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When the electric vector is parallel to the plane of incidence, the I'’s give 
the ratio of the reflected to incident magnetic vectors (Stratton 1941, p. 492): 


ee ki COSA) — Ro cosh, 
(4) fa = -ta= k, cos0y+ko cosd; 


For tan 6 <1, the approximations 


(5) t= Dod (145 , 
Rid _ it er J (1 r )] 
cosh; M etyyp) ti ap) | 
where M = V(e—sin’6o), 


can be introduced. After some algebraic manipulations it can be shown that 
the transmitted energy (7°7*) and the reflected energy (RA*) are given by 


-*[(l—a)? +6") 


inna coucmineaaaian a cient 
(6) f? 1—2¢7 *[a co! a +8 sinbo]+ (a ii en =? 
a a ee i (=f? cosbot+e* : 
7) RR* = 1—2e ° [a cosho +8 sinbo|+ (a°+,° * a 
where - 


= _ ie ) ’ 
~ M 2m?/** 


a kod é 
2 


by — 2kod M. 


a and 6 must be determined separately for each polarization. 
I. Electric vector perpendicular to plane of incidence: 


(p01. —- ;: pat {es +1)’ 


Li > 


(i+. 4° . 
€.} 
( ) B _ 2A (po, +1) — *) 
. (1+-A*)* 
where 


= sin (01-90) 
me sin(6;+09) ’ 


, 
€ 


~ 2M(M-+cos60) 


II. Electric vector parallel to plane of incidence: 


9 


(poy +1 14 2)” — (pop: | 2—A 1) 
(1+4,°)° 


2 (poy +A 1A. 2) (poy 2 — 1) 
(I+.12°)° 


ay, = 


(9) 
B = 
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where 
tan (0)—41) 
pa ™ tan (@o+6;) ’ 


€’ (cos0y— 1/2.M) 


€ cos6y-+ M 
_ €(cosdo+1/2M) 
: € cosOo + VW 
For most practical cases 
A?, A\Ao, A? K1, 
e <i, 
then 
Qa = POi-s ay = Por’, 

= 2p0,A(po.+1), By = 2poy(eoy42— Ai), 

B.*, By = 0. 
Equations (6) and (7) become 

ae € *(1—po’)” Ae ee 

1+po'e °—2e *(po cosbh +8 sinbdo) ’ 


| «= po (1—2e™* cosbote™) 
a aa I+po'e °—2e *(po cosby+8 sinbdo) ' 


(10) TT* = 


where the appropriate po, a, 8 are selected according to polarization. When 
the factor 8 sin bd) is omitted in the denominator of equation (11), the result 
for the reflected energy is in agreement with those of Wolin (1952), who used 
the concepts of impedance (Stratton 1941, p. 511) in approaching the problem. 

The exponential term e~® in equation (10), however, is at variance with 
this author’s result, our exponent ® being {1 —[sin?@o/(e—sin?@)]} times that 
obtained by Wolin. This is believed due to the better degree of approximation 
inherent in this derivation. 

For a lossless medium (tan 6 = 0) equations (10) and (11) reduce to 


varias (l—pe)” 
oh 7 * = > 
\ * l | + po —2po0 cosbhy , 


2po (1—cosbo) 


(13 RR* = 4 
1+ po’ —2py° cosdo 
The energy absorbed by the lossy slab (ZL) is given by the difference between 


the incident energy and that which is transmitted and reflected. Thus: 


(14) L = 1—(TT*+RR*). 


III. DISCUSSION AND CALCULATIONS 


The transmitted energy (77*) as a function of angle of incidence (4) for 
a half-wavelength thick slab with electrical properties « = 8.0 and tan 6 = 0, 
0.01, 0.05 calculated from equation (10) is shown in Fig. 2. Fig. 2(a) illustrates 
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the results for the case of the electric vector orientated perpendicular to the 
plane of incidence and Fig. 2(b) for the electric vector parallel to the plane 
of incidence. In addition, 77* is shown for « = 4.0, tan 6 = 0.05. It is seen 
that when the electric vector is parallel to the plane of incidence, the trans- 
mission remains much greater than for the electric vector perpendicular even 


as the losses increase. 


@ 
2 









— €=80, tons=O 








to — €=80,tons=0 
= 60 F 60 ——e=80, tan 6=O.01 
oO: rane O 2 ‘ 
© e-80,tonse-OOL =‘ Yl €=80, tan s =005 
SG wanew €=80,tans=005 — — €=40,tans = 0-05 
ae | =O; 40 2 
E 40 €=40, tans =O05 d= 2K 
Do 
g 
 .20 20) (b) \ 
E parallel to plane of incidence | 
\ 
0 iene 
0 20 40 60 80 0 20 40 60 80 
ANGLE OF INCIDENCE (6,) ANGLE OF INCIDENCE (6) 


Fic. 2. (a) Variation of transmitted power with angle of incidence for the electric vector 
perpendicular to the plane of incidence. (6) Variation of transmitted power with angle of 
incidence for the electric vector parallel to the plane of incidence. 


The reflected power as a function of angle of incidence is shown in Fig. 3. It 
is found that at angles of incidence near zero, the reflected power is greater 
for a lossy material, but at larger angles of incidence, less power is reflected 
from a lossy slab. The effect is most significant for the #-perpendicular 
polarization but the same general conclusions hold for both polarizations. This 
can be understood by considering the reflected field as the result of two 








Be (a) oe (b) 
may ¢ perpendicular to plane of incidence / E parallel to plane / 
L / of incidence | 
sol —— €=80, tons=0.0) fi ed -— €=80,tons=00! | 
e -eoe- e=80, tans=0,05 ij! ee e*80,tan6=005 | 
Ef —— ergo, wmeraos fff tens009 | 
{5.60 II} 60 If 
= if 
2 I} 
/j! 
/i 
Bag /}; 
o d= \24 /// 
“ y 
a he y 
WwW / 
© 20} Lf 
ok 
70 80 90 
ANGLE OF INCIDENCE (e,) ANGLE OF INCIDENCE (e, 


Fic. 3. (a) Variation of reflected power with angle of incidence for the electric vector per- 
pendicular to the plane of incidence. (6) Variation of reflected power with angle of incidence 
for the electric vector parallel to the plane of incidence. Only the region from 70° to 90° is 
shown. 
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vectors, the directly reflected ray and the resultant of the other reflected 
rays which undergo at least one internal reflection. For a lossy dielectric 
this latter vector is smaller than for a lossless material. Hence, at near normal 
incidence, where the two vectors are nearly out of phase for a half-wavelength 
slab, the lossy material reflects more. As the angle of incidence is increased, 
these two vectors become more and more in phase and consequently an angle 
is reached after which the lossless material will reflect more of the incident 
energy. The reflected power is significantly greater for a lower dielectric con- 
stant slab at high angles of incidence than for one of a higher dielectric 
constant of equal loss tangent and electrical thickness. 

Some interesting results are obtained if the energy absorbed by the slab 
is plotted against angle of incidence. In Fig. 4(a), L = 1—(77*+RR*) 


is shown for the two polarizations for a half-wavelength slab of dielectric 


(a) aca (b) 
— E perpendicular to we ieee 








30) plane of incidence 3 o 
—— E parallel to plane — E perpendicular to 
of incidence f plane of incidence 
—— E parallel to plane 
20) of incidence Brewster ungle 


ABSORBED POWER (L) 





0 20 40 60 80 0 20 40 60 80 


ANGLE OF INCIDENCE (e,) ANGLE OF INCIDENCE (e,) 


Fic. 4. (a2) Power absorbed by a parallel slab of dielectric constant 8.0 as a function of 
angle of incidence. (6) Power absorbed by a slab of dielectric constant 4.0. The positions of 
maximum absorption have both been shifted to smaller angles of incidence than those for a 
slab of dielectric constant 8.0. 


constant 8. The absorbed energy is different for the two polarizations, reaching 
a minimum (as expected) at the Brewster angle for the electric vector parallel 
to the plane of incidence, but a maximum for this angle of incidence in the 
case of the E-perpendicular polarization. A maximum in the absorbed energy 
for the E-parallel polarization is obtained at an incident angle greater than 
the Brewster angle. In Fig. 4(6), the same general behavior is shown for a 
slab of dielectric constant 4.0. The maximum absorption for the electric- 
vector parallel to the plane of incidence is shifted further from the grazing 
incidence angle, while the minimum occurs at the Brewster angle. The maxi- 
mum absorption for the E-perpendicular polarization again occurs at the 
Brewster angle. The maxima and minima occur at the same angles of incidence 
for a particular dielectric constant and polarization independent of the loss 
tangent of the material. 
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SPHERICAL-HARMONICS METHOD FOR NEUTRON TRANSPORT 
THEORY PROBLEMS WITH INCOMPLETE SYMMETRY '! 


B. Davison 


ABSTRACT 


The paper deals with the application of the spherical-harmonics method to 
neutron-transport problems in which the geometrical arrangement possesses 
spherical or cylindrical symmetry, but the distribution of sources and/or the 
conditions at infinity lack such symmetry. Certain classes of cases of this nature 
are considered, and for these, in an arbitrary-order approximation, the solution 
is given, in closed form, in terms of a finite number of Bessel functions. 


1. INTRODUCTION 


In the present paper we consider the application of the spherical-harmonics 
method to neutron-transport problems in which the geometrical arrangement 
is either spherically or cylindrically symmetrical, but the required solution 
cannot possess such symmetry because the distribution of sources and/or the 
conditions at infinity lack it. 

In the cases with complete symmetry, i.e. both the geometrical arrangement 
and the conditions at infinity possess spherical or cylindrical symmetry, the 
solution, as is well known, is given in terms of Bessel functions (Davison, 
1957a, 6). When such symmetry is lacking one is led, generally speaking, to 
solving partial differential equations. The purpose of this paper is to show, 
however, that, when the lack of symmetry is only in the conditions at infinity 
and/or the distribution of sources, it is often possible, just as in the case of 
complete symmetry, to express the solution in terms of Bessel functions, thus 
avoiding the necessity of dealing with partial differential equations. This 
applies in particular to the following cases: 

A. A spherically symmetrical system embedded in an infinite non-capturing 
medium with the neutron flux at large distances given by a polynomial in the 
Cartesian coordinates; 

B. A spherically symmetrical system, either bare or embedded in an infinite 
capturing medium, with distributed isotropic sources in one of the media, the 
source strength per unit volume being given by a polynomial in the Cartesian 
coordinates; 

C. A cylindrically symmetrical system embedded in an infinite non-captur- 
ing medium with the neutron flux at large distances given by a polynomial in 
the two Cartesian coordinates perpendicular to the cylinder axis; 

D. A cylindrically symmetrical system, either bare or embedded in an 
infinite capturing medium, with distributed isotropic neutron sources in one 
of the media, the source strength per unit volume being given by a polynomial 
in the two Cartesian coordinates perpendicular to the cylinder axis. 


‘Manuscript received December 13, 1957. 
Contribution from the Computation Centre, McLennan Laboratory, University of Toronto, 
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The analyses in the cases A and B are quite similar; they are so also in 
cases C and D. Thus it is sufficient to discuss the cases A and C, the modifi- 
cations required in the other two cases being self-evident. 

For the reader’s convenience we first summarize, in Section 2, the main 
results relating to the application of the spherical-harmonics method to un- 
restricted geometrical arrangements and put them into a more readily usable 
form than has been done so far. In Section 3 this is applied to the cases where 
the geometrical arrangement is spherically symmetrical, i.e. to the cases A 
and B, and in Section 4 to the cases where the geometrical arrangement 
possesses cylindrical symmetry, i.e. to the cases C and D. 


2. SPHERICAL-HARMONICS METHOD FOR UNRESTRICTED 
GEOMETRICAL ARRANGEMENTS 
2.1 The Differential Equations Governing the Spherical-Harmonics Moments 
In the one-velocity theory, with isotropic scattering in the laboratory 
system of coordinates, the neutron angular distribution is governed by the 
well-known equation 


(2.1) Q-grad ¥(r, @)+Y(r, Q)/I = (c/4nl) Jfve. Q’)da’ 


(Davison 1957a), where r is the position vector, Q the unit vector in the 
direction of neutron travel, ¥(r,Q) the neutron angular distribution (that is, 
yd VdQ is the probable number of neutrons in the volume element dV travel- 
ling in a direction within dQ), / is the total mean free path, and ¢ is the mean 
number of secondaries per collision. We have purposely omitted the source 
term from (2.1). 

In the spherical-harmonics method Py approximation one approximates 
¥(r,Q) by expanding it in terms of spherical harmonics and neglecting all 
harmonics of order greater than NV. In the cases of complete plane symmetry 
and of complete spherical symmetry (where y¥(r,Q) can involve only one 
component of each particular order), this amounts to expanding ¥(r, Q) in 
Legendre polynomials in the cosine of the appropriate angle. To expand 
Q- grad y(r,Q) in Legendre polynomials, one resorts to the recurrence 
relationships between those polynomials. In more general cases, where y(r, Q) 
may involve several spherical harmonics of the same order, a convenient 
procedure is as follows. 

Let U be the vector of arbitrary magnitude U and direction Q, 


(2.2) U = UQ, 


and let W, be 4rU"/(2n+1) times the contribution of all spherical harmonics 
of order 1 to ¥(r,Q). Then the spherical-harmonics expansion of ¥(r,Q) can 
be expressed as 


(2.3) (1/4) >> (2n+1)Ynlvas, 


where the factor (2n+1)/4m has been introduced to preserve similarity with 
the Legendre polynomial expansion used in the simpler cases. The quantities 








464 CANADIAN JOURNAL OF PHYSICS. VOL. 36, 1958 


W,, are functions of two vector variables r and U. In dealing with such functions 
we shall distinguish the operators grad, div, and V? taken with respect to 
r and U by corresponding subscripts. The fact that W,]y-1 is a linear com- 
bination of spherical harmonics of order m means that W, is a homogeneous 
polynomial of order » in U,, U,, and U,, satisfying the equation 

(2.4) Vyv, = 0. 

Conversely, if we have a homogeneous polynomial of order n in U,, U,, and 
U,, satisfying (2.4), its value for U = 1 will be a linear combination of 
spherical harmonics of order n. To expand Q - grad,(2”+1)¥,]y=1 in spherical 
harmonics, it is sufficient to notice the following. It follows at once from (2.4) 
that Vy?(divy grad,V,,) = 0, while, combining (2.4) with Euler’s identity and 
the fact that V, is a homogeneous polynomial of order in U,, U,, and U,, 
one also gets 


Vyu2{ (2n+1)U - —U* divy} grad, ¥, = 0. 
Thus, if we write 
Q - grad, (2n+1)V,]ya1 = {(2n+1)U - —U? divy} grad, V,Ju-1 
+ divy grad, V,]y-1, 


the first term on the right side consists of spherical harmonics of order (n+ 1) 
and the second of those of order (n—1). Thus, equating separately to zero 
the spherical harmonics of each order in (2.1) gives 


(2.5) divy grad, Wnait (2n+1—cbo.n)Vn/l 
+ {(2n—1)U- — U*divy} grad, V,_, = 0, 


where 69, is the Kronecker symbol. Terminating each of the series involved 
after the spherical harmonics of order N is equivalent to disregarding the 
equations (2.5) for » > N and supplementing those for  < N by 


(2.6) Wit = 0. 
We assume, as is usually done, that N is odd. 


2.2 The Formal Solution of the Equations 

To solve the system of equations (2.5) and (2.6) one can proceed, for 
instance, as follows. By combining (2.6) and (2.5), one can first eliminate 
Vyii, then Vy, then Vy_;, and so on. Thus for each 2 < NV one obtains an 
equation involving only W,, Y,-1,..., Yo; in other words an inhomogeneous 
equation for ¥, with the free term depending on W,_1, Vp-2,..., Wo. Using 
the latter equations one can first determine Wo; then, knowing Wo, determine 
W,; next, knowing Wp and W,, determine W2; and so on. The solution of the 
equation for ¥, will then be given by the sum of the particular solution, 
obtainable by applying the appropriate differential operators to the previously 
determined VY (k < m), and of the complementary function. 

We show now that in this way one obtains for ¥, the following expression: 


n  [}(N—k+1)] 
27) B= > 2 teed "Gace 


k=O j=l 
X Vi.4.(U-grad,; 4U’ Vz") =. ;(r, U), 
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where 
(2.8) Vy..(U-gradr; 3U° V,") 
_ (n — k)! 


P,(u) 
~ (n+k)! due 


i( ad 
y2 2 n—k) 
(UV, ) rae . 
du w= U. grad, /V(u'V ,*) 


the P,,(u) being the Legendre polynomials and the G, .(v) polynomials in 1 /y 
defined by the recurrence relationships 


2n + 1 — cdo, 


(2.9) (n +1 — R)Giuia(vy) + Gnx(v) + (n + R)G,-14(v) = 0 


together with 

(2.10) Gix(v) = 1, Gg-1.4(v) = 0. 

The »,,; for a given Rk and various 7 are the roots of the equation 
(2508) Gyiin(v) = 0. 


Finally Z,,,(r, U), for given & and j, is a homogeneous polynomial of order 
kin U,, U,, and U, satisfying the equations 


(2.12) Vurzx.s = 9, 
(2.13) divy grad, =k j= 0, 
2.14) €¥.* = Vz, 97/1?) gs = 0. 


The factor (n—k)!/(n+h)! in the definition of V,, and the factor (v;,,)~“-” 
in (2.7) have been introduced in order that the definition of G,x.(v) would 
agree with that used in the author's earlier paper (Davison 19576). The roots 
of (2.11) come in pairs—if v = »%,; is a root, so is vy = —v,,; However, only 
one root of each such pair has to be taken into account since (2.14) and (2.7) 
involve only the even powers of vx, ;. 

We denote by [¥,],,; the contribution to (2.7) arising from the term with 
a particular k and a particular 7. To prove the formula (2.7) it is sufficient 
to prove that [W,],,;, for any given k and j, satisfies (2.4) and (2.5). The 
equation (2.6) will then be automatically satisfied in virtue of (2.11). Since 
we are dealing with linear equations with coefficients independent of r, the 
differential operators 0/dx, 0/dy, and 0/dz can be manipulated as if they 
were numbers, as long as we do not substitute the explicit expressions for 
=,,,(r, U). Thus if we denote /0/dx, 10/dy, and 10/dz by qz, q,, and q- re- 
spectively, it is sufficient to verify that* 


(2.4’) Vu { Vax (U-q; 1U” q’) Ex. (8, U)} =0 


and 


*If the Fourier transform of =, ;(r, U) with respect to r exists, the equations (2.4’), (2.5’), 
(2.13’), and (2.14’) are equivalent to the Fourier transforms of the equations (2.4), (2.5) 
(applied for a particular k and 7), (2.13), and (2.14). Our more elementary argument for intro- 
ducing qz, dy, and q: shows, however, that the equations (2.4’) etc. are equivalent to the corre- 
sponding equations (2.4) etc., whether or not the Fourier transform of Zx,; exists. 
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ery _ 1—k r 72 2y nme 

(225°) VE,j ial “C. +1,4(%x, )Q-Srady { Vi+i..(U-q; 1U q ng (F, U)} 
—(n—k) i ie aa 7 
+ (20 + 1 = Cho,n)¥e, 5 Ga. Me, 3) Vn.e(U-G; $U"Q*) Ex. ,(r, U) 
—(n—1-k) 
+ %.3 i "G10 (Ye, 9) [(2n ~ 1)U-q 
— U'q-grady]{ Vi-1..(U-q; $U'9*)| E.,(r, U) = 0 
while the equations (2.13) and (2.14) are expressed now as 
(2.13’) q-grady =,,,(r, U) = 0 
and 
(q?— ve, 97) Sx, 9(F, U) = 0. 
In transforming (2.5) and (2.13) into the forms (2.5’) and (2.13') we have 
also made use of the obvious identity divy grad, = div, grady. In virtue of 
the fact that =,.,(r, U) is a homogeneous polynomial of order k in U’,, l’,, 
and U, by Euler’s identity, it also satisfies the equation 
(2.15) U-grady =,.;(r, U) = k =,.,(r, U). 

To verify that (2.4’) and (2.5’) are indeed satisfied we can proceed as follows. 
Notice that V,,, depends on U,, U,, and U, only via U-q and $U°q? and 
put for brevity 
(2.16) U-q =, 3U°¢? = n, 
so that 

grady ¢ = q, grady 7 = 7U. 
Working out now the left sides of (2.4’) and (2.5’) and using (2.12), (2.137), 
and (2.15) to eliminate the derivatives of Z,,;(r, U) one can readily see that 
(2.4’) and (2.5’) will be indeed satisfied, provided the following relationships 
hold: 
(2.4”’) (Vu° + 2kq’d/dn) V4 (E, 0) = 0 
and 
(2.5"") is a 1a (Ye. 9 18 Vins 1.4/0 St é0 Vast.e/ On] 

2 (2n =r 1 ie C50.) Mx, j rE chi sh Ps k 

—(n—1—k ‘ . 
+ Vx, j "'Gn—1.4 (Ye. s)[(2n or L)EVn-1,% 
"a 2nd Vi-1.%/9 aig 2€0 V»-1,x/9n] = (). 

On the other hand one can easily show, using the definition (2.8) and the 
properties of Legendre polynomials (see Appendix), that V,,,, does indeed 


satisfy the equation (2.4’’). Similarly one can show that 


r 


(2.17) IVn4 1% / OE + EIV 41 e/On = (n+1—k) Vix 
and 
(2.18) (2m —1EVn-14e — 2nOVy-14./ OE — WEIVi-1,4/0n = (N+R)Vi«, 


while (2.14’) implies that, insofar as [W,],,; for a given k and 7 are concerned, 


we should equate qg? to »,,,°. But substituting these into (2.57) we obtain 
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an equation automatically satisfied in virtue of (2.9) and (2.10). This com- 
pletes the proof that (2.7) gives a solution of the equations of Subsection 2.1. It 
can also be shown (Davison 1957a) that (2.7) is the most general solution 
of these equations. 

The question is thus reduced to determining the functions =,,,(r, U). The 
latter should, of course, be determined in such a manner that the resulting 
expressions for V,,(r, U) satisfy the boundary conditions 


(2.19) W,(r, U) continuous 


at the interface of any two adjacent media, as well as the appropriate con- 
ditions at infinity and/or at the free surfaces, etc. 


3. THE CASE WHERE THE GEOMETRICAL ARRANGEMENT 
IS SPHERICALLY SYMMETRICAL 

3.1 The Form of the Solution 

We concentrate now upon the case referred to in the Introduction as 
“problem A”, 1.e. the case of a spherically symmetrical system, embedded in 
an infinite non-capturing (¢ = 1) medium, with the neutron flux at large 
distances given by a polynomial in the Cartesian coordinates. Since the 
equations to be solved are linear and since any polynomial can be presented 
as a sum of homogeneous polynomials, it is obviously sufficient to consider 
the case where the conditions at infinity are 


(3:5) Wo(r) ~ Pm(x, ¥, 3) forr—> &, 


Pm(x, ¥, 2) being a homogeneous polynomial of order m. Further, since in a 
non-capturing medium the flux at large distances should satisfy the equation 
Ve? Vo(r) = 0, the polynomial p,,(x, y, z) appearing in (3.1) should satisfy 
the same equation 

(3.2) Ve7Pm(x, ¥, 3) = 0. 


In the case described, for each of the successive layers of materials, we try 
to determine the =,.,(r, U), introduced in (2.7), in the form 


k 


(3.3) =,.,(r, U) = = Ag. js()Ve [Pm (x, ¥, 2) Ve(r-U; 1°U*)] fork < m, 


=..(r,U) = 0 fork > m, 
where .1,,;-,(7) are some functions of the radial coordinates only, while 
(3.4) V,.(r-U; 4r2U?) = (rU)* P,.(r-U/rU), 


P, being the Legendre polynomial, i.e. Vy is Vg,o of (2.8). 

Recalling (2.7) one can see at once that with Z,,,(r, U) of the form (3.3), 
the flux at any position is given by p,»(x, y, z) times a function of the radial 
coordinate only, and the condition (3.1) can then be readily satisfied. 

Next, it follows at once from the properties of Legendre polynomials that 


(3.5) Vu V.(r-U; Ay? U?) = 0 


and hence (3.3) automatically satisfies the equation (2.12). 
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Finally by substituting (3.3) into (2.14), using Euler’s identity, and equating 
separately the coefficients accompanying V?"(p,V,) for each particular s’ we 
get the following equations for A,,;-5(7): 


(3.6) Fig ee + 2(m+k+1—2s)A's, y:2/7 + Ax,y0-1 = 1, PAr,9:3/l, 
One sek. 


where A,,;,-1(7) should be replaced by zero, since (3.3) does not contain any 
term with s = —1. The general solution of the system of equations (3.6), as 
can be readily verified, is given by 


€ P { @*t Tk — i r l 
(3.7) Axss(r) = y ee bs a | ates ay Fie 4 


T_(min—s— +9) (7Vx,9/L) 
_ —(m+k—s- I 
été .3f °F eg , 


/y\mM+k—s—t+4 
(rv, s/L) 


where the J’s are the modified Bessel functions and at,,;-, and a~,%,j;, con- 
stants of integration. 

It remains to choose these constants so as to satisfy the equation (2.13 
and the boundary conditions (2.19). 


3.2 The Supplementary Equation 

We consider first the equation (2.13), which may be referred to as the 
“supplementary equation’. To facilitate the application of this equation we 
start by rewriting (3.3) in a slightly different form. Put for brevity 


(3.8) r-U =é, 12U2 = 9, 


and notice that, if f(é, 7) is any function depending upon r only via & and 9, 
the following commutation relations hold between 0/0, 0/0n, and V,”: 


a 


(0/0&) Ve = Vr 0/0 + 2U" aEdn’ 


(3.9) 
a 2 Bs 9772 o- 
(0/0n) Ve =Vr 0/dn + 2U 3. 
On 
Using these commutation relations ws the fact that, similarly to (3.5), 
V,.(é, 7) satisfies the equation V,?V;.(é, 7) = 0, one can easily prove by in- 
duction that 


1 
. ry, 20" + B+ A 


(3.10) Voce = -2aet+eyl 


dE" An” B+1 - 


dE*An’ 
Then using (8.2), (3.10), and Euler’s identity one can also easily prove by 
induction that 


(3.11) Vie" din(x, 9,2) Vilé,n)] 


28s! LUO nt ae)! 
(U-grad,)*p,,(x, y, 3)|] saz-z 
(m — s)! d, al(s —a)! (U -geade) 7 ‘ Fy On 


and hence the formula (3.3) for k < m can be rewritten as 
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(3.3) =x, 4(r, U) 
ek ate (a + 8)! (m—a)! 728 
re dX 4 alB! (m-a-— gyi Ae sata (r) U 
x [(U-Brad,)"pa (x,y, 2)] 2a" 
r Pm 131% dE" On’ * 
Substituting this into (2.13), expanding V,(&, 7) and its derivatives in powers 
of & and equating to zero separately the coefficients accompanying 
£((U- grad,)* p,,.(x, y, z)] for each particular a and y, we obtain the following 
set of equations: 
PE ¢ 28 
(2/r) § (@+6+1)! 47+40-2 
9 ea a es ak St) I Be ol ae 
® £.¢-s-pei@—«-— si r A asm) 
(a + B)! , 
ae oo [(k + a)rA's, j:a+8(7) 
+ {2k(m+k+1) — (Qy + 2a — 1)(m+hk — y — a) J Ax, j:0+0(7)] 
r’ (a+pB-—1)! 


2 (m —a— B+ 1) AB — alm — a — B+ Ae sare-s(7)¢ = 0, 


+ 


which have to be satisfied for all a and y in the ranges 
(3.13) 0 <a, boy, a+t2y < k+l. 


The equations (3.12), however, are not mutually independent. On the con- 
trary, if Za,¢,y, say, is the expression appearing under the sign of summation 
in (3.12), one can easily prove by induction that 


y-1 y—1 
Gia 2 Tea we Faas 
p=0 B=q 


(@ + g — 1)1(2/r)** 


* (y—q—1)!g!(m -—a—qt 1)! 


E + q)(m ON See ot 1)A’, jatg(7) 


2 q-1 
ci “om +k-a- a-+- An ses) | re 2s, Ser eh ap Cate Ns 
- a’ 
where Soa’ are certain coefficients, readily determined, but not of direct 
interest for our present purpose. The identity (3.14) shows at once that if 
(3.12) are satisfied for y = | and all a in the range 0 < a < k—1 they will 
be automatically satisfied for all other y within the range (3.13). 

On the other hand substituting in (3.12) with y = 1 from (3.7) and using 
the recurrence relations between the Bessel functions involved, one can readily 
verify that the resulting equations will indeed be satisfied, provided the 
constants of integration a*,,;., appearing in (3.7) satisfy the equations 


(3.15) (R+t—1)Pat,, 5-4-1 + 2t(2t—1)(m4+1—2) yy, Pat... = 0,0 <t < bk, 


i.e., provided 


‘ PVG =— 1's — 4+ 1)! 
(3.16) a> nj: = [ :) (kR+1¢ (m + 1) a ue. CEE 
V;} 


(2t)'(b — 1)'m! 
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where we have assumed, of course, that k does not exceed m. Thus, with 
k < mand a*,,;-, given by (8.16), the equation (2.13) will be satisfied, and 
the final expression for =,.,(r, U), obtained by combining (3.3), (3.7), and 
(3.16), will depend on two arbitrary parameters, a*,,j.0 and a7, j-0. 

The situation would, however, have been rather different had we tried to 
apply the first formula (3.3) for k > m. In that case the equations (3.15) 
would have to be satisfied up to and including t = m+1, and this would 
require that all a*,,,-, for £ << m should vanish. Then by (3.16) ¥e Ax, 3:5 (r) 
for s < m would also vanish, while, for s > m, V,"*[pPm(x, y, 2) V(r + U; 4r?2U?)] 

vanishes. Thus, had we tried to apply the first formula (3.3) for k > m, we 
would have ended with =,,,(r, U) = 0. This is the reason why for k > m 
we took &, ,(r,U) = 0 from the start. 


3.3 The Final Expression for the Moments 

We still have to show that %,;(r, U) as given by (3.3) will allow us to 
satisfy the boundary conditions (2.19). For this purpose, however, we should 
first work out explicitly the expressions for [V,],,, with Z,.;(r,U) of the 
form (3.3). It is obviously sufficient, however, to work out the expression for 


(3.17) Vinx (U-grad,; 3U°V,?)z,,;(r, U) 


(see. (2 . We now show that, with =,., as given by (3.3) and k < m, the 
expression aan is equal to 


min(n,m) 


(3.18) Do Cack.g:s(0)V re" [Pm(x, ¥, 2) Va(t-U; $r°U*)], 


s=0 
where V, is defined by (3.4) with k replaced by m, min(m, m) is the smaller 
of the numbers ” and m and the C,.,.;:.(7) are given by 
K3i19)) Caz y-sl7) 
_ (2k)'(n — k)\(n — 5) J al oo ace 
(n + k)!k! s’=max(0, s+k— my (R- ee ‘\(s EN gi (n a —k- 28. 
7 ia k—s+s' “An eg a 
d(r’/ ~ —k—s+s" - 


gsts’ —k 


<* 


To prove that (3.17) equals (3.18), evaluate each of them ical 
paying particular attention to the term in which U,, U,, and U; enter only 
via r-U, rather than via U? = U,?+ U,?+U,7. Evaluating (3.17) one gets 


(2n)\(2k)'(n — k)! > Sar ae 
9k 


"niin + kik! <4 2k — s)¥l(n —k — 2)! d(r’/2)"* 
x (U-grad,)**' p(x, y, s) + terms divisible by U’, 
which can also be written as 
(2n)! Qk ! oe k ! min(n,m) aie “ 
aT Sage kh) !p! a » (r-U) [(U-grad,) Pm(x, y; z)] 


. s il d"- er Ay. é wr) /d(3r°y"* k—s+s’ \ 
Manele aa "(Rk = s')'(s == si)! (n-—k—- sts’)! 


+ terms divisible by U* 
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On the other hand by evaluating (3.18) one gets 


(2n)! min(n,m) V(r ¢ uy 


—— 2 <~ Ca:t.j:3(r) (U-Sradr)* Pm(x, y, 2) 


2"n! <= = (mn — ss)! 
+ terms divisible by U?. 


Comparing the last two expressions one can readily see that, if the C,:x,;-5(7) 
are as given by (3.19), the entire difference between (3.17) and (3.18) is 
divisible by U*. On the other hand Z,,; has been constructed so as to make 
the expression (3.17) satisfy the equation Vy?(3.17) = 0; (3.18) also satisfies 
Vy"(3.18) = 0 in virtue of (3.5). Their difference then satishes the same 
equation. But it is known that any polynomial, P, in U,, U,, and U,, divisible 
by U?+U,?+U, and satisfying Vy?P = 0 vanishes identically. Hence it 
follows that the expression (3.17) is equal to (3.18) with C,x,j:.(r) given 
by (3.19). 

Once the quantities (3.17) have been evaluated, the evaluation of the 
complete expression for ¥,(r, U) presents no further difficulties (see (2.7)) 
and one can see at once that it is given by . 


min(n,m) 


(3.20) ¥i(t,U) = D> Fas(r)Ve [Pm (x, y, 2) V,(r-U; 


s=0 


r°U’)], 


tol 


where 
min(n.m) [}(N—k+1)] 


(3.21) Fy, s(r) = » xX (L/vx, 5)" nk VE.) Cnt, j:3(1)- 
mt, 

3.4. Some Remarks on the Boundary Conditions 

We are now in a position to discuss the last remaining question, namely 
whether, while taking =,,,(r,U) in the form (3.3), one would be able to 
satisfy the usual boundary conditions (2.19). We restrict ourselves here to 
comparing the number of conditions to be satisfied with the number of arbitrary 
constants available for satisfying them. One should consider separately two 
cases: the cases where the order of the polynomial p,,(x, y, z) appearing in 
(3.1) is less than the order of approximation and is equal to or greater than 
it, ie. m < N, and m > N. We start with the former case. Recalling that 
for k < m the expression for =, ;, for each 7, involves two arbitrary constants, 
while for k > m it involves none, and denoting by [ ], as in (2.7), the integer 
part of the number inside the brackets, we get for the total number of arbitrary 
constants available per medium 


m 


(3.22) » 2i(N — k+1)/2] = (m+ 1)(N +1 — m/2) — [(m + 1)/2]. 
On the other hand it follows from (3.20) that, in order to safeguard the 
continuity of all components of V, for a given m at some interface, it is sufficient 
to ensure the continuity of min(z+1,m-+1) functions of r (namely of the 
functions F,,,(7)) at this interface. If, then, we wish to ensure the continuity 
of all the spherical-harmonics moments, other than those of the highest 
retained order, we have to satisfy 


N-1 


(3.23) 7 min(z + 1,m+1) = (m+ 1)(N — m/2) 


n=0 
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conditions per interface. This is [(m-+2)/2] less than the number of constants 
available per medium. Thus one should be able, using (3.3), to satisfy the 
conditions of continuity of all the moments up to and including those of 
order (.V—1), and furthermore, to satisfy at least half of the conditions 
necessary for safeguarding the continuity of the moments of the order J, i.e. 
of the highest retained order. On the other hand it is known (Davison 1957a) 
that this is as much as one can possibly achieve in the Py approximation. Thus 
(3.3) does in fact give the most general expressions for =,,,(r, U) for the 
problems under consideration. 

If m > N the analysis and conclusions are essentially as above; only m 
in (3.22), in (3.23), and in the sentence following (3.23) should be replaced 
by N. 


4. THE CASE WHERE THE GEOMETRICAL ARRANGEMENT 
IS CYLINDRICALLY SYMMETRICAL 

,.1 Preliminary Remarks 
We turn now to the case referred to in the Introduction as 
i.e. to the case of a cylindrically symmetrical system embedded in an infinite 
non-capturing medium, with the neutron flux at large distances given by a 
polynomial in the two Cartesian coordinates perpendicular to the cylinder 
axis. As-before, it is sufficient to consider the case when the polynomial in 
question is homogeneous. Thus, if the cylinder axis is taken for the z-axis, 


‘problem C”’, 


the conditions at infinity are 

(4.1) Wo(r) ~ p(x, y) for |r| > &, 

where ~Pp(x, y) is a homogeneous polynomial of order m in x and y. Similarly 
to (3.2) it should satisfy, of course, 


( 1.2) Ve7Dm it. y) = 0. 


In the preceding section, where ~,, was a polynomial in three variables, it 
was preferable to treat all the three coordinate axes on the same footing. 
Accordingly, we worked there in general quasi-vectorial notations. Now, 
however, it is more convenient to work in terms of coordinates. We use polar 
cylindrical coordinates both for the position vector r, and for the vector U 
defined by (2.2), putting 


(4.3) x = rcos8, y = rsin 6, 
U, = U; cosa, U, = U; sin ow. 


Notice that in these notations 
(4.4) divy grad, = cos(@—w) {0?/drdU, + (1/rU,)02/d0dw} 
+ sin(@—w) {(1/U,)d2/draw — (1/r)a2/a00U;} + d2/azaU:.. 
Notice also that the definition of p,,(x, y) appearing in (4.1), together with 
(4.2), implies 
(4.5) Pm(x, vy) = (r/ro)™ cos m(6—4o), 


ro and 60) being the constants specifying the polynomial p,,(x, y). 
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4.2 The Solution for the Case Stated 
Under the above conditions the expression for the =,,;(r, U) introduced 
in (2.7) is given, as we verify below, by 


k 8 


(4.6) =,,(r,U) = > — cost (0 — 0) + (k — 2s) (0 — w)] 


x By te—2e1 (U,, U.) {ax, jTm+x—2s (ve, 3/L) + by, Ko mtx os (1Vx, j Ty}. 


Here the J's and K's are modified Bessel functions and the B,..(U,,U.) are 
polynomials of order k in U, and U,, satisfying the equation 


(4.7) Vu" {cos awB,a(U;, U:)} = 0 


and normalized so that the coefficient accompanying U,* is unity. Notice 
that these conditions specify B,.(U,,U-) completely. Finally a,,; and d,,; 
are the two arbitrary constants to be determined from the boundary con- 
ditions (2.19). 

Before proceeding to verify (4.6), we should add a few words regarding the 
polynomials B;.«(U,;, U:). These have been considered in the author’s previous 
paper (Davison 19576) and it has been shown there that they satisfy the 
recurrence relationship* 


(4.8) (0/0U,;ba/U,)Bra(UnU2) = (ke) By-1,0¢1(Ur, U2). 


Also they are closely related to the Legendre polynomials, namely 


Bra = {[d°Px(u) /du"|-0} 'UA(U,? + U2) [@°P. (u) /du" avs sue)- 

We proceed now to show that (4.6) satisfies all the required conditions. The 
equation (2.14) is satisfied in virtue of the properties of Bessel functions, the 
equation (2.12) in virtue of (4.7), and, finally, using (4.4), (4.8), and the 
known recurrence relationships for the Bessel functions, one can easily verify 
that (2.13) is also satisfied. The fact that (4.6) will allow us to satisfy the 
condition at infinity, (4.1), follows at once from (4.5). It remains to verify 
that (4.6) will allow us to satisfy the interface boundary conditions (2.19). 


}.3 The Expressions for the Moments 
For this purpose, we start, as in Subsection 3.4, by evaluating 


V, «(U-grad,; }U°V,?)2,,;(r, U), 


where V,,, is defined by (2.8), and Z,,,; is given by (4.6). Notice that, similarly 
to (4.4), 


U-grad, = U,{cos(@—w)d/dr — sin(@—w) (1/r)d/00} + U,d/dz. 


Hence, using the Bessel functions’ recurrence relationships, one can readily 
see that on applying U-grad, to any expression of the form 


_ *It should be noted that the paper referred to contains a misprint: the right side of the 
first formula (2.10) of that paper should read 
(n+m+1)By mzi( U,, U2) 


and not 
(nm —1)By m#i( U,, U2). 
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n’ 


(4.9) 2; cos[m(@ — 00) + (n’ — 2s) (0 — w)]Imin—2s(rve, /l) Hn, s(U;, U2), 


s=0 


where H,,,.,(U;, U,) are any polynomials of order n’ in U, and U,, one would 
necessarily get an expression of the form 


n’'+1 
(4.10) > cos{m(6 — 6) 


+ (n’ + 1 — 2s) (0 =, @) | Dintn’ +1 25 (1x, 5 /L) My +1, s( (U r U,), 
where 


(4.11) H' 4i,e(U,U2) = (U,/2) {Hy (U,U0) +A -1(UrU2)}, 


i.e. they are some polynomials of the order (v’+1) in U, and U,. At the 
same time operating on (4.9) with V,* will simply multiply it by »,,;?//?. 
Similar remarks apply also to (4.9) with J replaced by A. In this manner 
we see that 


(4.12) (U/v_,))"—* Va (U -grad,; $U°V,2) Ex, ;(r, U) 


for our present =,,,(r, U) will in any case be given by an expression of the 
form 


n 


(4.13) > cos[m(@ — 00) + (nm — 2s)(@ — w)]H*,....(U,, l 


s=0 
, -k - / 
x {Qx, 1mtn 25(1Vx, ;/1) a = 1)” by, jKmtn os(rvx,;/1)}, 


where the factor (—1)”~* in front of A arises only when the Macdonald 
definition of these functions is used, and H*, ,,,(U,,U.) are some polynomials 
of order ” in U, and U,. It is clear, incidentally, from the above derivation 
that these polynomials are independent of m. On the other hand, since 
=,,,(r, U) satisfies (2.12) and (2.13), while V,, satisfies (2.4’’) it follows 
that (4.12), and consequently also (4.13), should satisfy the equation 
Vy2(4.13) = 0. Taking now some m > n, multiplying the last equation by 
cos(m+n—2so), and integrating over all 0, we get 


Vu'lcos(n—250)oH*,x,2(Ur,U:)] = 0, 


which determines /J*, ..,,(U;,U:) to the accuracy of a multiplicative con- 
stant. Comparing this last result with (4.7) shows that it is in fact pro- 
portional to B,,jn-2s5\(U,,U:). Thus (4.13), i.e. (4.12), is equal to 


(4.14) OS D,.x..cos[m(@ — 00) + (n — 2s)(0 — w)]Bn,in-2s|(U;, U2) 


s=0 
/ n—k 'd 
x (x, Im+n—25 (7%, j/l) + (= 1) be. jKmin—25(1%,3/l)}, 
where D,,x,; are certain numerical coefficients, easily obtained for any given 
n, k, and s, though their expression for arbitrary m, k, and s is somewhat 
unwieldy and we refrain from quoting it here. 
Substituting (4.14) for (4.12) into (2.7) one obtains 


(4.15) ¥,(r,U) = : F,, ,(r)cos[m(8—60) + (n—2s) (@—w) |By, n—25)(U, U2), 








cr 
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where 
n  [4(N—k+1)] 


(4.16) F,,2(r) = d 2 Gn. (Ve,9)Dn.x, s{Qx, jTmtn—2s (1%. 3/L) 
7 


+ (—1)" “by, ;Kman—2s(rm,)/D 4, 


so that in order to satisfy the boundary conditions (2.19) it is sufficient to 
ensure the continuity of the functions F,,,(7). The comparison of the number 
of conditions to be satisfied with the number of constants available can be 
carried out exactly as in Subsection 3.4. This completes the proof that the 
=,.,(r, U) for the problem defined in Subsection 4.1 are indeed given by the 
formula (4.6). 
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APPENDIX 
The purpose of this appendix is to prove the relationships (2.4’’), (2.17), 
and (2.18), where V,,, is defined by (2.8) and — and yn by (2.16). Notice that 
the formula (3.10) has been derived from first principles, without any appeal 
to the conclusions of Section 2, so that it is permissible to appeal to (3.10) 
in deriving (2.4’’), etc. Now replacing in (3.10) U by q, r by U, k by n, a 
by k, taking 8 = 0, and recalling that V, is the same thing as V,.0, we have 


2 O°T'Vn.o(€, 0) 


O° Vn.o(€ 0) _ "Vn. o(€, 0 
dt" An 


(A.1) Va ae = — 2kq 
On the other hand it follows at once from the definition (2.8) with the notations 
(2.16) that 
: , _ (n — k)! OV, o(€, 0) 
Viale 9) = oe 
(n + k)! dé 
Combining this with (A.1) we immediately arrive at (2.4’). 

Next, the equation (2.4’’) together with the fact that V,,(&, 7) is a homo- 
geneous polynomial of order m in — and 1/7 defines V,,.(&, 7) to the accuracy 
of a multiplicative constant. On the other hand one can readily verify by a 
direct substitution that the left sides of (2.17) and of (2.18) satisfy (2.4’’). 
Thus the left sides in question should be proportional to V,x(é, 7). Deter- 
mining the coefficient of proportionality from the comparison of coefficients 
accompanying &"-* completes the proof of (2.17) and (2.18). 











THE ELECTRICAL BREAKDOWN OF SMALL GAPS IN VACUUM! 


A. S. DENHOLM 


ABSTRACT 


An investigation of the electrical breakdown in vacuum of small gaps subject to 
contamination by diffusion pump oil showed that the most consistent results were 
obtained when electrodes were conditioned by a hydrogen discharge. Even with 
this method of conditioning a few preliminary sparks were usually required before 
the breakdown voltage reached a plateau level, so that the vacuum discharge itself 
could affect the final values obtained. The circuit parameters which controlled 
the discharge current were found to influence the magnitude and consistency of 
the breakdown voltage appreciably. 

Direct, alternating, and impulse voltage tests showed that the time for which 
voltage was applied to the vacuum gap influenced the breakdown voltage, and 
curves are presented which give the insulation strength of the gap. Two possible 
explanations of the time dependence of the breakdown voltage are given. 


1. INTRODUCTION 


1.1 General 

Vacuum is considered to exist at pressures below about 10~-* mm.Hg where 
the mean free path of residual molecules is large compared with system dimen- 
sions. A charged particle drawn from one electrode to another at such a pressure 
is unlikely to collide with a residual gas molecule. 

Electrons can be extracted from cold metal surfaces by high electric field 
strengths. This phenomenon is known as field or autoelectronic emission, and 
in 1928 Millikan and Lauritzen showed that the electron current varied with 
the field strength according to an exponential law. A theory accounting for 
this emission was proposed by Fowler and Nordheim (1928) and it has been 
improved until it now agrees with empirical results within the limits of 
experimental error (Dyke and Trolan 1953; Lewis 1956). As the voltage 
between electrodes in vacuum is raised, the field emission current increases, 
and there comes an instant when a current increment of several orders of 
magnitude occurs. The current is then limited, in general, almost solely by the 
characteristics of the supply circuit. This sudden large increase of current is 
known as vacuum breakdown or the vacuum spark. 

The experiments described here were not intended as a study of the mecha- 
nism of breakdown, but, rather, to determine if consistent results could be 
obtained under engineering vacuum conditions.* Data which might be used in 
the design of equipment such as vacuum power switches were desired, and 
pressures in the range 10-° to 10-° mm.Hg were maintained. 


1.2 A Review of Breakdown Theories 
The first description of vacuum breakdown seems to be that of Wood (1897). 


‘Manuscript received September 19, 1957. 

Contribution from the Division of Radio and Electrical Engineering, National Research 
Council, Ottawa, Canada. The experiments on which this paper is based were made at the 
Royal College of Science and Technology, Glasgow, Scotland. 

Issued as N.R.C. No. 4669. 

*I.e. Pumping with oil diffusion pumps (without a cold trap) and with a glass—metal system 
using neoprene gaskets. 
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During the subsequent years many have studied the phenomenon using various 
electrode geometries—concentric cylinders, plane electrodes, and point-to- 
plane electrodes—but as the techniques used differed considerably it is difficult 
to correlate the results obtained. 

It is generally accepted that the vacuum arc occurs in vapor produced at the 
electrodes, but the cause of the transition from the field emission to the arcing 
state, that is, the breakdown mechanism, is still uncertain. The requirement 
is a local temperature sufficiently high to produce the necessary vaporization 
for the arc. This may be produced by particle bombardment, by an increase in 
field emission due, perhaps, to change in local geometry, or by both of these 
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=n 
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Fic. 1. Vacuum breakdown mechanisms. 


acting together. The several theories which have been proposed to account for 
breakdown are indicated pictorially on Fig. 1 and for brevity will be itemized 
as follows: 

(a) (Myers and Raatz 1955). Breakdown occurs in the vapor produced at 
the anode by bombardment with field-emission electrons from the cathode. 
Boyle et al. (1955) developed this theory to account for breakdowns at very 
small gaps. They proposed that there is a high yield of electrons at the cathode 
per positive ion produced in the gap, owing to enhancement of the field at the 
cathode by positive-ion space charge rather than by ion bombardment. This 
gives breakdown of the anode vapor in the gap when the multiplication by 
electron avalanche is below that in the more usual forms of Townsend break- 
down. 

(6) (Chambers 1934). Field-emission electrons from the cathode strike 
the anode and produce, by secondary emission, positive ions which are acceler- 
ated towards and bombard the cathode. When the potential and number of 
ions are large enough, a rupture of the cathode occurs which leads to break- 


down. 
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(c) (Ahearn 1936). Prior to breakdown, field currents flow from surface 
projections where the local field and mechanical force are greatest. There is 
local resistive heating, which depends on the size and geometry of these pro- 
jections and their thermal contact with the body of the cathode. As the field 
is increased, a rupture occurs at the projection where conditions of mechanical 
force, resistive heating, and tensile strength are most favorable. Breakdown 
follows. 

(d) (Van Atta et al. 1933). The production of charged particles becomes 
cumulative and breakdown results when AB+CD > 1, where the coefficients 
are as specified below: 

1 electron striking the anode produces A positive ions, 

1 positive ion striking the cathode produces B electrons, 
1 electron striking the anode produces C photons, 

1 photon striking the cathode produces D electrons. 

(e) (McKibben and Beauchamp 1948). A breakdown occurs when AB+EF 
> 1, where A and B are specified as in (d) and 

1 positive ion striking the cathode produces FE negative ions, 
1 negative ion striking the anode produces F positive ions. 

(f) (Cranberg 1952). Loosely adhering material (termed a ‘clump’ by Cran- 
berg) can be detached from an electrode by electrostatic repulsion. Breakdown 
is initiated when this ‘clump’ crosses the gap and strikes the opposite electrode, 
where, it has been shown, local temperatures are produced which are greater 
than any known boiling points. It can be assumed that production of such a 
condition would lead to breakdown. Cranberg suggested that breakdown would 
take place when the energy per unit area delivered to the target electrode 
exceeded a value C’—a constant for any given pair of electrodes. He showed 
that the breakdown criterion became simply VE > C, where V is the gap 
voltage, E is the field strength at the electrode where the ‘clump’ originated, 
and C is a product of C’, some numerical factors, and, possibly, a field- 
intensifying factor due to microscopic distortion of the field in the neighbor- 
hood of the clump during detachment from its parent electrode. For a uniform 
field gap the breakdown criterion becomes V > (Cd)! where d is the gap spacing. 
Cranberg suggested that the clump could come from either electrode, but the 
cathode seemed the more likely source in most cases. 

If mechanism (c) alone determined the breakdown point, the anode would 
have no effect on the sparking voltage. Of those who investigated the effect of 
the anode, Ahearn (1936) and Gleichauf (19516) found it had no effect, 
while the work of Bennett (1931, 1932), Beams (1932), Chambers (1934), 
Anderson (1935), Imachi (1947), Heard and Lauer (1952), and Rozanova and 
Granovskii (1956) indicated that it had. The experiments of Anderson (1935) 
and Trump and Van de Graaff (1947), which showed a total voltage effect, 
indicated that process (c) can not alone account for breakdown at gaps greater 
than about 1 mm. 

The mechanisms (a), (b), (d), (e), and (f) are consistent with the fact that 
the anode affects the breakdown voltage, and with the transfer of anode metal 
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to the cathode at voltages below that necessary for breakdown. This effect 
has been met by Bennett (1931), Anderson (1935), and Heard and Lauer 
(1952). Measurements in recent years (Trump and Van de Graaff 1947; 
Filosofa and Rostagni 1949; Webster et a/. 1952) on the coefficients involved 
in process (d) have indicated that they are not large enough to account for 
breakdown. Measurements under breakdown conditions have not, as far as is 
known, been made on the negative-ion emission coefficients involved in process 
(e). 

The curves presented by Cranberg (1952) give strong support to his ‘clump’ 
theory of breakdown under uniform field conditions, particularly at voltages 
above about 50 kv. Under non-uniform field conditions the evidence in support 
of Cranberg’s theory is not conclusive. The experiments of Chambers (1934), 
which would seem to be the classic breakdown study on coaxial systems, do 
not positively confirm this theory, although they show that the condition 
of the anode definitely influenced the sparking voltage. Dyke and Trolan 
(1953), Dyke et al. (1953), and Dolan et al. (1953) showed that the emission 
of current from a point could alone cause breakdown by resistive heating, but 
the nature of their experiments (microsecond voltage pulses) was such as to 
eliminate the other possible effects which would depend on the anode. 

From this survey and the results in the next sections of the paper the follow- 
ing mechanism is suggested as the cause of breakdown at comparatively small 
gaps (< 1mm.) and high surface gradients, if the positive-ion transit time is 
short compared with the time of duration of the gap voltage. With a gap voltage 
close to the breakdown value, there are, at microscopic projections on the 
cathode, high electrical and thus mechanical stresses. Owing to field emission, 
a high current density exists which produces resistive heating of the projections, 
the temperature being increased also by the bombardment of the point by 
positive ions produced at the anode by the electron stream. These conditions 
of tensile stress and high temperature lead to a fracture of the weakest pro- 
jection, which initiates breakdown. The extent to which the sparking potential 
depends on the current density is determined by the gap voltage (ion energy) 
and the number of positive ions produced by the electron stream, which is 
influenced by the material and condition of the anode surface. 

Consider now the experiments of the two workers who found no anode effect. 
Ahearn’s method of conditioning the anode was by the process of spark break- 
down itself, but it is likely that only a small fraction of the anode surface 
would be outgassed in this way, whereas the cathode could possibly be con- 
ditioned by the removal of projections by melting (see later). This possibility 
is strengthened by the fact that when a previously shielded part of the anode 
was placed opposite the cathode a slightly higher breakdown voltage was 
obtained, rather than the lower breakdown voltage expected by Ahearn if the 
anode influenced the sparking voltage. Gleichauf found that fusing a film 
of glass on the anode did not affect the breakdown voltage. This is not com- 
pletely incompatible with the proposed mechanism if the emission of positive 
ions was in that case influenced largely by the adsorbed gas. 
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2. APPARATUS: ELECTRODES AND CIRCUITRY 


The electrodes were disks of 1} in. diameter with curved opposing surfaces 
(4 in. radius) which gave approximately uniform field conditions over the 
central region where the sparks occurred. The analysis of the electrode 
materials is as follows: 


Steel—Fe, 99.518%; Mn, 0.380%; C, 0.052%; S, 0.05%; 
Copper— 99.5% pure, impurities unknown; 


Aluminum—Al, 99.64%; Si, 0.23%; Fe, 0.1%; Cu, 0.08%; Mg, 0%. 


The surfaces of the steel and copper electrodes were normally finished with 
a polishing buff, but no significant difference in their performance was found if 
they were finished with fine emery cloth. A bright surface finish was obtained 
on aluminum electrodes by polishing slowly in a lathe with smoothed 3/0 
emery. On the other hand, finishing with a liquid metal polish made these 
electrodes smoother, but duller, and presumably an oxide skin was formed. 
Both types of surface were examined and are designated ‘bright’ or ‘dull’ 
aluminum. Electrodes were always washed in ether before mounting in the 
vacuum system. 

In all the alternating, direct, and impulse voltage experiments a relay discon- 
nected the supply as soon as the gap broke down, so limiting the energy dis- 
charged through the gap. The discharge circuit was essentially a reservoir 
capacitor, series resistor, and the gap. The capacitor was 500 yyf. except where 
stated during the impulse study, and the series resistance was varied from 
0 up to 500 kilohms. The impulse generator used had only one stage and is 


shown in outline in Fig. 2. 


VACUUM 
HIGH GAP 
IMPEDANCE 
SUPPLY 





Fic, 2. Basic impulse circuit. C; = 0.01 uf., Co = 500 puf., C3 > Co. 


3. CONDITIONING OF ELECTRODES IN A VACUUM SYSTEM 
WITHOUT A COLD TRAP 

Metals normally have adsorbed on their surfaces gases which affect their 
emission characteristics, but these gases may be removed by suitable treat- 
ments which often also remove microscopic surface projections. Possible 
conditioning processes were tried using steel electrodes and alternating voltages 
(50 c.p.s.) up to 70 kv. peak. 

The feasibility of conditioning by the vacuum spark itself was examined, 
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but although the breakdown voltages eventually reached a plateau, they 
depended on the value of the series resistance and were erratic. The maximum 
values attained were still considerably below those with well-conditioned 
electrodes. 

Induction heating of electrodes to a temperature of 800° C. (bright red heat) 
in vacuum was tried, but the subsequent breakdown voltages took even longer 
to reach a plateau level than when there had been no conditioning. Apparently, 
oil vapor from the pump cracks on the hot electrodes leaving minute carbon 
deposits. 

The following hydrogen discharge technique, in which current densities 
similar to those of Anderson (1935) were used, gave satisfactory conditioning. 
With the gap set at 2 mm. and a few hundred volts across it, hydrogen was 
leaked into the evacuated system until the discharge, which initially filled 
the whole chamber, covered only the opposing surfaces of the electrodes (at 
about 10 mm.Hg). The voltage was then adjusted until the current was 1 
amp. (for steel; for other metals see below), which was maintained for 3 
minutes. 

This process was the best of the methods tried, and the figures obtained from 
a typical breakdown run are shown on Fig. 3. Normally, a plateau was reached 
after a few sparks, and standard deviations of the order of 4% could be 
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obtained with an occasional low breakdown voltage which was disregarded 
(less than 80% of the maximum value was considered a low freak). 
Discharge currents of 0.5 amp. for copper and 0.25 amp. for aluminum gave 
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satisfactory conditioning. A current of 1 amp. produced sputtering with copper 
electrodes and sufficient heat to melt the surfaces of aluminum electrodes. 

Figure 4 shows, for steel, copper, and aluminum, curves of maximum break- 
down voltage vs. gap spacing, which were obtained when alternating voltage, 
increasing at a rate of approximately 6 kv./second, was applied across the 
gap. In this paper, unless otherwise stated, electrodes were conditioned by a 
discharge in hydrogen. 


4. INFLUENCE OF THE DISCHARGE CIRCUIT ON BREAKDOWN VOLTAGE 

The constants of the discharge circuit by influencing the spark current can 
affect both the level and the consistency of the breakdown voltage. A large 
resistance in series with the gap produced erratic, or, sometimes, just low 
breakdown voltages and the most consistent results were obtained when the 
resistance was nominally zero. Two tests, one with, and one without, a series 
resistor are compared in Table I. With series resistance below about 50 kilohms 
the effects are not so great, but apparently the tendency is still there, as can 
be seen from the results plotted on Fig. 5. 


TABLE I 
‘THE EFFECT OF SERIES RESISTANCE 
(Steel electrodes, 0.4 mm. gap, applied voltage 50 c.p.s.) 


Maximum Mean Standard 
Series B.D:N., BDV. deviation, Low freaks, 
resistance kv. (peak) kv. (peak) ® mean B.D.V.  % total B.D.'s 
0 53 1S .6 1.11 oot 
100 kilohms 19.5 38.7 21.6 


st8) 


(KV) 


VOLTAGE 





BREAKDOWN 


SERIES RESISTANCE - OHMS 


Fic. 5. The influence of series resistance on the magnitude and consistency of the direct break- 
down voltage (steel electrodes, 0.5 mm. gap). X, mean breakdown voltage; @, standard 


deviation, ©; mean. 


Although reasonably consistent results were obtained with a_ reservoir 
capacitance of 500 wuf., it is not known that this is the optimum value, but 
impulse experiments showed that a considerably higher capacitance gave 
erratic results (later). The addition of a comparatively small capacitance 
(17.5 ppt.) immediately across the gap noticeably increased the scatter of the 


results. 
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i. 6. Cathode craters. 
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It is thought that the spark can be considered as having two main stages: 

(a) the discharge of the gap and associated capacitance; 

(b) the discharge of the reservoir capacitor. 
The discharge of the gap capacitance and the breakdown mechanism itself 
produce a certain roughening on the cathode, which can be smoothed by 
melting if the ‘follow’ current from the reservoir capacitor is large enough. 
The greater the cathode roughening, the more likely is a low breakdown 
voltage to follow. Fig. 6 shows typical cathode ‘craters’, one produced by a 
spark with a high resistor in series with the gap, and one with no resistor in 
series with the gap. The cathode roughening produced in tests with high 
resistance was confirmed by a pre-breakdown current several times greater than 
in those experiments with low series resistance (when pre-breakdown current 


was normally less than 50 pa). 


5. DIRECT-VOLTAGE OSCILLOGRAPH STUDIES OF THE SPARK CURRENT 
The voltage across a comparatively small resistor connected between the 
low potential electrode and ground gave a record of the spark current (Fig. 7). 
At breakdown, the capacitance to ground of the high-voltage electrode and 
connections discharged, and this was followed by the discharge of the reservoir 





Fic. 7. Spark current—5.0 amp. peak, 0.95 amp. chop, copper electrodes, gap 1 mm., R, = 11 
kilohms, breakdown voltage 58.2 kv. direct, timing oscillation 0.2 Mc./second. 


capacitor, which was controlled primarily by the series resistance. The dis- 
charge of the gap capacitance took place in 0.1 ysecond or less. The spark 
current extinguished about 0.5 amp., which is called the chopping value, and 
if the series resistance was increased to a value such that the maximurn current 
possible from the reservoir capacitor through the gap was of the order of the 
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chopping value, or less, the breakdown consisted only of the discharge of the 
gap capacitance. Occasionally, several such ‘suppressed’ breakdowns were 
obtained on the one oscillogram. 

The minimum, mean, and maximum chopping currents obtained with 
different metals and gap spacings are shown on Table II. There it is seen that 
for a given metal the variation of chopping current with spacing is not great, 
nor, in comparing different metals, is it consistent. It seems likely that, if 
sufficient measurements were made, the mean chopping current for each metal 
at different gap spacings would be the same (at least over the range considered 
here), and that the actual value at any one breakdown depends on the local 


conditions at the arc roots. 











TABLE II 
VACUUM GAP CHOPPING CURRENTS (RX; = 11 KILOHMS) 
Electrode Min. chop, Mean chop, Max. chop, No. of 
metal Gap amp. amp. amp. measures 
Steel 0.2 0.19 0.51 0.86 10 
Steel 0.3 0.14 0.52 LS 5 
Steel 0.5 0.13 0.31 0.78 15 
Copper 0.5 0.33 0.74 2.0 12 
Copper 1.0 0.48 0.67 2.3 14 
Aluminum 0.5 0.14 0.29 0.58 15 
Aluminum 1.0 0.17 0.35 0.72 15 





6. DELAYED BREAKDOWN AND THE DIELECTRIC STRENGTH 
OF THE VACUUM GAP 

Figure 8 shows the maximum breakdown voltages obtained with direct, 
alternating, and impulse voltage* applied across a steel electrode gap. The 
impulse voltage figures resulted from incremental testing, and were not de- 
rived from oscillograph measurements on the wavefront with the gap over- 
voltaged; frequently breakdown occurred well down the impulse tail. 

It is evident that there is a time effect, as previous investigators have noted. 
There were two obvious methods of studying this effect. One was by finding 
the breakdown voltage for different rates of rise of voltage, but results obtained 
in this way would be of little practical importance and fundamentally depen- 
dent on the figures supplied by the other method, which would determine the 
time between the application of a fixed voltage and breakdown. Consequently, 
the following procedure was adopted. 

A direct voltage, which was about 25% below that at which the first break- 
down was expected to occur, was applied for 1 hour, then increased by a small 
amount and left for a further hour. The gap voltage was raised in this way until 
a spark occurred, when a relay removed, then reapplied, the voltage across 
the gap. The time when this and subsequent sparks took place was noted, 
and the voltage increased in steps until the rate of sparking was about two 
per minute. As the rate of sparking increased, the time during which measure- 


*Direct voltage rising at 6 kv./second, alternating voltage peak rising at 6 kv./second, 
and impulse waveform 12/50 (the unusual front time was due to energy considerations, and 
will be discussed later). 
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ments at a given voltage were made was decreased, to prevent, as far as 
possible, excessive surface damage to the electrodes. 
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One hour was decided upon as the initial time for which the voltage should 
be applied, as any appreciably longer period would give tests lasting for days 
and was impracticable in the circumstances. The form of the curves of break- 
down rate vs. voltage subsequently obtained suggested that if the breakdown 
did not occur in | hour it was unlikely to occur in 1 day. 

Before each test the electrodes were resurfaced, and, initially, hydrogen 
discharge conditioning was used, but it was later found that, as far as the 
initial breakdown voltage was concerned, no improvement was obtained by 
conditioning in this way, although conditioning was reflected in a lower rate 
of sparking at the higher voltages. 

Figure 9 shows three rate-of-sparking curves obtained for steel electrodes. 
A spark could influence the succeeding rate of sparking by producing surface 
roughness. This was shown clearly in some of the tests when rate of sparking 
increased with time at a given voltage, and in the test corresponding to curve 
B, Fig. 9, it was so markedly increased after one particular breakdown that 
two values are given for the rate at 41 kv. Curve C (points 4 and 5) shows that 
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a rougher surface had produced a greater number of sparks per hour, but it 
also indicates that the voltage still had an appreciable influence on the rate. 
That is, the rate of sparking did not increase with voltage merely because of the 
increasing number of preceding sparks. 

Table III gives, for various tests on steel, copper, and bright aluminum, 
the highest voltage which gave no breakdowns in | hour, and the voltage 
above this which gave the first breakdown. The mean of these figures is con- 
sidered to be the insulation strength for that particular gap. On a few occasions 
a spark occurred at a voltage below that which subsequently gave no break- 
downs in 1 hour (e.g. Fig. 9, curves A and C), but this spark was neglected 
and considered merely as a form of conditioning. Fig. 10 shows the mean 
insulation strengths obtained for the three metals tested. 


TABLE Ill 
INSULATION STRENGTHS, STEEL, COPPER, AND BRIGHT ALUMINUM 














Max. with- Insulation 
Gap, stand voltage, First B.D., strength, 
Metal mm. kv. kv. kv. 
Steel 0.25 20 22 21 
20 2 21 
24 26 25 
20 22 21 
0.5 34.5 38.6 36.5 
4.0 38.6 36.5 
34.5 38.6 36.5 
30:0 40.7 38.1 
35.5 40.7 38.1 
0.8 45.9 51.0 48.4 
49.0 52.0 50.5 
51.0 52.0 51.5 
45.9 49.0 47.4 
Copper 0.5 25.2 30.3 27.7 
2a:2 30.3 ee 
20.2 S033 Sek 
1.0 58.1 61.3 58.7 
51.0 56.1 oa:0 
51.0 56.1 53.5 
Bright 
aluminum 0.5 32.4 36.6 34.5 
32.4 35.6 34.0 
32.4 36.6 34.5 
1.0 46.0 49.0 47.5 
40.7 42.9 41.8 
40.7 44.0 42.3 
0 


46. 48.5 47.2 


Aluminum electrodes with the ‘dull’ finish gave higher, but more inconsistent, 
insulation strengths than those with the ‘bright’ finish (Fig. 10), and with 
aluminum, in general, often the voltage at which the first breakdown took 
place was the highest that the gap would withstand—subsequent breakdowns 
occurring at lower voltages. 
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INSULATION STRENGTH — KV 





GAP (mm.) 


FIG. 10. Insulation strength. X, steel; @, copper; A, bright aluminum; R, range for dull 
aluminum. 


7. IMPULSE BREAKDOWN 

Initially, overvoltage impulse tests were carried out with breakdown occur- 
ring on the front of the wave. It was found that 0.1/50 waveforms gave break- 
down voltages about 25% above those given by 1/50 waveforms, but testing 
with fast wavefronts (i.e. small R,, Fig. 2) gave very erratic and low results 
owing to the excessive energy discharged through the vacuum gap by the 
0.01 yf. capacitor. 

Previous tests with alternating and direct voltage showed that fair results 
could be obtained if only a 500 ywyf. capacitor was discharged through the gap 
at breakdown, and consequently tests were made with Ry increased to a 
value (43 kilohms) such that the discharge of the 0.01 yf. capacitor through 
the gap was suppressed (Section 5). This gave, as before, a reservoir capaci- 
tance of 500 wyf., and at breakdown resulted in a repeated impulse waveform 
as shown in Fig. 11. Sometimes the gap broke down on this repeated impulse. 

Using the 12/50 wave given by R, = 43 kilohms, impulses of increasing 
amplitude were applied to the gap until it broke down, whereupon a relay 
disconnected the supply to the impulse generator. With steel and copper, 
approximately 25% of the sparks happened at less than 80% of the maximum 
voltage in a test, but with aluminum the first breakdown voltage was usually 
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Fic. 11. 12/50 Impulse waveform with breakdown. Timing oscillation 50 kc./second. 


the highest of that test, and on Fig. 12, which shows maximum impulse 
breakdown voltages, the highest values obtained for aluminum in all the tests 
made have been indicated by the dotted curve. The impulse tests confirmed the 
findings of Section 3 with regard to hydrogen discharge conditioning. 


8. CONCLUSIONS 
8.1 General 

In an oil-pumped vacuum system without a cold trap electrodes may be 
conditioned by hydrogen discharge to give reproducible breakdown voltages.* 
A brief comparison of the results presented here with those of previous 
investigators who also determined breakdown voltages under uniform, or 
near uniform, field conditions is given in Table IV. The differences are marked 
and, but for steel, in sense contrary to what might be expected, bearing in 
mind the fact that the author’s results were obtained with electrodes con- 
taminated by pump vapor. The use of a cold trap increased by 8% the 
maximum breakdown voltage for steel electrodes at a gap of 0.5 mm. 

If a vacuum gap breaks down, the voltage at which subsequent discharges 
take place is influenced by the circuit parameters owing to the effect of the 
current on the cathode geometry at the arc roots. In particular, a resistance in 
series with the gap, which is large enough to prevent burning of the arc, 
appreciably lowers the next breakdown voltage. 

The current at which the vacuum arc extinguishes is appreciable, as has 
been noted previously (Gleichauf 1951; Leader 1953). Measurements of this 


*It is thought that in the experiments described here, where the electrodes were subject to 
contamination by pump vapor, the major effect of the conditioning process was to remove 
microscopic projections from the electrode surfaces. 
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Fic. 12. Maximum impulse breakdown voltage (12/50 wave). X, steel electrodes; @, 
copper electrodes; - - - aluminum electrodes. 


TABLE IV 
.\ COMPARISON WITH THE RESULTS OF PREVIOUS INVESTIGATORS 


Breakdown 


Gap, voltage, 
Investigator Metal mm. Voltage form kv. 
Denholm Steel 1.0 50 c.p.s., rate of rise 6 95* 
kv./second 
\nderson (1935) Steel 1.0 Direct, rate of rise unknown =122 
Trump and 
Van de Graaff (1947) — Steel 1.0 Direct, rate of rise unknown 110 
Denholm Steel 0.3 Impulse (12/50) 59 
Leader (1953) Steel 0.3 Impulse (1.5/70) 14 
Denholm Copper 1.0 50 c.p.s., rate or rise 6 68.5 
kv./second 
Anderson Copper 1.0 Direct, rate of rise unknown 37 
Denholm Aluminum 1.0 50 c.p.s., rate of rise 6 78 
kv./second 
\nderson \luminum 1.0 Direct, rate of rise unknown 41 


“By extrapolation. 


extinguishing value give some tenths of an ampere for the three metals exam- 
ined and suggest that it is strongly influenced by the geometry at the arc 
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roots. It was proposed by Leader that if the relaxation oscillation experienced 
in his investigation with impulse voltage and high series resistance occurred 
with direct voltage, as was found here (Section 5), it could account for the 
conflicting results obtained by some previous investigators. 


8.2 Delayed Breakdown 

It has been shown that the voltage at which a vacuum gap breaks down 
depends on the temporal form of the voltage and that breakdown can occur 
after a voltage has been applied for some time. The time effect has been 
mentioned in the literature but it does not seem to have been examined closely, 
although some breakdown voltages have been given as the voltage at which 
sparks occurred at a fixed rate (McKibben and Beauchamp 1948; Rozanova 
and Granovskii 1956). 

The long times which may elapse before breakdown exclude the possibility 
that chain mechanisms (e.g. (d) and (e), Section 1.2) account for the break- 
down, which would appear to be initiated by a single event. This leaves the 
Cranberg theory and that postulated under (c), Section 1.2, both of which 
require a reason for the time delay in tearing a projection from an electrode 
(presumably the cathode). Two possible explanations are offered. 

It has been suggested (Section 1.2) how a small projection on the surface of 
a cathode can experience high tensile stress and temperature—these conditions 
leading to fracture and the initiation of breakdown. Under such conditions 
there occurs a process known as ‘creep’ which is due to viscous flow at grain 
boundaries or crystal slip and which describes a normally slow extension of a 
stressed specimen ending ultimately in fracture. This creep, the speed of which 
depends on the severity of the conditions, may be the process producing dela yed 
vacuum breakdown. Gossling (1932) has also suggested this reason for delaved 
breakdown and he cited the work of Goucher (1924) with single crystal tung- 
sten wires. Goucher showed that at constant temperature the time to fracture 
(ty) at tensile stress (S,) was given by an equation of the form 


K 
(1) igs 


ep St—m)/n? 


where A, m, and n are constants. 

Since the tensile stress at a surface varies with the square of the electrical 
stress (Attwood 1949) an attempt was made to correlate the time of application 
of voltage with the square of the electrical stress at breakdown. From the results 
given in Fig. 13 a logarithmic plot corresponding to expression (1) was made of 
the square of the maximum breakdown stress against time for which stress was 
above 90% of its maximum. We chose this fraction arbitrarily, bearing in mind 
the fact that the maximum stress would have the greatest effect. It was assumed 
that when alternating voltage was applied across the gap the conditions in each 
half-cycle were not appreciably influenced by phenomena in previous half- 
cycles, or in other words, the voltage form was considered to be a pulse of a 
single half-cycle. This is suggested by the fact that in mechanical creep testing, 
a specimen returns nearly to its original condition when the stress is removed. 
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On the above basis the times of application of voltage ranged from 14 yseconds 
in the case of the 12/50 impulse wave to 0.9 second in the case of direct 
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Fic. 13. Maximum breakdown stresses and insulation strength, steel: (i) impulse voltage, 
12/50 waveform; (ii) alternating voltage, 50 c.p.s.; (iii) direct voltage, rate of rise 6 
kv./second; (iv) insulation strength. 


voltage rising linearly. Obviously, the best method of determining the relation- 
ship between time to breakdown and electrical stress is to apply rectangular 
or unit function waveforms, but even with the rough approach used here the 
breakdown stress fitted to 5% the following expression, which is of a similar 
form to (1): 


; 14.2 
(2) tap = =———— seconds, 
e(Se”—9700)/8700a 





where S, is the electrical stress in kv./mm. and a = 0.51—0.88d (d is gap 
spacing in mm.). 

The following estimate was made asa check on whether the electrical stresses 
encountered in a vacuum gap could account for rupture in the correct order 
of time. A macroscopic electrical stress of 80 kv./mm. (cf. plot of insulation 
strength of Fig. 13) and an assumed intensification due to a projection of 12.5, 
which is produced at the tip of an ellipsoidal column of height/base ratio 
1.8 (Race 1940), give a mechanical stress of 650 Ib./in.2 Unfortunately, no 
stress rupture information on low carbon steel, such as was used for the elec- 
trodes here, has been found, but from data given by Orr et al. (1953) on a 
killed carbon steel, which presumably is stronger, a time to rupture of 1 hour 
is obtained at a temperature of 960° C. for a stress of 650 Ib. /in. 
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Another explanation for the delay in breakdown is based on the extrusion 
of projections by surface migration under the action of the electric field and/or 
ion bombardment. From equations given in Dyke and Dolan (1956) it can 
be shown simply that the rate of increase in length of a projection (dz/df) 
due to surface migration is given by an equation 


(3) dz_Gie%™ fee %,) 
dt ick a RS’ 
where S is the electric field, 
T is the temperature, 
R is the radius of the tip of the projection, 
y is the surface tension, 
Q is the activation energy, 
k is Boltzmann’s constant, 


and C,, Co are constants. 


Since the actual dimensions and geometry of the surface projections which 
produced breakdown in these experiments are not known, it is difficult to 
determine which of the above two explanations is the more likely to account 
for delayed breakdown. 
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CRYSTAL ORIENTATION IN ICE SHEETS! 
F. G. J. Perey? AND E. R. POUNDER 


ABSTRACT 


Crystal size and orientation in ice sheets frozen under one-dimensional cooling 
from melts of pure water or melts containing traces of organic additives were 
studied using a polariscope with a universal stage. The surface layer is found 
to consist of crystals with nearly vertical optic axes. Horizontal sections cut at 
various depths show a gradual change of the mean orientation of the optic axes 
towards the horizontal. The change to horizontal inclinations is almost completed 
at depths of 4 or 5cm. under the freezing conditions used, and occurs more rapidly 
with the melts containing additives. With additives the ice sheet shows another 
feature also, a transition layer about 0.5 cm. below the surface. In this layer the 
crystals are much smaller than in the surface layer above them and the number of 
crystals extending through this layer is small. 

An explanation of the observations is offered, in terms of preferred growth of an 
ice crystal in planes perpendicular to the optic axis. This provides a mechanism 
permitting more rapid growth of crystals with inclined axes. 


INTRODUCTION 


In a previous paper, one of the authors (Pounder 1958) reported the low 
mechanical strength of ice sheets frozen on the surface of a melt containing 
traces of alcohol and one of several organic compounds characterized by long 
chains of carbon atoms. The striking difference in the appearance of these 
impure ices from that of pure ice suggested that the lowering of the strength 
might be a result of differences in crystal size or orientation or both. The study 
of the crystal habit of ice here reported was thus undertaken in the hope of 
obtaining further insight into the role of impurities in ice formation. It may 
be said at the outset that few differences were found in the crystal structure of 
pure and impure ice. The mechanism of ice growth appears much the same for 
pure and impure melts, with the impurities merely accentuating the peculiar 
plate structure of ice crystals. The explanation of the weakness of certain 
impure ices almost certainly lies in the behavior of intercrystalline boundary 
layers, as previously reported, but the technique used in these experiments shed 
no light on the nature of these layers. On the other hand the results obtained 
are of some interest in connection with the general mechanism of ice formation. 

Ice consists of optically uniaxial, birefringent crystals. For light of wave- 
length 5500 A the indices of refraction of the ordinary and extraordinary rays 
are respectively 1.3106 and 1.3120. The birefringence is thus small, but it is 
sufficient to display the crystal structure readily when a thin section is exam- 
ined between crossed polaroids. Crystal sizes can be obtained with the simplest 
polarizing equipment and the orientation of the principal optic axis or c-axis 
of an ice crystal can be found using some form of universal stage with several 
degrees of freedom of rotation. A number of studies of this nature have been 
made on the ice of lakes, rivers, and the oceans by, among others, Leydolt 
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(1851), McConnel and Kidd (1888), Tamman and Dreyer (1934), Seligman 
(1937), Wilson, Zumberge, and Marshall (1954), Butkovich (1956), and Weeks 
(1957). The reports of these authors are frequently in disagreement but may 
be summarized as saying that ice crystals in a natural ice sheet have strongly 
preferred orientations of their c-axes. The preference may be for vertical 
orientation or for horizontal orientation, in different ice sheets, but never for 
intermediate angles. Explanations for preferred orientation have been offered 
by Bertin (1878) in terms of thermal gradient, and by Miigge (1895) in terms 
of the action of gravity on the initially floating single crystals. Neither explana- 
tion appears to cover both orientations. As the present paper will show that 
the direction of the preferred orientation is a function of depth, it should be 
pointed out that the reports listed above are rarely explicit as to the exact 
location in the ice sheet from which the sections studied were cut. 

Some difficulty exists in defining a single crystal of ice. If it is taken as one 
that shows a uniform color under crossed polaroids (and this seems the most 
useful working definition), a single crystal may have an internal structure 
consisting of parallel plates or layers. This layer structure has been reported 
by Faraday (1860), McConnel (1890), and more recently by Nakaya (1956) 
and Rigsby (1956). The thickness of the layers has been variously given as 
from “indefinitely thin’? to up to 2 or 3 mm. Nakaya observed a mean layer 
thickness of 0.061 mm. with values ranging from 0.042 mm. to 0.070 mm., 
and Rigsby, studying the deformation of a single crystal under mechanical 
bending, observed that the spacing of the gliding planes agreed with the 
thickness found by Nakaya. 


EXPERIMENTAL METHOD 

Samples were frozen in plastic tanks in a deep-freeze unit at an ambient 
temperature of —30° C. The sides and bottom of the tanks were heavily insu- 
lated with foamed plastic insulation of low thermal conductivity to ensure as 
completely as possible one-dimensional cooling. Heat-flow lines are mapped in 
the ice by lines of air bubbles and showed that, with one exception to be dis- 
cussed later, the flow of heat was essentially vertical except for a very thin 
layer at the edges and bottom of the blocks of ice. Samples were about 30 cm. 
by 30 cm. in area and from 6 cm. to 20 cm. thick. After the block of ice was 
frozen, it was removed from the tank and further work was done in a cold 
room whose temperature was regulated at —15° C. The center of the block 
was cut out with a band saw in a prism of square cross section of side 10 cm., 
and this prism was then cut into horizontal slices with saw cuts at intervals of 
about 0.5 cm. Each slice was then frozen on to a glass plate about 10 cm. square 
and 1 mm. thick. A moderately warm glass plate melts a very thin layer of ice, 
which then freezes, making an excellent bond. Fortunately, the water film 
freezes with the same orientations as the crystals in contact with it, if the 
freezing is done from the ice and not from the glass plate. The ice on the glass 
plate was next thinned down by melting with a warm metal plate and then 
shaving with a small plane until its thickness was about 1 mm. The surface 
was then ground with fine emery paper to a thickness of 0.5 mm. and finally 
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polished with a leather glove. The thin section was then photographed between 
crossed polaroids. 

Determination ol crystal orientation was limited to measuring the azimuthal 
and polar angles of the c-axis, although it might be useful to measure the secon- 
dary axes of the crystals as well. For this determination the section was 
mounted on a universal stage with four axes of rotation. This stage was designed 
by Dr. G. P. Rigsby and its four degrees of freedom permit the measurement of 
the orientation of any crystal regardless of its polar angle. Two methods are 
used. In the first the crystal is rotated until the direction of propagation of 
the light in the crystal is parallel to the c-axis. The crystal is then isotropic for 
the propagation of light, and, the polaroids being set for extinction, the crystal 
appears dark. In this position the crystal is said to be ‘‘standing up’. Total 
reflection limits this technique to polar angles of less than 45°, where the polar 
angle is the angle between the c-axis and the normal to the plane of the section. 
For larger angles the second method is used. The crystal is rotated until its 
c-axis is horizontal and parallel to the fixed horizontal axis of the stage. This 
procedure is referred to as “laying down” the axis. In the standing up case the 
polar angle as read on the stage must be corrected for refraction of the light 
as it goes through the thin section. This correction is readily obtained from 
Snell’s law. Our analysis shows that no refraction correction is needed for the 
angle observed in laying down the axis. It should be pointed out that there is 
some disagreement on this point. Rigsby gives a correction ranging up to a 
maximum of about 7° for a polar angle of 40°, in using the laying down method. 
This correction, if necessary, would alter our quoted results very slightly 
but would not affect the general argument. 


RESULTS 
Pure Water 

Three samples were frozen from distilled water to give ice sheets of thickness 
6 to 8 cm. As analysis of the samples showed changes in the distribution of 
orientation with depth, two further freezing tests were made. A deeper tank 
of distilled water was frozen to a depth of 20 cm., freezing taking place at an 
ambient temperature of —30° C. In these four tests cooling was essentially 
one-dimensional throughout the volume of the tank. Finally, a much larger 
tank (60 cm. by 60 cm. by 40 cm. deep) was filled with tap water and frozen 
(for 2 weeks) in the cold room at —15° C. The thickness of heat insulation 
was relatively much less in this test and, although the center part of the ice 
block which was analyzed showed a one-dimensional freezing pattern, there had 
been considerable freezing from the sides and bottom so that there must have 
been continuous convection currents in the melt. 

The first four samples were analyzed in detail for crystal sizes and orienta- 
tions. In most cases the observations for a particular section were first plotted 
on a Schmidt equal-area projection. Fig. 1 shows the Schmidt chart for the 
surface section of one sample. In this chart each point represents a single 
crystal, plotted with azimuth and polar angles of the c-axis as coordinates, 
the center representing a polar angle of 0° (vertical crystals) and the circum- 
ference one of 90° (horizontal crystals). From Fig. 1, and from many others 
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Fic. 1. Schmidt equal-area projection showing orientations of crystals in the surface layer 
of a typical ice sheet. 
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Fic. 2. Composite histogram showing the fraction of the horizontal area covered by crystals 
of a given orientation. The two histograms show mean results for four samples at the surface 
and at a depth of 5 cm. 
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plotted, it is clear that there is no preferred azimuthal direction. The Schmidt 
chart gives no indication of crystal sizes, so the next step was to count the 
crystals with polar angles in a given range, multiply each by its area, and sum 
to find the percentage of the area of the section covered by crystals with 
c-axes in each 10° interval. Histograms prepared in this way showed that in all 
four samples the crystals in the surface layer showed a strong preference for 
vertical optic axes, but that with increasing depth there was a strong shift in 
the statistical distribution in favor of crystals with optic axes nearly horizontal. 
Fig. 2 shows a composite histogram giving the average results for the four 
samples at the surface layer and at a depth of 5 cm. Although the distributions 
are shown for two levels only, a number of intermediate sections were studied 
in all cases and showed that the transition in the mean direction of the optic 
axis is a continuous one as a function of depth. The fourth sample studied was 
thicker, and histograms for the surface and for depths of 5 and 13 cm. are 
plotted in Fig. 3. Note that in the first 5 cm. of growth most of the vertical 
crystals are replaced by ones with axes in the range 40°—90°, but that in the 
next 8 cm. of growth there is considerable attrition of the crystals in the 
40°-80° range with strong growth in the 80°-90° interval. Apparently in the 
initial growth, after the surface layer is formed, any crystal with an inclined 
c-axis tends to grow at the expense of vertical crystals, but as the latter dis- 
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appear the more nearly horizontal crystals are favored over the ones with 
intermediate inclinations. 

The final large sample presents a quite different story. The surface crystals 
were large, about 8 cm. in diameter in a horizontal section, and vertically 
oriented. Horizontal sections to a depth of 25 cm. showed the same crystals 
extending down all the way from the surface, and no evidence of crystals 
with inclined optic axes was found. 

No consistent pattern of crystal sizes was found in the surface layers of the 
five ice samples studied. The crystals obtained in the fifth test were an order of 
magnitude larger than those in the other four cases, as would be expected be- 
cause of the slower rate of freezing. However, the first four tests showed a wide 
range of crystal sizes, despite apparently identical temperature conditions. 
In a given freezing test the surface crystals are relatively uniform in size, but 
the average dimensions were about four times greater in the sample showing 
the largest crystals than in the one with the smallest crystals. Fig. 4 shows the 
pattern of a surface layer with crystals of intermediate size. Evidence from 
other experiments in this laboratory suggests the existence of a quasi-station- 
ary, insulating layer of air over a freezing surface of ice. The thickness of the 
laver, which represents a density inversion, is not yet known but the phenome- 
non is being investigated further. It is probable that the ambient temperature 
as measured several centimeters above a freezing surface does not give a very 
good indication of the rate of heat transfer from the ice. 


Additives 

A number of freezing tests were made in which small amounts of ethyl 
alcohol, or ethyl alcohol and a solution of a solid additive, were added to the 
distilled water before freezing. In some tests the additives were applied in such a 
manner as to concentrate them in the surface and in others they were diluted 
uniformly throughout the volume of the liquid. The various solid additives 
included sodium carboxymethyl cellulose, soluble starch, polyvinyl methyl 
ether, sodium stearate, methyl! cellulose, sodium chloride, and lithium chloride. 
Results of this series of tests can be summed up briefly as follows. Alcohol 
changes the freezing pattern slightly from that of pure water, particularly in 
the surface layers. None of the solid additives caused any further observable 
change in the crystal structure. Surface concentration of the alcohol and volume 
dilution produced ice of similar appearance and crystal structure. 

In a typical test with alcohol as an additive, 10 ml. of alcohol were sprayed 
on the surface of the tank after the 6 liters of water in it had been cooled to 
below 4° C. After freezing, sections were prepared as described above. In this 
case, as in all tests with additives, the structure of the ice sheet showed a 
discontinuity at a depth of about 5 mm. Fig. 5 shows a vertical section through 
the ice sheet with this discontinuity clearly evident. Above this layer the ice 
shows the characteristic strong preference for crystals with vertical orientation. 
Surface crystal sizes are indicated in Fig. 6. The transition layer, shown in 
‘ig. 7, is so completely chaotic that it was impossible to make a study, with 
the universal stage, of the distribution of the polar angles of the crystals. 
Below the transition layer the number of crystals decreased with depth and, 
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as in pure water, the average direction of the polar angle of the optic axis 
swung steadily towards the horizontal. In fact the change to the horizontal is 
more rapid with additives present than in pure water. Analysis of Fig. 2 
shows that the weighted mean polar angle of the crystals is 63°, at a depth of 
5 cm. below the surface in pure ice. The ice samples containing additives were 
not analyzed in as great detail but for a single, and apparently quite repre- 
sentative, sample the corresponding figure for the mean inclination of the 
polar axis is 77°, at a depth of 3.5 cm. 

With ice samples containing additives, each crystal usually shows a sub- 
structure of thin parallel plates varying in thickness from about 0.5 to 1.0 
mm. The optic axis is always normal to the plane of these plates. Figs. 8 and 9 
show vertical and horizontal sections through a sheet of ice grown from a melt 
containing traces of alcohol and sodium carboxymethy] cellulose. The plate 
structure is quite evident. In these additive tests, if the ice sheet is removed 
from the melt, the under surface of the ice is not smooth but shows the edges 
of the plates clearly. They protrude different distances into the melt and some 
of the dendritic growths may extend as much as a centimeter out from the 
average lower surface of the ice. The optic axes of several of these dendritic 
leaves were measured and were invariably normal to the plane of the leat. 


DISCUSSION 

In all, about fifteen freezing tests were made using either pure water or 
water containing additives. The most striking feature of the tests is the transi- 
tion from mainly vertically oriented crystals in the surface to ones with almost 
horizontal inclination at a small depth. This was observed in all but one case, 
the test in which a large volume of water was frozen slowly with relatively 
poor heat insulation. 

The following explanation of the change in crystal orientation is suggested. 
The initial formation of ice is in the form of lacy ferns or dendrites which spread 
rapidly over the surface. The optic axis of a dendrite is always perpendicular to 
the plane of the dendrite. As freezing continues the interstices of the dendrite 
fill in until it becomes a flat plate. Further growth may be of two kinds, by 
attachment to the corners or edges of the dendrite, increasing its area, or by 
attachment to its plane surface, increasing its thickness. The former may be 
called edge growth and the latter columnar growth. Thermodynamically, edge 
growth will be favored because less energy transfer is involved in the transition 
from a “‘free’’ molecule in the liquid to a molecule bound toa few lattice mole- 
cules at a corner or edge than in the change from the liquid to a more firmly 
bound position on a completed plane of the lattice. Another way of looking at 
this is that the corner or edge presents more holes in the lattice for a new 
molecule to fit into than does the plane side of the crystal. 

When a body of water is below 4° C. and heat is removed through the surface, 
the top layer will become supercooled since no convection takes place. Even- 
tually some nuclei of condensation become activated. These nuclei will have 
random orientation but those which are vertical or near vertical will grow 
rapidly since the favored growth plane coincides with the supercooled layer. 
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This accounts for the large preponderance of vertically oriented crystals in 
the surface layer. After the surface layer is completely formed, vertical crys- 
tals can only increase by columnar growth whereas crystals with inclined 
optic axes may still grow by edge accretion. This favored growth will be 
retarded by the temperature gradient in the water but crystals of this type 


| can be assumed to project at least a few molecules farther into the melt than 
those extending by columnar growth. Fig. 10 shows how this process can 


result in inclined crystals growing in horizontal extent at the expense of ver- 
tical crystals. The dotted lines indicate the plate structure of the crystals. The 
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Fic. 10. Preferred growth of crystals with inclined optic axes, resulting in gradual extinction 
of a vertically oriented crystal. The lower sketch represents a later stage in freezing. 


lower part of the sketch shows a later stage in the development.§The middle 
crystal, B, is now cut off from the melt and can grow no more. Crystal A is 
growing more rapidly than crystal C because the dendritic planes of A are 
more nearly vertical. This mechanism explains the observation that the 
crystals with vertical c-axes gradually disappear, being in large part replaced 
at moderate depths by crystals with axes of intermediate inclination and these 
in turn by others with their c-axes essentially horizontal. 

This suggested explanation applies equally to the tests with ice frozen from 
melts containing additives but here an additional factor comes into play, 
namely the rejection of impurities by ice crystals. As impure water freezes, 
the growing ice sheet pushes a layer of impurities ahead of it. This is a type ot 
zone refining although freezing was so rapid in these tests that layers or 
pockets of impurities were trapped between crystals in the ice sheet. In the 
previous paper by one of the authors already referred to, the existence o! 
mechanisms for surface concentration of organic impurities was pointed out. 
When an ice sheet forms, this surface concentration combines with the impurity 
rejection by the ice to form a layer of sufficiently high impurity concentration 
as to inhibit freezing completely. As further heat is removed the water below 
this very impure layer cools to the freezing point and new crystal nuclei will be 
activated there, starting a second layer of ice. Eventually the impure water 
separating the two ice lavers will freeze as the temperature in the ice continues 
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to drop. It is suggested that this is the origin of the chaotic, transition level 
observed in all the samples containing additives. 

Below the transition level there are now two factors favoring the growth of 
crystals with inclined or horizontal optic axes. It is probable that a consider- 
able fraction of the impurities are trapped at the transition level. Hence the 
impure concentration layer ahead of the advancing ice surface will only be a 
thin film and the preferred edge growth of crystals will be still more highly 
favored because the edges of the dendritic planes will project through this 
film into less contaminated water than that in which columnar growth is 
taking place. The rapid transition of the optic axis towards the horizontal 
under these circumstances has been pointed out already. 

It remains to explain the exceptional case, the one in which vertically oriented 
crystals flourished. As noted already, the insulation on the sides of the large 
vessel used in this test was relatively scanty. Study of the heat-flow pattern, 
as mapped by the air bubbles, showed that considerable ice grew in from the 
sides of the container. This would result in extensive convection currents and 
thorough mixing of the water. This in turn would result in a very large tem- 
perature gradient in the layer of water immediately adjacent to the bottom 
of the ice sheet. It is suggested that this gradient was sufficiently great as to 
nullify any inherent advantage of edge growth so that inclined crystals grew 
no more rapidly than vertical crystals. The surface preponderance of the latter 
would thus be maintained throughout the ice sheet. 

If correct, this suggests that the transition of the optic axis towards the 
horizontal is a result of one-dimensional cooling and may not occur when 
convection currents or other causes result in mixing of the water. It is possible 
that this is the explanation of the fact that lake ice is found sometimes with 





vertical axes and sometimes with horizontal ones. 
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THE THERMAL-NEUTRON FISSION CROSS SECTION OF Pu*?! 


C. B. BIiGHAM 


ABSTRACT 


A measurement of the thermal-neutron fission cross section of Pu®*® is described. 
rhe value obtained, —0.8+0.7 barn, is compared with other measurements. 


INTRODUCTION 


Recent measurements of the thermal-neutron fission cross section of Pu*® 
are in serious disagreement. Hulet, Bowman, Michel, and Hoff (1956) have 
reported a value of 4.4+0.5 barns. The measurement of Leonard, Seppi, 
Friesen, and Kinderman (1956) at the 1 ev. resonance indicates a thermal 
value of about 0.05 barn. A low value is supported by earlier unpublished 
values of <0.2 barn (Fields, Friedman, and Bentley) and <0.5 barn (Fried, 
Bentley, and Pyle) quoted by Jaffey, Hibdon, and Sjoblom (1955). This 
paper reports a measurement which also supports a low value. 


METHOD 
Basis 


If isotopically pure samples of Pu?** and Pu*® were available, the ratio of 
the a-disintegration rates and the published ratio of the a-decay half-lives 
would give the ratio of the masses of the samples. The ratio of the fission rates 
in these samples, when irradiated in the same neutron flux, would then give 
the ratio of the fission cross sections. The accuracy would be limited to the 
accuracy with which the ratio of the a-decay half-lives is known. The situation 
in the present measurement is much less favorable in that the best available 
sample of Pu’ contained about 10° Pu*’*, which contributes a fission cross 
section of roughly 100 barns per atom of Pu*®. This means that an uncertainty 
of 0.6°% in the ratio of the Pu? and Pu**® @ half-lives introduces an error of 
about 0.6 barn (i.e. 0.6% of 100 barns) in the value obtained. 


Ex perimental 

In order to reduce the errors in the @ and fission counting well below that 
arising from the @ halt-lives, the following technique was employed. The 
fission rates were measured in a neutron beam that was uniform to better than 
+1% over an area somewhat larger than the sample area. Plutonium samples 
} in. in diameter were cut from large 0.016 in. aluminum sheets that had been 
uniformly coated with several layers of lacquer (cellulose nitrate dissolved in 
amyl acetate) containing plutonium nitrate dissolved in acetone. Experience 
with this technique indicated that the uniformity was better than LO%, so 
that with a beam uniformity of better than 1°¢, the error introduced by non- 
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uniformities was less than 0.1%. The samples were mounted on either side of 
a thin aluminum partition in the common cathode of a back-to-back gridded 
ion chamber. For the fission counting the chamber was operated at low pressure 
to permit complete discrimination between fission and @ pulses. The central 
part of the common cathode could be rotated to allow each source to be 
counted for an equal time in each side of the chamber; this eliminated the effects 
of neutron-beam attenuation and of any difference in the settings of the 
thresholds for registering pulses. For the a counting, the neutron beam was 
shut off, the chamber pressure doubled, and the amplifier gain increased by 
a factor of 16. The @ particles were also counted for equal times in each side 
of the chamber to eliminate differences in threshold setting. There was no 
handling of the sources between fission and a counting and therefore no change 
in counting geometries. The samples were surrounded by a simple ring col- 
limator of 0.75 in. internal diameter and 0.020 in. high. The threshold bias was 
set so that no pulses were recorded from oblique tracks terminated by the ring. 
This procedure does not entirely eliminate the back-scattering from the sample 
backing, but since the a energies are closely the same for the two isotopes being 
compared, the differential effects are extremely small. It does, however, reduce 
the effects of small irregularities in the surface of the samples. 

Experience has indicated that, over the range of sample thicknesses used 
here, measurements of the ratio of fission to a rates of two samples can be 
made with a reproducibility of about +0.2%. 


TABLE I 
ISOTOPIC COMPOSITION OF SAMPLES 


Atom % abundance 


Mass No.  Pu*#® sample Pu®® sample 
238 — — 
239 10.131+0.025 97.884+0.011 
240 89.365+0.026 2.043+0.010 


241 0.483+0.006 0.073+0.005 
242 0.018+0.005 = 





The isotopic compositions of the samples, measured by M. Lounsbury of this 
laboratory, are given in Table I. The amount of Pu** was too small to be 
detected by the mass-spectrometer analysis. The contribution of Pu*** to the 
a activity of the samples was measured directly by a-pulse analysis, using a 
gridded ion chamber similar to that described by Harvey, Jackson, Eastwood, 
and Hanna (1957). It was about 1.5% of the total @ activity of both samples. 
This procedure also corrects for the presence of any Am*' @ activity (from the 
B decay of Pu*'), which cannot be resolved from that of Pu’. The contribu- 
tions of Am*!, Pu?’8, and Pu®” to the fission rates and of Pu*! and Pu?” to 
the a rates are negligible and therefore have not been included in the following 
analysis. 


Analysis 
Since each sample was counted in both sides of the chamber, the over-all 
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counting efficiency and neutron flux were the same for each sample and the 
ratio of the fission rates is 


F’/F = (N’/N)(Xin's/dino), 


where the primed and unprimed symbols refer to the Pu? and Pu**® samples 
respectively, NV is the total number of atoms in the sample, and }ono = 
Ngogt+Nooo +01, Where a9, oo, 0; are the fission cross sections and M9, mo, m the 
atom percentages of Pu’, Pu®#°, and Pu*! (given in Table I). 

The ratio of a counting rates is 


a’/a = (N’/N)(Sn'd/T nd), 


where Sond = nsdAgtMgrAg+MoAo, and Xg, Ag, Ao are the a-decay constants and 
Ns, Ng, No the atom percentages of Pu**8, Pu*89, and Pu**®. 

The quantity A = (m9/m9’)(32n’A/>donX) was determined from the measured 
Pu’ contribution and the value given by Inghram et a/. (1951) for o/s, 
which is (24,400 years/6580 years) +0.6%. 

If R = (F'/F)/(a’/a), mn = o1/09, No = no/ng, and Ny = m/m, the ratio 
of the fission cross sections of Pu*#® and Pu? is given by 


oo/oy = {AR(1+nN,) — (1+) }/(No' — ARN}). 
















The value of ~ is sensitive to the energy spectrum of the neutron beam. 
Therefore the value used, 7, = 1.25+0.02, was determined in the same 
neutron beam by the method used here but with a sample containing about 
50% Pu™!, 










Cadmium Ratio 

This beam came from the graphite reflector of NRX, 9 in. from the calandria. 
The temperature of the thermal component, 180+20° C., was derived from a 
gold transmission measurement, and the epithermal neutron density, (1.0+ 
0.1)%, from a cadmium transmission measurement, both with a boron detector. 
The correction for fissions produced in the Pu**° of the Pu**° sample by neutrons 
with energies greater than 0.5 ev. was found by measuring the cadmium 
transmission of the Pu”® and Pu*!® samples. These were 0.44% and 0.52% 
respectively, which gives a correction of 0.08+0.01% to R, equivalent to 
0.09+0.01 barn in apo. 
















Checks for Contamination 

Because of the small value of oo, contamination of the samples by a emitters 
or fissile materials can be serious. The important case would be if contamination 
caused an underestimate of oo; this could happen if the Pu*? sample was 
contaminated with an a@ activity or the Pu*® by a fissile substance having a 
low specific activity. The a spectra were examined in the gridded ion chamber 
mentioned above. Except for the Pu**’ group already corrected for, there were 
no a groups larger than 0.1% of the total « activity outside the energy interval 
4.90 to 5.24, which was masked by the main group. Pa”*! is the only possible 
contaminant emitting an @ group in this energy interval which is not accom- 
panied by higher-energy a groups of similar intensity from daughter products 
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In this case the decay chain is held up by the 22 year Ac”’. However, since 
the sample had not been purified for at least 12 months before the measure- 
ment, accumulated daughter products would produce higher-energy @ groups 
totaling about 1/10 of the Pa**! @ activity. A search for these a groups was 
made using a low geometry chamber (Cruikshank, Littler, and Ward 1948) 
at a pressure sufficient to stop all except the Pu*** a particles outside the count- 
ing volume. This showed that Pa**! contamination was less than 0.02%. 

The only fission contaminant that would not have been detected by the 
mass-spectrometer or a-pulse analyses is U?**, Chemical analyses showed that 
the Pu**® sample contained about 0.8% uranium. In order to show that this 
uranium was not enriched in U?*, the uranium was separated from most of the 
Pu*** and counting samples were prepared from both the uranium and pluto- 
nium fractions. The fission/a ratios of these samples were compared to that 
of the original sample and were the same except for a redistribution of the 
Am! during the Chemical processing. This shows that the U?* content of the 
Pu’? sample was negligible. 


RESULTS 


The values of oo/o9 obtained are given in Table II]. Three measurements 
were made using two pairs of samples referred to as A and B. The internal 














TABLE II 
THE RATIO OF THE THERMAL FISSION CROSS SECTIONS OF Pu? AND Pu?3? 
Samples R = (F’/F)/(a’/«) Standard dev. (%) (o7(40) /os(49)) X 104 
A 0.03351 0.16 —9 341.9* 
A 0.03356 0.2 —7.642.3 
B 0.03345 0.2 —11.542.3 
Average 0.03351 0.12 —9. 341.4 





“The errors quoted are statistical errors from counting only. 


consistency of the results is illustrated by showing only the statistical errors 
arising from the measurement of RX. There are additional errors of +7.1 X 104, 
+3.1X10-', and +0.6X 10~ from errors in A» As, Mass-spectrometer analyses, 


and 01/09 respectively, and a mean value 
oo/o9 = (—9.3847.7) X 107 


is obtained. For this neutron beam, og = 900+40 barns (Hughes and Harvey 
1955), which gives for the thermal fission cross section of Pu**° 


oo = —0.8+0.7 barn. 


This supports the conclusion that the thermal-neutron fission cross section 
of Pu*’ is much smaller than the 4.4+0.5 barns reported by Hulet et a/. (1956). 
The average resonance spacing for neutron capture in Pu™? is about 11 ev. 
(Fluharty 1957), which means that the thermal cross section is mostly from 
the wing of the 1 ev. resonance. That is, the fission-to-capture ratio would be 
the same at thermal energies as at the resonance. The 0.05 barn value for the 
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thermal fission cross section (Leonard et a/. 1956), based on this assumption 
and a measurement at the 1 ev. resonance, is probably the best estimate 
available. It is hoped that further measurements with more nearly pure Pu**° 
samples will give a direct value.* 
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SOME OBSERVATIONS OF GEOMAGNETIC 
MICROPULSATIONS! 


H. J. Durrus AND J. A. SHAND 


ABSTRACT 


New observations are reported of the diurnal variation, direction, frequency 
spectrum, and geographical distribution of geomagnetic micropulsations (Pc 
and Pt). It is argued from a comparison of published data that these phenomena 
have primarily a solar time dependence, and may occur locally. At our two 
stations, separated by 60 longitude degrees, related signals are only occasionally 
observed. The phase characteristics of such related signals, together with differ- 
ences in directional characteristics reported by other observatories, put serious 
strain on recent theories of the outer atmospheric origin of micropulsations. 


A. INTRODUCTION 

kor over half a century observers have described and conjectured upon the 
origin of geomagnetic micropulsations—the name given to the as-yet-un- 
explained oscillations of the earth’s magnetic field having amplitudes of 
only a few gammas, and periods from a fraction of a second up to more than 
ten minutes. Fig. 1 shows a series of examples taken from a single day’s record. 

Three types of micropulsation are recognized. Typical of day are regular, 
long wave-trains of up to about one minute period having an average amplitude 
of one gamma or less. These, called Pc’s, are illustrated by the top two 
sections in Fig. 1. The most recent work suggests that they have a 27 day 
recurrence tendency, and that there is a clear correspondence between their 
occurrence and the appearance of Babcock’s unipolar magnetic (UM) regions 
on the sun (Kato and Akasofu 1956). Angenheister (1954) has shown that 
their diurnal variation is similar to the variation of the magnetic K in- 
dex. Typical of night are well-damped longer-period oscillations, usually of 
greater amplitude, called Pt’s. A sample is illustrated in the fourth section 
of Fig. 1. These too have a 27 day recurrence tendency, and their occurrence 
follows the K index. Jacobs and Obayashi (1957) have recently shown that 
micropulsations of long period (probably Pt’s) have periods which increase 
with magnetic latitude. It is predicted by the theory that they are caused 
by toroidal oscillations of the outer atmosphere. In spite of their preferred 
times of occurrence, both Pc’s and Pt’s can occur at any time of day either 
separately or together and in great variety. A third type, distinguished by 
large amplitude and regularity, is called the Rolf or Giant micropulsation. 
Because it is recognized only in or near the auroral zone, we have not recorded 
it. 

The electromagnetic background consisting of oscillations of slightly 
higher frequency (0.1-10 c.p.s.) has received scant attention, being very 
weak (0.01-1 y) and difficult to record. It is the subject of our current work 
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Fic. 1. Some examples of geomagnetic micropulsations recorded during a single day. 


and will be reported in a separate paper. Characteristic trains of oscillations 
“en perles” as described by Troyigkaya (1957) occur in this frequency range 
together with the lowest frequency components of lightning discharges and 
several other as-yet-unexplained types of disturbance. 

At higher frequencies Goldberg (1956) has reported measurements of the 
background in the 1-150 c.p.s. range, Holzer et al. (1957) in the 10-900 c.p.s. 
range, and Aarons (1956) in the 37-30,000 c.p.s. range. These workers were 
concerned principally with the part of the electromagnetic background 
arising from lightning discharges. However, Holzer and Aarons independently 
report a peak in the spectrum at 30 to 40 c.p.s., which Aarons attributes to 
the gyromagnetic frequency of ionized sodium. A second peak in the region 
of 8-10 kc./s., reported by Holzer and other workers, is attributed to the 
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energy spectrum of lightning discharges. It is apparent, therefore, that the 
spectrum of the electromagnetic background in the range 0.001 to 10,000 
c.p.s. is very complicated, encompassing the gradation from large geomagnetic 
disturbances to atmospheric electrical discharges together with the enormous 
man-made background of switching transients, power frequencies, and their 
harmonics. 

This paper reports observations on geomagnetic micropulsations in the 
range 0.001—0.1 c.p.s. There is considerable conflict in the literature about 
even the limited frequency band we report here. Some authors maintain that 
micropulsations occur simultaneously on a world-wide basis, while others 
find no coherence beyond a 300 km. radius. Some observers find most activity 
to occur about local midnight, while others report it to occur at noon. One 
of the main purposes of this paper is to examine the probable causes of these 
discrepancies. 

It is a natural result of the conflicting data that there should be a number 
of conflicting theories of the origin of these regular oscillations. St¢rmer’s 
early theory that micropulsations are caused by particles of solar origin 
following periodic orbits in the earth’s magnetic field has recently been 
revitalized by Bennett and Hulburt’s (1954) discussion of magnetic focusing 
of solar ion streams. Holmberg (1953) has suggested Alfvén waves in the 
lower ionosphere as a source, and Dungey (1954) has examined torsional 
oscillations of the outer atmosphere as a possible cause. Several Japanese 
authors—Kato and Akasofu (1956); Kato and Watanabe (1957); Watanabe 
(1956)—propose that Rolf micropulsations are caused by Dungey’s torsional 
oscillations. However, they conclude that the ordinary micropulsations 
(which are observed simultaneously at widely separated stations) are caused 
by poloidal oscillations of the outer atmosphere. They attribute the steady 
flux of daytime oscillations (Pc’s) to the incidence of turbulence and shock 
waves in the solar corona upon the sunlit side of the ionosphere. 

Our observations confirm the distinction made by other observers of the 
two regimes of oscillatory activity, Pc and Pt. We do not often find either 
type to be continent-wide in extent, and even when they are, we do not 
find them to have the same phase. In short, then, these observations do not 
support the “poloidal oscillation” theory in its present form. As no Rolf 
pulsations were recognized in the data reported here, and no simultaneous 
data from a different magnetic latitude are included, the torsional oscillation 
theory of Dungey cannot be tested. 

Our data reveal other interesting differences from what has already been 
reported. The directional characteristics of the micropulsation disturbance 
vector are found to be different from those reported earlier from Japan. 
The spectral distribution is found to be more like the continuous spectrum 
recently reported by Angenheister (1954) than like the several “‘line’’ spectra 
reported by earlier workers using equipment less capable of covering a wide 
frequency range. As a result of these and other known properties of micro- 
pulsations, it is possible to explain some apparent contradictions in the 
literature. 
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B. EXPERIMENTAL TECHNIQUE 
The magnitude and direction of micropulsations in the frequency range 
0.001 to 0.1 c.p.s. have been measured over about 150 days from 1954 to 
1957 at Victoria, B.C. Simultaneous measurements of the vertical component 
only were made on several occasions, using identical forms of equipment, 
at both Victoria and Halifax, N.S. 


(a) Detection 

The detecting equipment at Victoria consists of an array of three mutually 
perpendicular air-core induction coils mounted in a rigid wooden octahedron. 
Each coil consists of 1100 turns of No. 21 Formel-covered wire, and is 20 
feet in diameter, thermally insulated, and electrostatically shielded. (These 
large coils were chosen because of their low self-inductance, as work at higher 
frequencies was contemplated.) Deflection of the members of the octahedron 
in a gusty wind is an inconvenience because any movement of the coils in the 
earth’s field induces large noise signals. These, however, are generally of 
higher frequency than we were recording. Towards the end of the experiments 
a number of mu-metal-cored coils became available and it was possible to 
check the directions obtained on the air-core coil system with measurements 
made on a differently oriented mu-metal-core system. No difference was 
found between the results of the two detector systems, and we conclude 
that any small magnetic anomalies in the immediate vicinity of the detectors 
have not influenced our observations of directions. 

The signals from the coils are fed by buried, lead-shielded cables to d-c. 
chopper amplifiers driving Esterline-Angus pen recorders. Each amplifier 
incorporates an integrating circuit of 1-minute time constant. A sinusoidal 
magnetic field of known amplitude may be applied to all of the detector coils 
by means of a single calibration coil. The over-all system response of one 
of the detector-amplifier-recorder combinations is shown in Fig. 2. Although 
each system is subject to d-c. drift this is very slow compared to the periods 
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which we measure and is non-sinusoidal. Drift corrections were noted on 
the chart as they were manually applied. 

The detecting equipment at Halifax consists only of a coil with its axis 
vertical. Otherwise the equipment is the same. 
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(6) Analysts 

The paper-chart records were analyzed by visual inspection, counting as 
an oscillation any event with three or more equally spaced axis crossings. 
This criterion may be contrasted with that of Holmberg (1953), who took 
the average of five cycles. Fluctuations consisting of only a departure and 
return to the axis of the record as treated by Chernosky et al. (1954) were 
not included in the analysis. Although our criterion would pass over a wave 
beating with its second harmonic if it were of short duration, it is unlikely 
that the situation would go unrecognized in a wave train of any length. 
Most of the errors from this source will be for the long-period, damped Pt'’s. 

As there is a great danger that personal bias may slant the results of visual 
inspection, an electronic analysis was also carried out. The data were recorded 
on slow-speed (0.003 in./second) amplitude-modulated magnetic tape, re- 
recorded several times until a speed ratio of up to 80,000:1 was achieved, 
and then analyzed on a Panoramic Sonic Analyzer. Although good qualitative 
agreement was obtained between the two analyses, very few data in the 
frequency range discussed in this paper were treated by magnetic tape, 
as the technique was not sufficiently flexible for the analysis of other para- 
meters. Instead, the data from visual inspection were analyzed both flexibly 
and inexpensively by punching them on approximately 5000 I.B.M. cards 
which were subsequently programmed to give the histograms shown in the 
accompanying figures. 

A greatly improved tape recording facility is now available and we propose 
to make a further electronic analysis when sufficient data have been recorded. 


C. RESULTS 

Unfortunately we do not have an adequate theory of the origin of micro- 
pulsations. In searching for clues about their origin it seems reasonable to 
ask whether any particular frequency has an unusually large amplitude 
which might be caused by a resonance phenomenon, or whether the spectrum 
shows a smooth frequency dependence which might be caused by the con- 
tinuous change of a single variable such as ionospheric viscosity. An idea 
ot the stability and damping of the processes involved can be gained by 
examining how long the wave trains last and how this depends on their 
frequency. The number of times wave trains of given frequency are initiated 
during each time of day will suggest whether the origin of that frequency 
is associated with solar effects on the local ionosphere, or on the ionosphere as 
a whole. The dependence of the direction of the disturbance vector upon 
frequency and time of day, as well as its general behavior, should give an 
indication of the required position and extent of the source. A comparison 
of simultaneous records from two stations will, of course, reveal much more 
about conditions which a postulated source must fulfill. 

The following properties of micropulsations have therefore been presented: 
(i) frequency spectrum, (ii) probability of occurrence, (iii) diurnal variation, 
(iv) direction, (v) geographical distribution of phase and amplitude. These 
properties will be described separately and compared with the literature 


before an interpretation is made. 
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(1) The Frequency Spectrum (Fig. 3) 

Fig. 3 shows the relationship between average amplitude and the period, 
disregarding other factors such as time of day or duration of wave train. 
It was obtained from only the vertical component of the fluctuating magnetic 
field because part-time recording of this extended for some years prior to the 
recording of other components. 
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There is no enhancement of amplitude about a period of 0.4 to 0.6 minute. 
This is noteworthy considering Holmberg’s (1953) data for the U.K., and 
Chernosky’s et al. (1954) data for Tucson, Arizona. However, the Tucson 
data only cover periods from 0.1 to 0.5 minute, and the U.K. data were 
taken on an instrument having a response which was strongly peaked at a 
period of 0.2 minute and fell off sharply by 2.5 minutes. Both these sets of 
data may therefore tend to focus attention on the 0.3 to 0.5 minute range. 
Most of Holmberg’s data and all the Tucson data are for the vertical field 
component. 

The second feature of Fig. 3 is a much stronger peak in the spectrum at 
periods of 0.8 to 1.0 minute. This is outside the range of the Tucson measure- 
ments, but agrees very well with the U.K. results and with the reports of 
other observers, who, unfortunately, give no quantitative data, e.g., Terada 
(1917) in Japan. 

The third feature is a pair of peaks (probably representing only a single 
flat maximum) in the range 1.5 to 2 minutes period. This feature is also 
reported by Terada (1917) and agrees with the data of Holmberg. 

The general increase in amplitude at long periods is reported by 
Angenheister (1954), but his points are too widely spaced at shorter periods 
to resolve the peaks discussed above. 

It is fair to say that when the frequency response of their apparatus is taken 
into account, all authors report the features shown in the spectrum of Fig. 3. 

These spectral peaks at 1, 2, 3 minutes appear statistically to be har- 
monically related. Different observers agree that several frequencies may 
occur simultaneously but seldom as exact harmonics. It also appears that the 
wider the bandwidth of the detecting equipment, the less striking are the 
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spectral peaks. Their reality, however, is suggested by the similar distribu- 
tions found by different observers. 


(it) The Probability of Occurrence 

The probability of occurrence of micropulsations depends upon both 
frequency (thus, in the light of Fig. 3, upon amplitude) and upon the time 
of day. The probability of detecting a micropulsation at a given instant 
depends upon its probable duration and the probability of initiation. Fig. 4 
shows how the duration of the wave trains of the vertical component varies 
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with frequency. The resemblance to the variation of amplitude with period, 
Fig. 3, is very marked except in the case of periods of 0.1 to 0.5 minute. 
Thus, large-amplitude wave trains last a long time, but very-long-period 
oscillations, which generally occur at night, are damped out faster than the 
shorter-period Pe’s even when averaged over all recorded data. 

The average duration of wave trains having periods from 0.2 to 0.4 minute 
is less than those in the 0.1 to 0.2 minute band although their total recorded 
duration is larger. This is because, although all these frequencies tend to 
persist, several often occur at once resulting in an interference pattern. Nulls 
in the envelope may be included as separate initiations, and the average 
value of the duration, obtained by dividing the total duration recorded 
by the number of initiations, is therefore small. The impression gained from 
the original record is that the 0.2 to 0.4 minute band is very often present, 
as indeed is shown by its large total duration. 

Both total and average durations are large in the 0.8 to 1.0 minute band. 
The average duration in the 1.8 to 2.0 minute band is also unusually high, 
but its reality is doubtful as it depends on only 39 separate measurements. 
Nonetheless it corresponds to a peak in average amplitude in Fig. 3, and has 
some support from Terada’s (1917) work. 

The peak at 2.5 to 3.0 minutes is suspect again for lack of data, and is 
unsupported by any other evidence. It is to be expected that the average 
duration will increase with period because durations of less than one cycle 
do not occur in our analysis. 
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In order to gain a complete picture of the probability of occurrence, we 
must consider the probability of initiation of wave trains as well as how long 
they lasted. The total number of initiations of wave trains is given in the last 
column of Table I. In the absence of a recognizable distribution function, 
one does not know how to divide up the spectrum into appropriate bandwidths 
for these numbers to be meaningful. Hence the rather large counts at the long- 
period end of the spectrum must be judged in the light of their great band- 
width. The results would look quite different if plotted in terms of a bandwidth 
of 1 c.p.s. as is often done. Nevertheless, there is marked similarity between 
the average duration and the probability of initiation. We have already 
distinguished between the two types of micropulsations, and these data 
substantiate the picture of occasional, damped, long-period wave trains 
which appear at night, and the continuous, long-duration, short-period 
oscillations which appear in large numbers during the day. 

The product ‘average amplitude X average duration”’ is related to the 
electromagnetic energy of the micropulsations. We term it the “‘activity’’. 
It is plotted in Fig. 5, and displays quite prominently the preferred periods 
0.8 to 1.0 minute and 1.8 to 2.0 minutes, but shows only a very slight maximum 
from 0.2 to 0.6 minute. 
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(211) Diurnal Variation 

We shall now examine in detail the daily course of micropulsations. Table 
I shows how the number of wave trains initiated in Victoria at each hour of 
the day depends upon the frequency. It shows that the most probable time 
for the initiation of any frequency is around noon, but that the lowest 
frequencies have a high probability of initiation about midnight. Fig. 6 
shows the diurnal variation of activity (not occurrence) covering the whole 
frequency band as deduced from our observations of the vertical component 
only of the fluctuating field. 

It was apparent from large variations between half-hourly means that 
Victoria data are far too few for rigorous treatment. Nevertheless, certain 
trends may be noted. 

The most important feature of Fig. 6 is that there is a maximum at noon 
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and at midnight. The very large peak at about 1800 is contributed by a 
relatively few long-period large-amplitude oscillations and has such prominence 
because of our definition of activity. Indeed the quietest time of day is normally 
between 1800 and 2200 hours. 
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The peak centered on 1200 is contributed principally by the 1 to 3 minute 
band, whereas the activity having less than 1 minute period appears to peak 
at 1000 and 1300 with a slight minimum at noon. 

As has been noted, there are conflicting reports in the literature about 
the diurnal variation. Any reasonable theory of the role of the sun and 
ionosphere must be based upon a firm knowledge of the diurnal variation. 
Hence we shall try to clarify the picture of the world-wide diurnal variation 
of activity. Fig. 7 shows most of the published data. 

The curves of van Bemmelen (1903) and Lubiger (1935), Fig. 7(d), show a 
maximum probability of occurrence toward midnight. There are several quota- 
tions from van Bemmelen’s work to be found in the literature to the effect that 
the frequency of occurrence is greatest at midnight and also that the frequency 
is highest (period least) at midnight. This is correct only when referring to 
his observations taken from 9 p.m. to 6 a.m., local time. He was prevented 
from making daytime observations at the time of one of his reports by inter- 
ference from electric railways. In general he observed the frequency of 
oscillation to be highest during the day. Nevertheless, the greatest number 
of micropulsations were detected near midnight, and we shall return to this 
point. 

Troyi¢gkaya’s (1955) curves were originally plotted on G.M.T. and purport 
to apply to the whole world, but in Fig. 7(e) we have replotted them on a 
local time scale appropriate to Alma Ata (77° E.). Only Angenheister’s (1954) 
and Terada’s (1917) results, Fig. 7(@) and (0), are plotted in frequency bands, 
and it will be noted that different frequencies have quite different diurnal 
variations. 

Before discussing Fig. 7 in further detail, it will be worth while to examine 
the directional properties of the micropulsation disturbance vector. We 
shall see that, between the frequency dependence of the diurnal variation 
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refer to the months in which the observations were taken. 


(b) Angenheister’s (1954) data for Gottingen, 1952-53. The axes are similar to those of (a). 
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Fic. 7. (c) Holmberg’s (1953) data for the U.K., 1932-33; and Chernosky, Maple, and 
Coon’s (1954) data for Tucson, Arizona, 1947. 

(d) Van Bemmelen’s (1903) and Lubiger’s (1935) data. 

(e) Troyigkaya’s (1955) and Kalashnikov’s (1955) data from earth current (curve A) and 
fluxmeter (curve B) measurements at Alma Ata, U.S.S.R., 1951. 

(f) Lazarenko’s (1950) data from magnetometers, U.S.S.R., 1937-48. The Zi-Ka-Wei data 
may have been taken from the same source as (d). 


and the directional properties, it is possible to rationalize the apparently 
unrelated curves of Fig. 7. In addition, however, we shall have to introduce 
the conclusions of Vestine et al. (1947) that the amplitude of micropulsations 
is strongly latitude dependent. 


(iv) Direction 

In the total time during which directional observations were made, only 
a small proportion of the record was considered suitable for measurement. 
The reasons for rejecting most of the record are: 

1. During quiet periods the signal amplitude is too low for reliable measure- 
ment. 

2. During active periods the signal may be too complex to measure by 
visual means. 

3. Strong winds disturbed the record. 

4. Most of the time when moderate and readable signals are obtained, 
phase differences occur between two or more of the directional components. 
It is most important to note that our data on directions apply only to 
oscillations which are in phase in all three components. Only when the 
components are in phase is the locus of the vector in time a line in space 


for which unique angles can be defined. 
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In spite of these limitations on the selection of data, we are confident 
that the directions plotted are representative. Whatever the level of activity 
or the phase relation between components, we consistently operate one 
detector system at one-third the sensitivity of the other two orthogonal 
detectors. Thus the limits of directional variability remain about the same 
as the limits shown by the scatter diagrams (Figs. 8, 9, 10). Exceptions 
to this general rule are occasionally observed, especially during magnetic 
storms when for a short time a considerable variation in direction can be 
expected. Another notable exception occurs when steeply inclined short- 
period signals are received. 
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The signals whose directional components are not in phase show some 
regularity in the following respect. The angular difference between East-West 
and vertical signals is frequently about 7/4, but is usually less between the 
North-South and vertical signals. It is interesting to compare this with 
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Terada’s (1917, p. 60) statement ‘‘... van Bemmelen found AZ/AX (the phase 
lag) in Batavia invariably insignificant. On the other hand (in Japan) the 
most regular waves with periods less than about 4 minutes show as a rule 
remarkable?phase retardation of the vertical component after the horizontal 
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Fic. 10. (a) 8 to 35 seconds period. (b) 40 to 155 seconds period. (c) 160 to 1000 seconds 
period. 


ones, generally greater than 2/4’’. Our observations do not agree with this 
description of events in Japan, but may be like van Bemmelen’s observation. 
That is to say, in Japan the horizontal components are in phase and the 
vertical component lags them. In Victoria the North-South and vertical 
component are in phase, and the East-West component leads or lags. 

The plots of inclination and azimuth show that the range of azimuth exceeds 
the range of inclination. From Fig. 10 and Table II there is a marked trend 
for smaller inclinations to be associated with smaller azimuths and greater 
inclinations with greater azimuths. (The mean does not lie in either the 
geographic or the magn: +: meridian.) Because short-period signals are 
most frequently encountered between 1000 and 1600 hours local time, the 
maximum ranges of these angles are observed then. We have discovered 
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no other diurnal variation of either inclination or azimuth. It will be noted 
from Figs. 8 and 9 that there is a considerable spread in azimuth and dip 
for the higher frequencies and less for the lower frequencies. 


























TABLE II 
Azimuth in degrees true Inclination in degrees below the horizon 
Periods Mean Standard Number of Mean Standard Number of 
in sec. azimuth deviation observations inclination deviation observations 
O- 20 324.3 26 3l 55.7 22 30 
20— 40 aléia 35 69 47.9 19 69 
40— 60 336.6 31 19 39.7 16 19 
60— 80 338.4 27 31 S08 12 30 
80-100 335.5 22 29 34.3 10 29 
100-120 344.7 12 14 oie 9 16 
120-140 ga0.2 8 eH 35.8 10 26 
140-160 320.6 38 8 30.5 ll yg 
160-180 337.5 19 t 
180-200 332.5 22 8 38 y 10 
200-220 312.8 26 9 a1 i7 12 g 
220-240 
240-260 320.0 24 14 39 14 15 
260-280 
280-300 
300-320 311.3 20 15 31.3 ll 16 
320-340 
340-360 


360-380 317.9 26 12 33 12 10 


The standard deviations are shown for the inclinations and azimuths in 
Table II. (It is apparent from Figs. 8 and 9 that the distributions are not 
Gaussian.) These deviations show that the spread of inclinations decreases 
rapidly with increasing period, but the spread of azimuths is more nearly 
independent of period. 

Terada’s observations of the directions of micropulsation vectors are the 
only ones with which we are familiar, although a description of the direction 
of a few Rolf micropulsation vectors has been published. Terada (1917) 
shows only the average dip as a function of frequency, without indicating 
the spread in observations. However, if the Victoria data are treated in this 
way in Table II, a qualitative difference emerges, viz., in Victoria the lower 
the frequency the less the mean angle of dip—-in Japan the lower the 
frequency the greater the angle of dip. Terada’s results were obtained from 
data for the entire year of 1913. The Victoria results were obtained from 
only about six months’ sampling, during which, however, no change in 
direction was noticed. 

An unknown but possibly very important factor in this study is the effect 
of earth currents. Bondarenko (1953) has shown experimentally that some 
earth currents are induced by the vertical magnetic field changes, and 
Cagniard (1956) discusses the relationship between the horizontal electric 
and magnetic components given a homogeneous, layered earth. 

The difference in the inclination vs. frequency curves for Victoria and 
Japan is difficult to understand on the basis of earth currents, since these 
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would be expected to produce the same sense of frequency dependence of 
inclination at both places, although the rate of change of inclination with 
frequency might be different because of different earth conductivities. 

The work of Burkhart (1955) shows that for high frequencies there is a 
large inductive reactance as well as the inherent ohmic resistance of the 
induced earth current circuit. This means that the secondary magnetic field 
from the earth current system will be nearly 180 degrees out of phase with 
the primary inducing field, and will give rise to the sort of in-phase components 
we deal with. It would also imply that the vertical components of short-period 
oscillations would probably be relatively smaller than those of long-period 
oscillations, which is not what we observe. It explains Terada’s observations 
of the dependence of the inclination on frequency, but does not explain why 
the phase lag of the Z component behind the X component approaches zero 
at long periods. 

This bears upon possible explanations for the differences in diurnal variation 
in Fig. 7. First of all, Terada’s (1917) curves in Fig. 7(a@) show that below 80 
seconds’ period there is a peak of activity near local noon. From 80 to 200 
seconds there is less variation, but for longer periods the maximum appears 
toward midnight. 

Terada is reporting the horizontal field fluctuations, and his short-period 
oscillations are more nearly horizontal. Thus despite their small amplitude 
he will have a good chance of detecting them toward solar noon when all 
observers (including Troyig¢kaya) agree that the periods of the oscillations 
are shortest. The longer-period oscillations usually appear toward midnight, 
and although they have a smaller horizontal component, yet they have an 
easily detectable amplitude in conformity with the general rise in amplitude 
with period in the spectrum of micropulsations shown in Fig. 3. 

In Fig. 7(c) neither the U.K. results of Holmberg (1953) (10-100 sec.) 
nor the Tucson measurements of Chernosky et al. (1954) (5-30 sec.) extend 
to a sufficiently long period to display the high probability of occurrence of 
micropulsations at midnight. 

Our own results (Fig. 6) tend to emphasize the noontime contribution 
of the higher frequencies because their direction favors our vertical induction 
detector and the low-frequency, high-amplitude, horizontal oscillations are 
discriminated against. 

Returning to the older work of van Bemmelen (1903) and Lubiger (1935), 
Fig. 7(d), we note that their results deal with oscillations of period pre- 
dominantly over 90 seconds, and the distribution resembles Terada’s curve 
for this sort of period. 

On these grounds then the curves of Fig. 7(d) for Samoa, Batavia, 
Buitenzorg, and Zi-Ka-Wei may be compatible with those for Victoria, 
Tucson, the U.K., and Japan. The curve for Potsdam, however, is still 
unexplained, as are the curves of Troyickaya (1955), Kalashnikov (1955), 
and Lazarenko (1950), of the U.S.S.R. 

The curves of Troyickaya (1955) are for earth current measurements, 
which Bondarenko (1953) shows are at least partly due to induction from 
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the changes in the vertical component of the earth’s magnetic field; those of 
Kalashnikov are from the measurements of the vertical component by means 
of a fluxmeter. Unfortunately Lazarenko (1950) does not state to which 
component his curves apply. At any rate, his data (plotted to local time) 
do not bear out Troyi¢kaya’s contention that the maxima in Fig. 7(e) are 
related to G.M.T., and, in particular, to the transit of the sun across the 
North and South geomagnetic poles. Indeed, the differences between the 
diurnal variation curves of Terada (1917) for different years suggest that it 
would be most difficult to assign such close time limits to the maxima, even 
if a suitable frequency band could be found. 

From Lazarenko’s description of the observations at Keles, Slutsk, and 
Uelen, one may guess at a reason for the apparent discrepancy in Fig. 7(/). 
One should note first of all that it is the high-latitude stations which show 
a maximum about local noon. Furthermore Lazarenko states that 20,000 
hours having micropulsations were observed during 10 years at Slutsk but 
that a mere 2000 such hours were observed at Keles during 13 years. As 
Vestine et al. (1947) showed the micropulsation amplitudes to be greatly 
enhanced in auroral zone latitudes, the foregoing suggests that many of the 
short-period (sic, small-amplitude) oscillations were below the threshold of 
the Keles magnetometer (1 gamma/mm.). As we have seen, the Zi-Ka-Wei 
data‘ of Fig. 7(d), reproduced by Lazarenko, pertained to periods greater 
than 90 seconds. 

Finally, in examining Lazarenko’s curve from Keles, which does not show 
a maximum of occurrence about noon, we should note the geographical 
proximity of Keles to Alma Ata, where at least some of Troyicgkaya’s data 
for Fig. 7(e) showing a maximum at noon were obtained. One must conclude 
that different frequency bands, different directional components, or different 
selection criteria were used. 

Troyigkaya (1955) shows two curves, distinguishing two classes of micro- 
pulsations in the bands 10-40 seconds and 50-60 seconds period, respectively. 
It is disturbing that periods as short as the latter produce a maximum about 
midnight, for, as is easily seen from Table I, our results do not show such a 
maximum. However, this may be due to the increase of period with latitude 
noted by Jacobs and Obayashi (1957). Combining Troyi¢gkaya’s two sets of 
data gives the curve shown in Fig. 7(e) with maxima about noon and mid- 
night. The enhancement of the noon maximum due to higher frequencies 
for the fluxmeter data (B) suggests that, as in Victoria, these high frequencies 
are more nearly vertical in direction than the long periods, and hence are 
detected preferentially by induction coils having their axes vertical. We 
cannot properly comment upon the properties of the earth currents, but 
one suspects from Bondarenko and Cagniard that they are induced by Pc’s 
in the day, and Pt’s at night, and that they have the properties reported 
by other observers; i.e., are not dependent upon G.M.T. 

All the diurnal variation curves except that for Potsdam have now been 
assessed, and possible explanations offered for their differences. The only 
reasons one can offer for the lack of daytime activity at Potsdam, lying 
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as it does at a comparable magnetic latitude to Slutsk and Eskdalemuir, 
are either (1) that only one component was measured, and that this com- 
ponent was deficient in higher-frequency oscillations, e.g., horizontal, (2) that 
different selection criteria were applied, or (3) that only long periods were 
detectable because of low sensitivity or slow paper speed. The third is not 
unreasonable in view of the early date of this report. Furthermore, Angen- 
heister’s (1954) data for the horizontal component should be applicable 
for Potsdam too. 

All these results bear out the supposition of a local time dependence of 
frequency of occurrence, amplitude, duration, and period. Troyi¢kaya (1955), 
however, states that initiation is simultaneous over the whole of the U.S.S.R. 
and is, indeed, world-wide, and that its probability of occurrence depends 
upon universal time. Yet Holmberg (1953) found no coherence beyond a 
radius of 300 km. 


(v) Geographical Extent 

In order to examine more closely this question of the geographical extent 
of micropulsation disturbances, simultaneous observations of the vertical 
component were made at Victoria and Halifax using equipment identical 
in design. In general there is little correlation between the records at the 
two stations. Micropulsations of different frequencies may occur at each 
station, or micropulsations of large amplitude may occur at one station 
when nothing is detected at the other. However, there are a number of 
occasions upon which micropulsations of the same period occur at both 
stations. When this happens we observe a very marked tendency for the 
two signals to occur either exactly in phase or exactly out of phase, although 
there are times when the signals seem to jump discontinuously from in-phase 
to out-of-phase in a random fashion. In order to substantiate these obser- 
vations the phase was plotted as a function of the time of day, and a very 
strong indication was found of a diurnal variation. More data are required 
to show this conclusively, however, and it is mentioned here as only a tentative 
result pointing out a direction for future work. The implications of this 
conclusion are very important in justifying the poloidal oscillation theory 
of Kato and Akasofu (1956). 


D. DISCUSSION 


We have concluded that frequency and probability of occurrence of 
micropulsations vary with local solar time and latitude. Our observations 
show that micropulsations usually are detectable over only a portion of the 
earth’s surface at any given instant, although a few disturbances may be 
world-wide. There may be 0 to 180 degrees phase difference between the 
vertical components of what appears to be the same micropulsation distur- 
bance within a distance of 3000 miles. It has been shown that the distribution 
of directions of the perturbing vector is different at different times or places. 
Neither do the phase relationships between components appear to have the 
same characteristics everywhere. 
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The distinction between Pc and Pt remains, however, and the careful 
separation of these two types of micropulsation may hold the key to their 
origin. We also hope that the study of the higher-frequency disturbances 
(up to perhaps 10 c.p.s.) which usually accompany magnetic storms will 
lead to a clearer picture of the initiation of ordinary micropulsations. The 
fact remains that none of the present theories of the origin of micropulsations 
predict the finer details of their behavior. 
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